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SUMMARY 


This final report describes the physical models em- 
ployed in the NASCAP ( NASA Charging Analyzer Program) code, 
and presents several test cases. A NASCAP User's Manual is 
available under separate cover. 

NASCAP dynamically simulates the charging of an object 
made of conducting segments which may be entirely or partially 
covered with thin dielectric films. The object may be subject 
to either ground test (electron gun) or space (magnetospheric) 
user- specified environments. The simulation alternately 
treats (1) the tendency of materials to accumulate and emit 
charge when subject to plasma environment, and (2) the 
consequent response of the charged particle environment to 
an object's electrostatic field. 

NASCAP is applicable when the Debye length in the 
plasma environment is long compared with body dimensions. 

Then particle trajectories are determined primarily by the 
charge on the satellite and only trajectories which begin or 
end on the object need be considered. For those cases in 
which photosheath effects are important, a first order ex- 
plicit sheath calculation is provided. 

NASCAP contains an object definition language which 
facilitates construction by the user of complex objects built 
of the basic cube, wedge, and tetrahedron elements, and 
allows specification of surface materials. The object 
definition output provides the information required by the 
conjugate gradient potential solver and by the various 
charging and emission routines. 

Parameterized formulations of the emission properties 
of materials subject to bombardment by electrons, protons, 
and sunlight are presented. Values of the parameters are 



suggested for clean aluminum, A^O^, clean magnesium, MgO, 
SiOj, kapton, and teflon. A discussion of conductivity in 
thin dielectrics subject to radiation and high fields is 
given, together with a sample calculation. 

Results of test cases run with the NASCAP code are 
presented. The test cases include bare aluminum and 
dielectric-coated plates under test tank conditions, alumi- 
num and dielectric-coated spheres under space conditions, 
and the SSPM (aluminum plate with four material samples) 
under both ground test and space conditions. 
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1. INTRODUCTION 


This final report describes the work performed during 
the past year at Systems, Science and Software under Contract 
NAS3-20119 to study the electrostatic charging in materials. 

The end result of our study is a computer code, NASCAP ( NASA 
Charging Analyzer Program) which can simulate three dimen- 
sionally the dynamical charging of complex objects in either 
laboratory or space environments. This report describes 
the physical models employed in the code and presents the 
results of the first NASCAP calculations. A NASCAP User* s 
Manual , which fully describes the operational details of the 
code, is available under separate cover. 

NASCAP is able to predict how an object made of conduct- 
ing sections which may be entirely or partially covered with 
thin dielectric films responds to a specified charged-particle 
environment. The environments of interest are those found in 
the earth's maghetosphere and in a ground-based test chamber 
designed to simulate spacecraft charging effects. NASCAP 's 
approach is to divide the spacecraft charging problem into two 
sections: (1) the tendency of materials to accumulate and 

emit charge when subject to plasma environment, and (2) the 
consequent response of the charged particle environment to 
an object's electrostatic field. NASCAP treats both these 
sections in sufficient detail to simulate the charging of a 
complex satellite. 

The objective of the inaterials study portion of our 
program has not been to break new ground in understanding 
fundamental material properties, but to review the existing 
literature and determine which processes are the most impor- 
tant ones to consider for the charging analyzer program. 

Much study has gone into the selection of which material pro- 
cesses are important so that the accuracy of any calculation 
would not be impaired due to the neglect of a dominant 



mechanism. This study has also exposed voids in literature, 
where important relevant properties have not been adequately 
measured, and where there are no good theoretical values for 
necessary material parameters. In spite of such voids, we 
believe that the material properties are sufficient for NASCAP 
to be a useful dynamical charging model. 

The response of the charged particle environment to an 
object's electrostatic field requires the calculation of 
the electric potentials on and near complex objects, and the 
determination of how those potentials influence charged 
particle trajectories. The electrostatic potential about the 
satellite or in the test tank is calculated by NASCAP using a 
finite element formulation of Poisson's equation. Under 
magnetospheric conditions, the Debye length X D « (kT/4Trne 2 ) ly/2 
is typically hundreds of meters, so that space charge can be 
ignored, except for a positively charged satellite which may 
develop a photoelectron sheath. The computational space 
consists of an arbitrarily large number of nested cubic 
meshes. The resulting set of several times 10^ linear equa- 
tions is solved using the Conjugate Gradient technique. The 
satellite or test object is defined within the innermost mesh 
and may have surfaces normal to any of the twenty-six cubic 
symmetry directions. It consists of one or more conductors 
which may be covered with thin dielectric layers. The con- 
ductors may be floating, held at fixed potentials, or biased 
relative to one another. 

The net charge accumulation by each surface cell of the 
satellite is calculated in the presence of the electrostatic 
and magnetostatic fields about the satellite and specified 
environmental characteristics. In the ground test case, the 
incident flux is provided by a monoenergetic electron gun of 
specified beam profile. In the space case, the incident flux 
of electrons and ions at surfaces is determined using the 
reverse trajectory sampling method. The ambient plasma may 
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be isotropic and Maxwellian or may be represented by any of 
several sets of data from ATS-5 prepared by MAYA Development 
Corporation for S^. Alternatively, a spherical probe approxi- 
mation may be used. Optionally, a first-order photo-sheath 
calculation may be performed, but typically, because magneto- 
spheric Debye lengths are large compared with spacecraft 
dimensions, space charge is neglected. 

Section 2 of this report describes the potential and 
flux models used in NASCAP. The material property descrip- 
tions are discussed in Section 3 and Appendix E, while test 
case results are described in Section 4. Appendix A contains 
a description of the ambient space environment. An experimental 
plan for NASCAP verification is contained in Appendix Bt, Ap- 
pendices C and D contain details of the finite element poten- 
tial solver. Finally, a paper describing the charging of a 
materially complex spacecraft using NASCAP comprises Appendix 
F. This paper was presented at the IEEE conference on SPACE/ 
RADIATION Physics held July 14-17, 1977, in Williamsburg, 
Virginia. 



2. NASCAP PLASMA MODEL 


2.1 INTRODUCTION TO NASCAP 

Because it is bombarded by ions and electrons, a satel- 
lite will accumulate, emit and redistribute charge, as well as 
undergo material degradation. In addition, the flux of par- 
ticles to the satellite will be substantially influenced by 
the satellite's own electromagnetic field. To describe this 
process, we have developed a three-dimensional dynamical com- 
puter code, NASA Charging Analyzer Program. This computer 
code, NASCAP, simulates the charging of geometrically, materi- 
ally, and electrically complex objects in both ground test and 
magneto spheric environments. A block diagram of the code is 
shown in Figure 2.1. In this section, we discuss the physical 
models currently in the NASCAP code. More details of the 
NASCAP model can be found in Reference 1. 

A summary of typical length and time scales involved 

is given in Table 2.1. A quick calculation shows that a one 

meter satellite subject to the full incident current of 

10 amperes/m will charge at ^10 volts/sec. If we suppose 

the satellite sufficiently near steady state that the net 

charging current is at least two orders of magnitude below 

the incident current and that ten volts per timestep will give 

sufficient accuracy, then the timestep for our computation will 
—3 

be ^10 seconds. During the differential charging process, a 
satellite is even closer to steady state. In such circumstances, 
NASCAP has been shown able to take timesteps of 100 seconds or 
longer . 

A timestep of one millisecond or longer is long compared 
with any time characterizing the plasma, which can therefore be 
treated by the time-independent Vlasov equation. This allows 
us to calculate the flux to the satellite by the reverse tra- 
jectory technique described below. The electrostatic potential 
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Table 2.1. Orders of Magnitude Characterizing Satellite 
Charging in the Magnetosphere 


Object dimensions (L) 

Transit times: 

Light 

Hot electrons 
Hot ions 

Plasma frequency (w ) 

Sr 

Collision frequencies 

Incident current 

Debye length 

Larmor radius 
Hot electrons 

Hot ions 

Time for 

Circuit element breakdown (arcing) 

Charging of bare conductor 

Differential charging of thin 
dielectric overlayer 

Charge redistribution in a 
dielectric 

Change in environmental conditions 
Material degradation 


1-10 meters 

.O 

10 seconds 

-7 

10 seconds 
10” 3 seconds 
10 5 -10 6 sec” 1 

<< 0 ) 

p 

-5 2 

10 amperes/m 

10 2 -10 3 meters 
4 

10 meters 
, -6 

10 meters 

10 ® seconds 
10 3 seconds 

10^-10 3 seconds 

£ 10^ seconds 
10^-10® seconds 

7 

> 10 seconds 
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is governed by Poisson's equation. Since the timestep is long 

compared with charge redistribution times associated with 

metallic conduction, conducting portions of the satellite 

may be treated as equipotentials. Differential charging, 

however, is governed by the large capacitances and low conducti- 

-3 +3 

vities of dielectric films, so that a timescale of 10 - 10 

seconds is suitable for its study. 

While an entire satellite has a dimension of meters, 
adequate representation of its geometrical and compositional 
complexity requires a spatial resolution of 10 cm or better. 

4 

Even this coarse resolution requires ^lO mesh points in the 
immediate neighborhood of the satellite alone. NASCAP there- 
fore uses a finite element, nested-mesh technique in order 
to achieve an accurate potential near the satellite while 
covering a substantial region of space with appropriately 
reduced resolution. NASCAP generates automatically a series 
of grids within grids as shown in Figure 2.2. All of the 
object is constrained to be within the innermost grid. 

Objects are constructed from the four basic building 
blocks shown in Figure 2.3. These can be assembled in almost 
arbitrary fashion to generate the variety of objects in 
Figure 2.4. By adding these octagons, wedges, cubes, etc. 
together, one can construct an object as complicated as the 
primitive representation of ATS-6 shown in Figure 2.5. Some 
limitations of the present object definition routines are 
apparent in this figure. First, no object can be treated as a 
thin sheet. As a result, the solar panels, for example, are 
modeled to be a full cell in thickness. Secondly, booms must 
also be a cell in width. 

The surface of the object may be covered with a thin 
layer of dielectric material, or may be exposed metal. The 
code automatically handles the electrical coupling between 
dielectric surfaces and underlying conductors. An example of 
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Figure 2.2. Cross-section of grid, showing first four embedded meshes. 




Figure 2.3. Four shapes of volume cells to be considered by 
NASCAP code: (a) empty cube; (b) wedge-shaped 

cell with 110 surface; (c) tetrahedron with 111 
surface; (d) truncated cube with 111 surface. 






Figure 2.4. Six objects constructed by the NASCAP code. The top 
figure is drawn without hidden line elimination, and 
the lower is a full shadowing treatment showing indi- 
vidual surface cells. 
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13 


Figure 2.5. A model for a geometrically complex object resembling ATS-6. This model 
illustrates the geometrical capabilities of the current NASCAP code. 







the variety of material surfaces possible is shown in Figure 

2 . 6 . 


NASCAP is applicable when the Debye length in the plasma 
environment, A Q , is long compared to body dimensions. In these 
circumstances, particle trajectories and the associated charg- 
ing currents are determined primarily by the charge on the 
satellite and are nearly independent of the distribution of 
space charge around the satellite. In the absence of space 
charge effects, only those particle trajectories which begin 
or end on the object relate to spacecraft charging. Moreover, 
invariance under time reversal and the conservation of phase 
space volume along particle trajectories greatly simplifies the 
determination of the current of particles intercepted at object 
surfaces. Finally, when object dimensions, L, are small com- 
pared to Ap, asymptotic potentials of the form Q/r, where Q 
is the object charge, are attained at distances r < A^; conse- 
quently, the computational mesh required for numerical solu- 
tions is established by object dimensions rather them the much 
larger Debye length. 

From the computational point of view, the long Debye 
length approximation reduces both storage requirements and 
the expense of the iterative process that would be required in 
the self-consistent determination of the potential and space 
charge distributions. Thus, there is a substantial economic 
incentive for NASCAP 's utilizing the long Debye length ap- 
proximation where the physical circumstances justify it. 

The neglect of the effect of space charge on the poten- 
tial distribution is valid in those circumstances in which the 
satellite is differentially charged to large negative potentials 
in the range of a few hundred to a few ten thousands of volts, 
such as can develop during magnetic substorms. However, secon- 
dary or photoelectrons can form sheaths with dimensions in the 
range from several centimeters to a few meters near surfaces. 
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Figure 2.6. Surface material composition of a model satellite. 
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When the surface is at a large negative potential, the emit- 
ted electrons will be swept away; for positively charged 
surfaces, or for negative potentials in the range of a few 
volts or tens of volts, however, the space charge in such 
sheaths will strongly influence potential fields around the 
object. 

The charges in the space surrounding a body consist 
not only of ambient particles, but also of charged particles 
emitted from surfaces as the result of impact of electrons, 
photons and positive ions, and in some cases, of electrons 
emitted prior to and during the process of dielectric break- 
down. Emission processes and the trajectories of emitted 
particles play a prominent role in establishing the levels 
of charge and potential on both dielectric and conducting 
surfaces. Particle emission is taken into account by a 
first-order explicit sheath routine. 

The equations that describe the plasma in the 
neighborhood of the satellite are the equilibrium Poisson- 
Vlasov equations: 


V^<{> = -e J’t d v/e Q 

NASCAP uses an implicit equilibrium finite difference 
analogue to these equations: 

V 2 ^ = -e d 3 v/e Q 
v* vf t+1 - | = o . 

in ^ 

3v 


( 2 . 1 ) 

( 2 . 2 ) 


( 2 . 3 ) 

( 2 . 4 ) 
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The time dependence occurs in the boundary condition. In 
particular, the potential (Equation (2.3)) reflects charge 
buildup on the satellite in response to currents from the 
external plasma. 

In the following sections of this chapter are pre- 
sented the techniques used in NASCAP to solve the potential 
Equation (2.3) and the Vlasov Equation (2.4). In the final 
section, we briefly outline the time sequence of the code. 
Further details of the code can be found in the NASCAP User's 
Manual. ^ 

2.2 POTENTIAL CALCULATION 

In calculating the potential in three dimensions around 
an arbitrary object, a gridded method must be employed since 
the specification of the surface is far too general for ana- 
lytical or multipole techniques. Since satellites are the 
order of meters in length, we need at least 10 cm resolution 
as an upper bound in the vicinity of the spacecraft. However, 
for determining particle orbits, the fields hundreds of meters 
away must also be known. In order to keep storage down to a 
reasonable level, some type of variable gridding must be em- 
ployed. This precludes the use of any straightforward Fourier 
transform technique. One technique for achieving high reso- 
lution in the region around the object and still being able to 
handle vast quantities of space is through local mesh refine- 
ment. Finite difference approaches, however, have difficulty 
in mesh transition regions, especially when grid lines are 
terminated and generally lose an order of accuracy in such 
regions. 

As a result of this, NASCAP employs a finite element 
approach using right parallelepiped elements and blended 
linear univariate edge interpolates. This permits the same 
degree of accuracy over the entire mesh, even though the mesh 



elements differ in size. It results in the standard trilinear 
interpolation scheme for each element. 

The fundamental approach is to solve Poisson's equa- 
tion 


eV 2 <p = -p/e Q (2.5) 

by solving the associated variational principle 

0 = §?(/ dV [f (V(p)2 + f^]J ♦/(♦ V(J>-dS + dsj. (2.6) 

The first term in the integrand corresponds to the Laplacian 
operator. The second term is the volume spacecharge contribu- 
tion. The remaining terms are surface contributions, refer- 
ring to the surface charge and electric field, respectively 
(see Appendix C) . 

In the variational calculation, we use locally defined 
basis sets, that is, trilinear interpolants within each cube- 
like element. Finite mesh volumes are given the correct vari- 
ational weight, ensuring the maintenance of accuracy through 
mesh transition regions. The problem of local mesh refine- 
ment is approached by having grids within grids, that is, a 
Chinese doll-like hierarchy of grids shown schematically in 
Figure 2.2. The theory of this technique is discussed in 
References 2 and 3. In order to have high computational 
speed, the linear equations resulting from the variational 
principle [Equation (2.6)] in the interface region were 
coded up explicitly in a series of thirteen subroutines. 

These same routines are used for interfacing any pair of the 
meshes. 


The NASCAP code does not require that objects be com- 
posed solely of rectangular parallelepipeds; it allows surfaces 
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normal to any cubic symmetry direction, i.e., 110 and 111 as 
well as 100. This requires the treatment of surface elements 
shaped as right isosceles triangles, equilateral triangles, 
and /2 * 1 rectangles, as well as squares, and three new 
shapes of volume elements in addition to the cube (0<x<l, 

0<y<l, 0<z<l) (see Figures 2. 3-2. 4). 

(1) l<x+y<2 , 0<z<l 

(2) l<x+y+z<3 

(3) 2<x+y+z<3 

To meet the requirement that the potential be continuous, 
we have adopted the convention that the potential be bilinear 
on a square surface element, and linear on a triangular surface 
element or face. This results in a proliferation of "special 
cells" which have one or more square faces divided into two 
right triangles by the presence of surfaces in neighboring 
cells. Nonetheless, for any cell, we can write 

|V<0| 2 dv = 2 W i;j <J>. (2.7) 

cell ^ 

where i and j index the vertices of the cell. The coef- 
ficients are derived by "linear blending" techniques 

(see Appendix D) . Thus, any "special cell" can be fit into 
the finite element scheme. NASCAP provides for a sufficient 
number of "special cells" to give the user reasonable flexibil- 
ity in objeet specification. 

A general derivation of the linear equations resulting 
from application of the finite element formalism to a system 
of charged conductors is given in Appendix C. In this sec- 
tion, we specialize this method to a cubic mesh with plane 
surfaces composed of thin dielectric layers. 
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NASCAP requires an object constructed so that mesh 
points lie on its dielectric surface. For each surface ele- 
ment, define nodes on the conductor immediately below each 
surface node. Taking the potential within the dielectric 
element to be a linear interpolant of the potentials of its 
eight (for a rectangular element, six for a triangle) vertices, 
(C.25) becomes: 


WPCOND (I ,n ) 
c 



( 2 . 8 ) 


where I is a surface point, the sum runs over those surface 
elements overlaying conductor n of which I is a vertex, 

w 

and A^, e^, d^ and n^ are, respectively, the area, dielec- 
tric constant, thickness, and number of surface points associ- 
ated with the k*"* 1 surface element. Terms proportional to the 
dielectric thickness are ignored. Equations (C.22a) can now 
be written 

- WPCOND (J,n ) <t> + WPCOND (J ,n )<J>(J) 

L J C J n c J c 


* P(n c ) + Q(n c ) (2.9) 


where Q(n ) is the charge on conductor n and 

C C 


p(n o> - S \ ° k ^ 


( 2 . 10 ) 


where the sum runs over surface elements associated with n , 
and and x^ are the surface charge density and mean depo- 

sition depth. (Similarly, the charge associated with a surface 
point J is 


P(J) 





d 


k 


( 2 . 11 ) 
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( 2 . 12 ) 


Equations (C.22b) are written in the form 

27 

WPCOND(I,n ) <J> + 2 W(J,I) <p r = p(I) 

C n c J=1 J 

where the index J runs over I and the 26 points surrounding 
I in the 100, 110 and 111 symmetry directions. The coeffici- 
ents W are given by 

W(J,I) e k J VN 1 • VN J (2.13) 

^ cell k 

where k runs over all volume cells common to points I and J. 

The coefficients and interpolation functions associated with 
cell types implemented and the linear blending techniques 
used to obtain the weights and interpolation functions are 
described in Appendix D. 

The resultant system of linear equations is solved 
using a conjugate gradient technique. This technique is a 
very efficient method for solving large systems of linear 
equations. It is necessary to use an iterative technique 
because of the large number (tens of thousands) of unknown 
potentials. The implementation of the conjugate gradient 

method in NASCAP is discussed in detail in the NASCAP User's 
Manual . 

2.3 FLUX CALCULATIONS 

The solution of the Vlasov equation is approached using 
the fact that phase space density is constant along a particle 
orbit. The numerical method used for this purpose is the time- 
symmetric, reversed trajectory method which enables us to bring 
the information in the velocity space from the undisturbed 
outer boundary to the spacecraft surface. Since <J>^, the 
value of the potential distribution <J> at the nth stage of the 
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computational process is known, particle trajectories at the 
nth stage can be computed. The incident flux on any given 
surface element is determined by tracking time reversed orbits 
starting from an object surface element and invoking the 
principle of phase space volume invariance along particle tra- 
jectories in a collisionless plasma. If a particle starting 
at r 0 » v 0 on an object surface reaches r,v, then 

f <r o ,v o ) - f(r,v) (2.14) 

where f is the phase space density. Thus, if r = r(r Q ,v o ,t) 
is- a point remote from the object, the orbit is one which a 
particle incident from remote distances can follow to the 
object; therefore, the flux of plasma particles incident on 
the surface can be determined from 

f<V* 0 > * (2.15) 

where f (V) is the distribution function of particles in the 
undisturbed plasma. 

The procedure followed in connecting the spectrum of 
particles in a given element at r Q to the unperturbed plasma 
spectrum is 


1. Track particles with specified velocities v 

from their point of origin r Q along time reversed 
trajectories until they reach the outer boundary 
at r p of the computational mesh, or until it is 
clear that the particles will never reach the 
outer boundary. Here, the outer boundary is 
chosen so that the plasma beyond it is not sig- 
nificantly perturbed by the presence of the 
satellite. 
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2 . Identify the incident particle spectrum at r Q as 

f( r o/ v o ) = f(V) (2.17) 

where V is the velocity with which particles on 
their time reversed trajectories reach the outer 
boundary at r p . 

This process is illustrated in Figure 2.7. If v fc is the 
terminal velocity at the outer boundary due to the backward 
trajectory tracing of a particle emitting with an initial 
velocity v^ at the surface, then 

f(rv.) = f(v t ) (2.17) 

is the distribution function value at the point v^. By scan- 
ning these sampling points throughout the velocity space and 
completing their backward trajectory calculations, we can con- 
struct the whole distribution function f (r ,v ) . 






Figure 2.7. To construct the distribution function f(r o ,v Q ) at 
a spacecraft surface. 
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In the general situation, we have to integrate the dis- 
tribution function to obtain the charging current density for 
each species, that is 

2tt 1 ® 

3 C = n / I I ev^f o (v o ,y,4»)dv o pdyd<|) (2.18) 

J o J o v m (y,«) 

where y = cos0 , v m (y,$) is the minimum escape speed of an 
electron emitted from the surface element with an angle (6,<f>) 
with respect to n, and where f Q (v o ,y,<|>) includes anisotropy 
due to the magnetic .field. 

2 . 4 DYNAMICAL MODEL 

The dynamical model then consists of three parts: 

1. The calculation of surface charging currents 
given potential distributions. 

2. The calculation of dielectric, surface and 
space charge distributions given potentials 
and charging currents. 

3 . The calculation of potentials given charge 
distributions and boundary conditions. 

In the typical one- or two-dimensional equilibrium code, the 
steps are iterated until a self-consistent solution is reached. 
However , the additional complexities introduced by the three- 
dimensional nature of the problem and the sophistication of 
the material properties treatment make iteration for each 
timestep prohibitive. Therefore, the solution sequence is a 
timestepping procedure from one quasi-static state to another. 
Initially, all potentials and charge distributions are speci- 
fied* The dynamical parts of the problem are driven by charge 
accumulation on the body from external sources (ambient, plasma, 
electron gun, etc.), charge depletion (surface emission, etc.) 
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and conduction in dielectrics. Each timestep includes, in an 
explicit fashion, a fully three-dimensional electrostatic 
potential calculation time-staggered with a procedure in which 
incident charged particle fluxes, leakage currents, emission 
currents and emission current induced space charge effects are 
found according to the derived quasi-static equations. 



3 . MATERIAL PROPERTIES 


3.1 GENERAL CONSIDERATIONS 

The interaction of charged particles at the surface of 
a geosynchronous satellite and the redistribution of charges 
in surface layers of the body are determinants of the electric 
stresses in satellite materials. The most important pro- 
cesses for satellite charging are charged particle deposition, 
photoelectron emission, secondary electron production by 
electron impact, electron backscattering and electron pro- 
duction by proton impact. The last process is generally of 
less importance than the others; in eclipse, however, where 
at electrical equilibrium electron and proton currents are 
nearly in balance, protonic and electronic emission processes 
are of comparable importance. NASCAP includes formulations 
for all of these processes. Other processes, such as electron 
emission by the impact of naturally occurring He and 0 ions, 
are not considered by NASCAP, but could readily be included 
if later judged to be important. 

It is difficult to determine the particle emission 
properties of a given material accurately by laboratory ex- 
periments on well-characterized surfaces. For this applica- 
tion, the situation is still worse, since the surface proper- 
ties of satellite materials are not carefully controlled; 
moreover, surface properties can be changed substantially by 
exposure to the magneto spheric environment. These factors 
should be viewed at this time as constituting a fundamental 
limitation on the quantitative accuracy with which electric 
potentials can be predicted. 

Charged particles impacting a surface not only cause 
emission of other charged particles but also deposit charge 
beneath the surface. While charged particle deposition is 
not an important consideration for conducting materials, the 
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depth of deposition and the processes affecting it can have 
a profound effect on the strength and distribution of field 
within dielectric materials. Among the parameters affecting 
the internal electric fields are the intrinsic and field- 
and radiation-enhanced conductivities of the materials. 
Conductivity enhancement may occur by electron irradiation 
in a layer of thickness of the order of the electron range, 
by solar illumination to a depth depending on the optical 
absorption characteristics of the material, and by the pro- 
duction of charge carriers in strong electric fields 
(£l0^ volts/cm) . Accurate determinations of the internal 
fields are rendered difficult not only because of the limited 
amount and quality of relevant data, but also because a com- 
plete and unambiguous theoretical description based upon the 
underlying physical processes is lacking for the materials 
of interest. Even if such a description were available, its 
incorporation into a three-dimensional computational scheme 
could seriously limit the efficiency of that scheme. 

Based on the uncertainties inherent in material 
properties, as well as the high premium associated with the 
efficient operation of a three-dimensional code, a pheno- 
menological approach to the charged particle transport within 
and emission from satellite materials is fully justified. 
Section 3.3 describes NASCAP's approach for the determination 
of electron emission resulting from electron, proton and 
photon bombardment. The approach described is comparable 
to state-of-the-art methods for estimating emission cur- 
rents, and in addition it extends those methods to permit 
emission estimates over extended regimes of energy and angle 
of incidence. The parameters required to estimate emission 
currents from clean aluminum, Al^O^, clean magnesium, MgO, 
SiC^, kapton and teflon are presented. 



We show in the following sections that secondary emis- 
sion due to electron bombardment and electronic charge depo- 
sition are closely related to the range of the incident electron. 
A substantial simplification in the description of deposition 
profiles results if the energy and angle dependence of the 
incident electrons can be described a priori; such a simpli- 
fication permits useful estimates of upper bounds on the 
electric field in the dielectric at its vacuum interface. One 
possibly useful description of the deposition is described in 
Section 3.2. The general validity of such simplified treat- 
ments of charged particle deposition as we consider here can 
best be established by calculations of incident particle spectra 
under a variety of charging conditions. 

For conductors covered by a thin dielectric film, the 
region in which we seek solutions of Poisson's equation in- 
cludes the dielectric. Boundary conditions are applied on 
conducting surfaces, but the charge on the conductor is deter- 
mined not only by plasma currents, but also by leakage currents 
through the dielectric films. Except for extremely thin 
dielectric layers (<0.1 mil), the charge density and potential 
distributions within a dielectric into which charge is injected 
at one surface vary strongly as a function of position. More- 
over, the injected charge is redistributed by charge transport 
processes which are not well understood, and in any case, are 
difficult to quantify. 

We anticipate that in many practically occurring cir- 
cumstances, the potential in the vacuum region depends only 
weakly on the charge distribution within the dielectric layer. 

In Section 3.6, we express the coupling between a conducting 
surface and vacuum through a dielectric film by means of an 
effective boundary condition at the dielectric surface. The 
basic approximation that leads to a computationally simple 
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boundary condition is that the thickness of the layer contain- 
ing the injected charge is small relative to the film thick- 
ness. This approximation leads not only to a simple boundary 
condition, but also to the result that the electric field and 
leakage currents are independent of position throughout most 
of the depth within the dielectric. Under these circumstances, 
it is unnecessary to "zone" the dielectric. Therefore this 
simplified treatment is used by NASCAP. 

While seeking to simplify the treatment of dielectric 
films, it is recognized that there are reasons for knowing 
the potential distribution in the dielectric and for under- 
standing the processes which affect it. First, one should 
understand the limitations on our treatment of dielectrics in 
terms of an effective boundary condition. Second, knowledge 
of the electric field and charge density structure within the 
dielectric is likely to be of importance in the determination 
of conditions marking the onset of breakdown through the film. 

Section 3.7 describes a simple field-dependent bulk 
conduction model for the transport of charge in a dielectric. 

There it is argued that the potential drop across a thin sec- 
tion of dielectric is insensitive to the detailed processes 
which govern charge migration and carrier production within 

the dielectric. This argument is considered in Section 3.8, 
the subject of which is a detailed one-dimensional treatment 
of charge migration that takes into account the production, 
mobility and trapping of a single species of charged carrier. 

Should the conclusion regarding the insensitivity of the 
potential drop have validity beyond the examples that demon- 
strate it, then separation of the calculation of surface potentials 
from the detailed calculation of field distributions within di- 
electric films would be valid. One can then invoke the surface 
potentials calculated in three-dimensions as boundary conditions 
for one-dimensional calculations of internal field distributions 



that incorporate physically reasonable descriptions of the 
kinetics and transport of charged species. 

3.2 CHARGE DEPOSITION 

A good insulator will, with exposure to the space en- 
vironment, develop a charge density profile similar to that 
shown in Figure 3.1. Within a few tens of angstroms from the 
surface a positively charged depletion layer forms due pri- 
marily to emission of secondary electrons and photoelectrons. 
Superimposed on this is a distribution of negative charge due 
to the stopped incident electrons. As indicated in Appen- 
dix F.2, this distribution can adequately be represented by 
a simple exponential, at least for the case of an isotropic, 
Maxwellian plasma: 

N(x) « (1 -A q )/x exp (-x/x) (3.1) 

where A q is the net albedo (see Section 3.4 below), x the 
mean deposition depth, and the energy dependence is based on 
Feldman's range formula. 1 For keV plasma temperatures, x 
is a few hundred angstroms (see Table 3.1). The derivation 
of Equation (3.1) can be easily generalized to a non- 
Maxwellian, non-isotropic plasma. 

Neglecting charge migration due to conduction pro- 
cesses (see section on dielectric properties), and treating 
the depletion layer as a surface charge, we have (q « 
electronic charge) 


p(x,t) * p(x,o) + 


t 

/ dt'js(t f )5(x)lq| 


1-A 


“^|q|f e (t') exp (-x/x ( t * ) ) +|p- Pi (t') 


x(t') 
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(3.2) 
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Table 3.1. Albedo and Range For Electrons in Materials 


Material 

A 

o 

n 

xT** n (A-keV 

Kapton 

0.187 

1.505 

560 

Teflon 

0.235 

1.63 

400 

Si0 2 

0.262 

1.69 

316 

ai 2 o 3 

0.262 

1.69 

219 

MgO 

0.262 

1.69 

233 

A1 

0.299 

1.77 

313 
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for the charge density in the dielectric, where S is total 
secondary yield (due to electron and ion impact and photo- 
emission) , f is the electron flux at the surface, and p. 
is the charge density due to ion deposition. The NASCAP 
potential solver is in principle capable of handling the O*"* 1 
moment, a, and the first moment, 6, of this distribution: 

t 

a(t) » a(o) +| q\f dt 1 (S (t ' ) -ir (1-A ) f (t ' ) +tt (1-A . ) f . (t ' ) ) 

/ o e oi i 

o 


cr(t)6 (t) 


o(o)6 (o) 



dt'[(l-A 0 )x(t')f e (t') 


- (l-A oi )x i (f)f i (f)] . (3.3) 

where we introduce f^, A q ^, and x\ for the flux, albedo and 
mean deposition depth of ions. Because the ionic flux is 
small, an accurate representation of the ionic charge 
deposition profile is not essential. Furthermore, we anti- 
cipate that in most, and perhaps all interesting circum- 
stances, the first moment 6 will have only a negligible ef- 
fect on the external potential. The present version of 
NASCAP ignores the first moment, 6. 


3.3 SECONDARY EMISSION 
3.3.1 Electron Impact 

The emission of low-energy electrons upon electron im- 
pact is one of the most important factors determining the sign 
of a spacecraft's charge. Unfortunately, the pool of experi- 
mental data characterizing this phenomenon is far from ade- 
r 5 — pi 

guate. 1 1 It is therefore necessary to calculate secondary 
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emission using formulae which relate it to better character- 
ized material parameters and which are general enough to 
accept future modification. Furthermore, the formulae must 
be applicable to non-normal incidence, for which the amount 
of available data is miniscule. 

NASCAP uses a formulation based on the range and 
energy loss rate of the incident particles. Nearly all of 
the energy lost by an incident electron goes into electronic 
excitations, and we assume the probability of an electronic 
excitation resulting in an escaped secondary varies exponen- 
tially with depth. We then have 


5 



dE 

as 


-c^xcose 
e dx 


(3.4) 


where 6 is the number of emitted secondaries per primary 
incident at angle 0, and the range and energy loss rate are 
related by 



If the range function is known, the above expression can, in 
principle, be evaluated with the constants c^ and c 2 deter- 
mined from the energy of maximum yield, e^, and the corre- 
sponding yield 6^. 

For a general range expression. Equation (3.4) can be 
evaluated by assuming a constant dE/dx. This is reasonable 
since most of the secondary electrons originate in a thin 
surface layer. The upper limit of Equation (3.4) must be 
set to give the correct total energy loss (set c. = 1 and 
c 2 ® 0) . We then have 
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dx 


6(E,0> 



e 


-C2XCOS0 


o 

[1-exp (-c 2 cos0EdR/dE) ] 
~ C2Cos0dR/dE 

The angle averaged yield then becomes 


6(E) « 2c 1 E (Q-l+exp (-Q) )/Q 2 


(3.6) 


(3.7) 


where 

Q * c 2 EdR/dE . (3.8) 

NASCAP, in fact, evaluates (3.4) by assuming dE/dx 
is linear in x: 


dE 

dx 



(3.9) 


The range is represented by the sum of two exponentials: 


R 


n. 


n. 


r l E 


+ r 2 E 


(3.10) 


For such data as is available this gives a good representation 
of the range for 100 eV < e < 100 keV. (See Figure 3.2.) 

The upper limit, R u , on the integral Equation (3.4) is taken 
as the lesser of the solutions of 


dE 1 
dx* 


R 


u 


0 


(3.11a) 
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(3.11b) 


/ 


u 


Igl ax = E 


Then, letting Q = R u cos we have 


6 (E, 0 ) = c 


p (jH L\ -1 l-exp(-Q) 
u ldE Q ) Q 


n 2 d 2 R / dR \~ 3 1- (Q+l) exp(-Q) 
U dE 2 \**o) Q 2 


(3.12) 


6 (E) = 2c- 


i 4s* 

u\dE_ 


k -l 


^Q-l + exp(-Q)j/Q" 


+ ZR 


2 /d 2 R\/ dR \ ~ 3 

U \dE 2 / 


(3.13a) 


where, in (3.13) Q is evaluated for normal incidence and 


z - f udu iziemy. ysE. <to«> 

-/ Q Z u 


(3.13b) 


which may be expanded for large and small values of Q. 

NASCAP uses this formulation of secondary emission to 

evaluate the constants c^ and Cj from user input parameters 

<5 , E . 
m m 

Figures 3.3 through 3.7 show illustrative curves 
generated by the electron-secondary subroutines. These curves 
are based on the parameters given in Table 3.2* 
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gure 3. 6. Secondary emission coefficients for teflon using Equation (3.9) 



(keV) 


w 


4 J >1 

r 0 

~ Q) 

o c • 

• *H ► 
co a3 ► 
^ -P W 

C 0 M-i 

o o 

-H 0) 

■P P 0) 
Id 0) D 
P 

D 1 <d 
w w > 
0) 

tT' > P 
C P -H 

•h a mh 
cn u 
3 o 
T3 P 
C 0) 

O *C in 
*P ffl r* 
Cu id • 
03 Q rH 
M 

O 

P • P 
O — 

0 } O 

<D in 
> • 
P H 
d 

<U O g 
•H O 

o *a p 

•H *H M~l 
<W »H 
M-I O P 

owe 
o ~ a) 
o 

CN 

c • 

o CO 


03 

P 

C 


c 

o 

a 

x 

<D 


to a> 
WHO 
•P A O 

id 

Em 


s 


c 

id 

p 


O' 

c 


>1 g 
P o 
03 P 
*a M-I -H 
C tP 

o id c 
o p <o 
a> id us 

CO T3 o 


CO 

a> 

p 


•H 




3.3.2 Proton Impact 


Secondary emission dye to proton impact^ is 
treated by NASCAP in similar fashion to the electron case. 

A difference is that the energy loss is well represented by 

if | = cE l/2 /(l + E/ Emax l , (3.14) 

where E max ^ 50 keV. Proton secondary emission can be large 
in the 10-100 keV energy range. , Below 1 keV, the "potential 
emission" process comes into play. However, the potential 
emission coefficient is seldom as large as 0.1 and can probably 
be ignored. Because energetic protons travel long distances 
in straight lines, and because emission by ions is generally 
less important, Equation (3.5) is used in the formulation, 
and the angular dependence is taken to be simply sec 0. 

Secondary emission by aluminum on proton impact is indicated 
in Figure 3.8. 

3.3.3 Energy and Angle Distribution 

The energy distribution of secondary electrons is 
peaked at a few volts. Either a Maxwellian or a uniform 
distribution provides an adequate representation for space- 
craft charging purposes. 

A small emitting surface emits secondary electrons 
into a unit solid angle at 0 at a rate proportional to cos 
0* This results in an isotropic flux of secondaries above 
an extended emitting surface. 

3.4 BACKSCATTERING AND REFLECTION 

3.4.1 Albedo for Electrons 

Backscattering of electrons is discussed in Appendix 
F.l. We describe a large-angle scattering theory similar to 
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Cousinie* 17 ! 



Figure 3.8. Secondary emission by aluminum for proton impact at normal incidence 
Curve is Equation (3.10). Experimental points as indicated. 



MO 1 

that of Everhart, but generalizable to arbitrary angles of 
incidence. (Large-angle -scattering approximations are known 
to be superior to diffusion approximations for low z materials.) 
For normal incidence, and assuming the Rutherford scattering 
cross-section and the Thomson-Widdington slowing down law 
this theory can be integrated to yield 


dE/dx a E , 


n » 1 - <f> 


(3.15) 


where a renormalized exponent a = 0.037Z gives backscattering 
coefficients in good agreement with experiment. This result 
is expected to be valid for 10 keV < E < 100 keV. 

The large-angle scattering theory, together with Monte- 
Carlo (ELTRAN) data and experiments by Darlington and Cosslett, 
indicate that the angular dependence of backscattering is well 
described by 


n^» n(0) exp [r^ (l-cos6) ] 


(3.16) 


where the value of tk is, within the uncertainty in the data, 
what would be obtained by assuming total backscattering at 
glancing incidence, viz. * -log n Q . The net albedo for an 
isotropic flux is then 


= 2 [1 - 


n o (1-log 


n Q ) 3 / (log n Q ) 


(3.17) 


As the energy is decreased below 10 keV the backscattering 

[ 20 ] . .. „ 

indicate an increase of 


increases. Data cited by Shimizu 
about 0.1, almost independent of Z. 
component of backscattering by 


NASCAP approximates this 


6n o = 0.1 exp[-E/5 keV] 


(3.18) 
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At very low energies the backscattering coefficient becomes 
very small and, below 50 eV, backscattering and secondary 
emission are indistinguishable. NASCAP takes account of this 
by a factor of log (E/50 eV) 0(E - 50 eV)/log (20). The 
formula for energy-dependent backscattering, incorporating 
these assumptions, is then 

n Q = {[log(E/O.O5)0(E - 0. 05)0(1. 0-E) /log (20) ] + 0 (e- 1.O)> 
x [0.1 exp (-E/5) + 1 - (2/e) * 037Z ] (3.19) 

where energies are measured in keV. Resultant curves for the 
net albedo are shown in Figure 3.9. 

3.4.2 Energy and Angle Distribution of Backscattered Electrons 

Monte-Carlo (e.g., ELTRAN) data as well as several 
approximate theories indicate that the energy and angle distri- 
butions of backscattered electrons are smooth and exhibit 
surprisingly weak dependence on incident angle. The mean 
energy of backscattered electrons in the 10 keV - 100 keV 
range is about two-thirds the incident energy, while the 
angular distribution is roughly cosQ for all but the most 

glancing angles. Accordingly, the energy and angle distribution 
of backscattered electrons may be approximated as 

00 

f _ (E , 0 ) * 2K I dE.f (E. ) (E/E?)cos0 (3.20) 

o O / X X X 

E 

where f(E^) is the angle-averaged incident flux at energy 

3.4.3 Reflection of Protons ^^'^^ 

There is little data on reflection of protons from solids. 
Indications are that the net reflection from low-z materials is 
no more than 10 percent. However, many protons are reflected 
as neutrals or negative ions, so that the charge reflection 
coefficient is surely small. NASCAP neglects reflection of 
protons . 
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3.4.4 Sputtering 
. T121 

Kaminsky indicates that sputtering coefficients for 
metals under proton bombardment seldom exceed 0.1 atoms/ion. 
Furthermore, the charge distribution of sputtered particles is 
unknown. In view of the low proton fluxes, sputtering should 
be a negligible factor for spacecraft charging, except insofar 
as it results in surface degradation. NASCAP does not treat 
sputtering. 

3.5 PHOTOEMISSION 

In a sense, photoemission is a relatively simple process 
to treat, since the spectrum of incident particles is unique 
and fairly well characterized, viz. the solar spectrum. ^ 
However, few measurements of photoelectric yield or optical 
properties have been made in the vacuum ultraviolet. Further- 
more, the photoyield is strongly dependent on surface condition. 

NASCAP expects the user to enter a number (based, say, 

on the work of Feuerbacher and Fitton^^) for the photoyield 
2 

in amps/m under normally incident sunlight. The program will 
then correct for angle by assuming a constant yield per incident 
photon . 

If more data were available, the yield could be calcu- 
lated by: 


Y(0) 


r, 


de y (e)f (e)h(e,0)cos0 


(3.21) 


where y(e) is the yield (el/photon for normally incident mono- 

2 

chromatic light at energy e,f(e) is the solar flux (photons/cm - 
sec), and h(e,0) is the ratio of the monochromatic yield at 
angle 0 to the normally incident yield. As discussed in 
Appendix F.3, h(e,0) depends only on the optical properties of 
the material, and is given by 



(3.22) 


h(e,6) 


sec 6 


|E(0) I 2 .. l+2a (o) L 
| E (o) j 2 1+20 (B)L 


where E(0) is the electric field just inside the solid calcu- 
lated using the usual electromagnetic boundary conditions, 
a(0) is the E-field attenuation coefficient normal to the 
surface, and L is the escape depth for photoelectrons, which 
may be taken as the inverse of the constant c 2 appearing in the 
preceding discussion of secondary emission. 

3.6 EFFECTIVE BOUNDARY CONDITIONS IN THIN DIELECTRIC MATERIALS 

For conductors covered by a thin dielectric film, it is 
convenient to express the potential drop across the film by 
means of an effective boundary condition at the dielectric 
vacuum interface. The desired relation between the potential, 

<f> Q , at the vacuum-dielectric interface (x = 0) and the potential, 
♦ , at the dielectric-metal interface (x = d) follows from the 
definition 

d 

♦ c “ - ~f E - (x)dx (3.23) 

o 

and the boundary condition 

s 

k(o)E + (o) - E_(o) = (3.24) 

o 

by integrating Poisson's equation 


SS (kE) = f- (3.25) 

o 

through the dielectric material. Here, 

E o * 3?? x 10 " 7 farad /">eter , ( 3 . 26 ) 
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E + (o) and E_ (o) are the electric fields at x = 0 + and x = 0~, 
respectively, « (x) is the relative static dielectric constant 
at the position x, p is the charge density, and a S is the 
density of surface charge. For a single dielectric layer, the 
effective boundary condition takes the form 


n*V<{> - < 


(<f> -*• ) 


- ±- <a 5 + a v (l - J)> 


(3.27) 


where 


CJ 


V 



p (x)dx 




X 



xp (x) dx 


(3.28) 


(3.29) 


and n is the unit normal vector directed from the dielectric 
into vacuum. It follows that if x << d, that is if the excess 
charge injected into the dielectric remains near the vacuum - 
dielectric interface, it is a good approximation to consider 
the net injected charge as a surface charge. 


cr 


eff 

s 


_s 


+ a 

V 


This conclusion, however, applies only for the purpose of compu- 
tation of potentials in the space outside the dielectric. 

Clearly, the determination of the electric field within the 
dielectric, particularly near the vacuum— dielectric boundary, 
requires a knowledge of the charge distribution in the dielectric. 
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3.7 SIMPLE MODEL FOR LEAKAGE CURRENTS AND FIELDS IN THIN 

DIELECTRIC FILMS 

The mathematical treatment of charge leakage through 
an ohmic medium is extremely simple. Consider a beam of 
electrons incident onto one face of a thin dielectric slab, 
and denote by the current of beam electrons within the 
medium. The charge deposition profile is given by q^j^/Sx), 
where q is the electronic charge. If p is the charge density, 

E the electric field, a the constant conductivity of the 
medium and e the dielectric constant, the equations governing 
the electrical behavior within the dielectric are 

3 E 

e = P (3.30a) 



(3.30b) 


(3.31) 


For a deposition profile of time invariant shape 

j fa (x,t) = j£(t)f(x) (3.32) 

the spatial dependence of the charge distribution of an ini- 
tially uncharged dielectric is also time invariant. Moreover, 
if the range of electrons in the dielectric medium is less 
than its thickness, the electric field in the region beyond 
the range of the incident electrons is independent of posi- 
tion. Typically, the thickness of the layer of dielectric 
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adjacent to the "vacuum" is greater than 1 mil (2.5 x 10~ 3 cm) 
whereas the range of 10 keV electrons in teflon is only about 

_4 

10 cm, so that the electric field varies spatially only in 
a small region near the surface at which electrons are inci- 
dent. In such circumstances, one can suppose that the 
deposited charge resides on the surface so that the electric 
field is spatially uniform through the dielectric layer. The 
approximation would be a good one for the calculation of the 
electric fields everywhere except within the deposition zone. 

As long as deposited electrons are treated as a sur- 
face charge, the electric field is spatially uniform through 
the dielectric even for a field dependent conductivity. If 

the conduction current is j = j (E) , then 

c c 


3p „ 3j c _ 3j c 3E = p 3j c 

5t F 5T5x sJT 


(3.33) 


Thus if p vanishes initially it vanishes for all times and 
3E/3x =0. On the other hand, within the deposition zone 
the form of the charge density profile would be modified by 
the effects of conduction. 

As long as we maintain the assumption that any excess 
charge in the dielectric remains in a surface layer, the 
problem of calculating dielectric leakage currents is mathe- 
matically straightforward. It still remains, however, to 
describe the manner of dependence of conduction currents on 
electric field strength. 

Many authors have advanced models for electrical con- 
duction by dielectrics at high fields based on the classical 

r 95 1 r 2 6 1 

ideas of Schottky 1 J and Frenkel. 1 1 These models have 

f 271 

been summarized by Adamec and Calderwood. The latter 

authors have also proposed a model for polymeric insulating 
materials which yields a relationship between conductivity 


53 



fisld strength in good Agreement with experimental values 
for a number of polymers (including polyamide at 250°C) over 
a range of field strengths ranging from 10 4 to 10 6 volts/cm. 
According to their model, the field dependence of conductivity 
may be expressed as 


o q (T) 


(2 + cosh (S F E 1/2 /2 kT)) 


(3.34) 


fields less than about 10 volts/nu Here/ in ntks units f 



(3.35) 


is the Frenkel parameter, < the relative dielectric constant 

at high frequency, k is Boltzmann's constant, T the absolute 

temperature, e * 1/36 x 10" 9 farad/meter and a is the 
u o 

intrinsic low field conductivity, which may depend on the 

temperature . 

The precise values of parameters to be used in Equation 
(3.34) are not at all certain. Nominal values of resistivity 
have been given for kapton, ^ 28 1 teflon* 29 ^ and fused silica; ^ 30] 
presumably, the measured nominal values of resistivity are 
determined from a resistance which is obtained by dividing 
an applied voltage by the observed current through the sample. 
The values of resistivity so obtained depend, in general, on 
the applied voltage, sample thickness and temperature, and 
quite possibly the measurements also reflect electrode and 
space charge effects. There is also some uncertainty con- 
cerning the behavior of the conductivity at high fields. 

Equation (3.34) gives a high field conductivity which varies 
as exp ( 8 f E 1//2 /2 kT) . Models proposed by Johnscher ^ 33 ^ and by 
Mead (32] give field conductivities varying as exp (8„E 1/,2 /kT) , 

IT 
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that is, with the exponents having an argument twice as 
large as that in Equation (3.34). 

Let us consider the magnitude of field enhanced con- 
ductivity. For a material with k * 3.5 at T * 300°K, the 
quantity 8p/2kT has the numerical value 

8 F -4 

2kT ~ x ^ 

when fields are given in volts/meter, then 
a * F a o ex P ( 7 * 8 x 10 ” 4 E ly/2 ) 

q 

z 407 a Q for E * 10 volts/meter. 

Such enhancement of the conductivity in kapton by fields in 
the megavolt/cm range are consistent with the measurements 
of Hoffmaster and Sellen. For kapton, the nominal value 

of o q at T = 300°K is of the order of 10~^ - 10”^® mho/m. 

g 

At fields of order 10 volts/meter, the dielectric is 
able to support steady state currents in the range from a few 
hundredths to a few tenths of nanoamps/cm . It is also 
interesting to observe that a dielectric relaxation time 
defined by 



is reduced from days to minutes in going from low fields to 
fields of one megavolt/cm. In Section 3.8, we will include 
field enhanced conductivity in the determination of the 
field structure in the interior of a dielectric. 



3.8 KINETIC DESCRIPTION OF PROCESSES IN A DIELECTRIC 

Numerous experiments over the years* 34 ' 35 ^ have demon- 
strated that polymers such as teflon conduct small, but sig- 
nificant currents with applied fields of the order of 10 5 v/cm, 
and are capable of storing charge for long periods of time 
(days at low temperatures). Moreover, for high applied fields, 

>10 5 v/cm, currents increase in a roughly linear manner on a 
~ 1/2 

Schottky (log I versus V ' ) plot. These facts indicate the 
presence of a species of "free" carrier with some mobility, 

U, the presence of deep traps, and probably the existence of 
field assisted excitation of charge from trapped sites. Be- 
low, we construct one simple model which incorporates these 
important features of insulating materials used in space 
applications. For simplicity, we assume that electrons are 
the only carriers of electricity. 

Our purpose here is to examine properties that are 
required to model the charge migration and field buildup in 
a dielectric. Such a description will allow 

1. A better understanding of the material para- 
meters influencing charge transport. 

2. A basis for assessing the approximations of 
the simple phenomenological description of 
leakage currents given in the preceding 
section. 

For the present, we consider bulk effects in the absence of 
radiation-induced effects. 

The temporal evolution of the charge and field distri- 
bution in a dielectric is governed by Equation (3.30a) and 

§£ = - y *3 (3.36) 


i 
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where p is the charge density and j the current density. To 
complete the description, we must relate the current density 
j to p and E. 


The particle current density j (and electric current 
density j ) is the sum of the plasma electron current pene- 
trating into the dielectric, j^, and a drift current 
j c * -n c ME 

j = qj p - q(j b + j c > • (3.37) 

Here n is the density of free carriers and y their mobility. 

In addition, the capture and release of charge from traps is 
governed by^^ 

3n c 

j£- = -n c <ov> [N fc - n t J + vn fc - V*j c (3.38) 

where <av> is the capture rate, v is the trap release fre- 
quency, is the trap density, and n^ is the density of 

occupied traps. For low fields, release from traps is 
thermally activated and 

Ae 

0 

v “ v e (3.39) 

o 

where Ae is the depth of the trap measured from the "conduc- 
tion band" and 6 the temperature in energy units. For the 
trap population, we have 

3n t 

= n c <av> [N t - n fc ] - vn fc . (3.40) 

The charge density given by 

P - q(n c + - n co - n to > . < 3 - 41 > 
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where n cQ and n fcQ are the thermal equilibrium values of free 
and trapped carrier density, satisfies Equation (3.36). 

Typically, n fc in space insulating materials is sub- 
stantially less than N fc , at least before the onset of di- 
electric breakdown. The density of free carriers is small 
compared with the trapped charge density, both in thermal 
® < I u ilikriura and for a dielectric containing a space charge, 
and the free time, t ■ [N t <av>]~ 1 , is short (<10” 9 sec) com- 
pared with the time scale of macroscopic variation. Then it 
is a good approximation to set 3n /3t =0. If additionally, 
V*j is neglected, one obtains a component of current of the 
form py eff E, but with a mobility y gff j* y. Although neglect- 
ing V*j c in Equation (3.38) may not be. a valid approximation 
for all conditions of interest, it is nonetheless useful to 
examine its consequences. 

We find 


n 

n 

C Z \)T = 

CO 

n t 

n to 

p 5 q(n t 

o 

c 

1 


(3.42) 


(3.43) 


j ■ qjp = qn c yE = qn°yE + p(yvt)E 

= a Q E + py eff E . (3.44) 


Thus, transport of excess charge is proportional to the net 
charge density, but with a trap modulated mobility 

U eff = (UVT) « y , (3.45) 
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a quantity frequently appearing in the literature on charge 
storage , and the object of many experimental observations. 

[371 

Gross, et al., have used Equation (3.44) in the analysis 
of experiments on the transport of charge injected into a 
dielectric. In the context of modeling the electrical be- 
havior of insulators in the space environment, it is 
important that we understand the significance of the ex- 
perimentally observed mobilities. 

The preceding theory also permits treating field en- 
hanced thermal activation of occupied traps. Assuming that 
trapped sites are neutral when occupied and invoking a Poole- 
Frenkel ionization mechanism, ^®^^e required modification 
consists in replacing v q by v o e aEl where the coefficient 
a is related to the dielectric constant and temperature of 
the medium. Provided that space charge effects are not 
pronounced, this leads to the frequently observed linear 
relationship between Jln(J) and E 1//2 at high fields. 

The theory, as elaborated so far, does not admit the 
effects of conductivity induced by electron or solar irradia- 
tion. Typically, under constant irradiation with a dose rate 
D, the material acquires a steady state conductivity follow- 
ing Fowler's law^ 3 ^ 

a - k(D/D Q ) A (3.46) 

where D Q is a reference dose rate, usually taken as 1 rad/sec, 
and A is a material parameter with values between 0.5 and 1 
depending on the energy level distribution of traps in the 
material. A simple model with a single trapping level is not 
expected to be a good model for radiation-induced conductivity 
(RIC) . Nevertheless, RIC can be simulated by adding to Equa- 
tion (3.38) and subtracting from Equation (3.40) a term Rn^, 
where R determines the rate of trap ionization by the impressed 
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radiation field. The effect of RIC on field structure within 
the dielectric is determined in this manner in Section 3.10. 

To place the foregoing considerations into the context 
of the entire spacecraft charging model, it is important to 
have some notion of the relative time scales involved in the 
problem. The important characteristic times for the dielec- 
tric are t and 1/v, as previously defined, the dielectric re- 
laxation time 


t 


d 



(3.47) 


where a is the conductivity, and the transit time t fcr across 
the thickness of dielectric. These times should be compared 
with body charging time t ; for an initially uncharged body 
of dimension R, the latter time scale is estimated by 

4TTR 2 jt c = C0 (3.48) 

where j is the one-sided plasma electron current density, 

© is the plasma temperature in volts and C the capacitance 
relative to the zero of potential. For R of order one meter, 

(C ^ lO” 1 ^ farad), j 10*" 5 amps/m 2 , 0 'v 10^ volts, 

-3 

t =10 sec . 
c 


Dielectric relaxation times for good insulators are much 
longer than this, even for levels of conductivity that may be 
induced by the radiation levels at geosynchronous altitudes. 
The trap residence time is highly variable, depending on 
temperature and field strength, and can be larger or smaller 
than t . Transit times are very long, 

C 


*tr 


= L/p 


eff 


60 



_2 

for a thickness L *\» 10 cm, a trap modulated mobility 

**1 0 2 * 

U eff ^10 cm /volt-sec and a field strength E ^ 10 5 v/cm. 


t tr » 10 3 sec . 


The carrier free time T f is much shorter than t c , so that 

the free carrier concentration relaxes instantaneously to a 

quasi-steady value 8n /at z 0. 

o 

The important consequence of the foregoing considera- 
tions is that while the body charges to a quasi-equilibrium 
characterized by a vanishing net current in the space around 
the body, charges in the dielectric hardly move at all. This 
suggests that the problem of spacecraft charging separates 
into two rather distinct parts; one being the overall charging 
equilibrium, followed by redistribution of charges on the 
body. 

Before attempting numerical calculations based on the 
carrier kinetics, it is worthwhile to relate the kinetic de- 
scription to the simple model described in Section 3.7. First 
we observe that at low fields and in the absence of buildup of 
excess charge, n c and n^_ have their thermal equilibrium values 
n co and n to' res Pectively, which are related by 


n co <av> [H t - n tQ ] 


vn 


to 


with a dark conduction current 


j 


o 


n co y qE * 


Thus experimental knowledge of the dark conductivity 


a = n u e 
o co o 
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constrains the product Field enhanced conductivity 

is introduced by allowing the trap ionization rate coeffi- 
cient v to depend on the electric field. The Adamec and 
Caldwood model of Section 3.7 may be obtained by replacing 
v q in Equation (3.39) by 

y- [2 + cosh ( 3 f E 1//2 /2 kT)] . 
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3.9 ONE-DIMENSIONAL MODEL FOR DIELECTRIC EFFECTS IN A 
CHARGING ENVIRONMENT 

Calculations have been performed on the charging of a 
one-dimensional system consisting of a conductor coated with 
a thin layer of dielectric. The problem geometry is schema- 
tized in Figure 3.10. The conductor is assumed negligibly 
thick. 



conductor 


Figure 3.10 


Here, j is the undisturbed plasma current, part of which is 
reflected if the surface potentials V or V + are negative; 
in that case, the current incident on the dielectric surface 
is 



where 0^ is the temperature of the assumed Maxwellian plasma. 

A similar result applies for the current incident on the con- 
ductor. The photo-currents emitted by the dielectric and con- 
ductor are denoted by j and j*, respectively. Secondary 
emission caused by impact of a single electron is represented 


63 



by 6 ± . Here no attempt is made to account for the energy 
dependence of 6 ± , even though it would not be difficult to 
do so. 

The total particle currents incident on the slab are 



where Q & is the "temperature" of the emitted electrons, and 

f (y) - 1 , y > 0 

5=5 exp (y) , y < 0 . 

In addition to the kinetic equations describing charge trans- 
port in the dielectric, we have (by differentiating Equation 
(3.30a) with respect to t and integrating with respect to x) 


e (x,t) + qj (x,t) = J(t) 


(3.50) 


with 

-19 

q = -1.6 x 10 coulomb 
j(x,t) - j“ x < 0 

- jp x > l (3.51) 

* j d (x,t) = j b (x,t) + j c (x,t) 0 <_ x <_ i 

where j. is the current of plasma electrons penetrating into 

the dielectric, j is the dielectric particle current, and 

c 


64 



J(t) is an integration constant. In the plasma , we neglect 
space charge so that if the plasma "ground planes" are at 
x = -D a nd x * D + i, then the space fields are independent 
of position: 


E(x,t) = E~ = - — - x < 0 

+ v + 

E (x, t) = E = — x > l . 


(3.52) 


Integrating Equation (3.50) over space from x = -D to 
x = D + l, using the boundary condition V =» 0 on the plasma 
ground planes, gives 


J = q( 2 D + i) 1 |o(j+ + j") + i J j d (x,t)dxj (3. S3) 

where k = e/e Q is the relative dielectric constant of the 
* material . 

The body charging time t discussed earlier is in 

W 

general very short in relation to the time scale for develop- 
m ment of potential differences across the thickness of di- 

electric. In the present circumstances t is determined by 

c 

the capacitance per unit area of the dielectric metal slab 
relative to the plasma and is proportional to 1/D. If the 
« time t elapsed after commencement of charging is much greater 

than t c , then we expect that the plasma electric field will 
vary slowly, so that it is a good approximation to neglect 
the vacuum displacement currents e Q 3E/3t in Equation (3.50), 

0 and obtain from Equations (3.49) - (3.50) algebraic relations 

between the surface voltages V ± and the circuit current J. 

Thus, for example, if V* < 0, 


* 
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+ K - j /q 

V = e £n — — 

p 3 p (l-6 + ) 


J - qj 


__ 3J 

+ at 


J + gj 


at 


Now, integrating Equation (3.50) over the dielectric gives 


It < v ‘ 


r A -j 

- V + ) = i Jl - f j d (x,t)dx 


or, using Equations (3.55) - (3.56), 


tJ - qj d ] [J + qj J U - qj v l 
q(j* + j“) 


Xt 

J J j d (x,t)dx 


Equation (3.57) shows that a steady state for the sys- 
tem occurs when J/q = Equation (3.58) shows that, if 

|j| << | q j “ | 'v qj^ this steady state is approached with a 
time constant 
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(3.60) 


t = e9 p 
d *qj v * 

This is the time constant associated with development of a 
potential difference across a dielectric. With e * 3 x 10~^ 
farad/meter, 0^ ^ 10^ volts, qj v 'v 10 ^ amp/meter^, and 
i = 10 * meters 

t^ z 300 sec. 

Such differential charging time scales enter together with 
the time constants associated with the kinetics of the charge 
carriers in the determination of the electric fields within 
the dielectric medium. 

A further useful relation is obtained by using the 

field equations to eliminate j* and j~ from Equation (3.53); 

P P 

<e a . 

J 53 “ — j— (V - V ) + q3^ . (3.61) 

The field in the dielectric then satisfies 

e • 3 d> - nr It (v+ " v "> • »•«> 
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3.10 NUMERICAL METHODS FOR ONE-DIMENSIONAL DIELECTRIC 

CALCULATIONS 

3.10.1 Methods 

The results of numerical calculations reported in this 
section are obtained from difference approximations to equa- 
tions developed in the preceding two sections. The difference 
equations will not be displayed, but we will discuss briefly 
means for efficiently obtaining solutions through several 
minutes of charging time. The basic limitations on a prac- 
tical numerical scheme are imposed by the large disparity 
between the various time scales that occur in the defining 
equations; the potentially most stringent limitation is 
associated with the drift current in Equation (3.38). The 
code used to perform these calculations has not been in- 
corporated into NASCAP, but has been delivered separately to 
NASA/LeRC . 

Equations (3.36), (3.38) and (3.41) are taken as de- 
fining equations in the dielectric. In Equation (3.38) how- 
ever we neglect 3j_/3x, which in the examples to be con- 

C 

sidered is small in comparison with the remaining terms on 
the right hand side of the equation. One extremely useful 
consequence of this approximation is that it permits a sub- 
stantially larger time step in the numerical scheme than 
would otherwise be possible. 

The algorithm for the particle conduction current at 
the grid point k(x = kAx) is 

j c (k) = j vE(k) (n c (k+1) + n c (k)) 

-|u|E(k)| (n c (k+l) - n c (k)] . (3.63) 
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If n and E are staggered on a uniform spatial grid, this 

C 

algorithm is first order accurate in Ax; without the term 

proportional to |E(k)| it would be second order accurate. 

The latter term, together with a Courant restriction on the 

time step, which is required for stability, assures that the 

difference equations maintain the particle densities as 

positive quantities. If the term 3j /3x were retained in 

c 

Equation (3.38), the condition 


At < At 


_ 1 Ax 
1 2 pE 


would be sufficient for stability. Neglecting 3j /3x, the 

c 

sufficiency condition is relaxed to 


At S At2 


1 

7 


Ax 


eff 


E * 


Using y • 10 cm 2 /volt sec, y -- ■ 10 0 cm 2 /volt sec, 

-3 5 err 

Ax * 2 x 10 cm and E * 10 volts/cm, gives 


At^ * 10 ^ sec 
2 

Atj =* 10 sec 


A time step limitation as small as At^ would be impractical 
for calculations which extend over several minutes of 
charging time. 

Other time step limitations which could occur in an 
explicit time-differencing scheme are for all practical pur- 
poses removed with an implicit scheme. For example, the 
time scale 

t » N fc <av> 



does not limit the stability of the difference scheme when 
the first term in the difference analogue of Equation (3.38) 
is evaluated at the advanced time. 


3.10.2 Results 

Three sets of calculations were performed with the 
following parameters 

D * 10 cm 

l = 10" 2 cm 

-10 2 

jp = 10 amp/cnr 

= 0.75 x 10 ^ amp/cm 2 

- -10 2 
j = 0.50 x 10 amp/cm 

<5~ = 0.1 

6 + = 0 

N t » 10 18 cm” 3 

4 

0 =10 volts 

P 

0 =2 volts 

e 

k = 2 

The dielectric was divided into 50 spatial zones, each 

-4 

having Ax = 2 x 10 cm. The beam current profile within 
the dielectric was assumed to be linear, dropping to zero 

_3 

in a distance of 10 cm. 

The three cases considered were (1) a passive di- 
electric, (2) a dielectric with natural and field enhanced 
conductivity, and (3) a dielectric identical to that in 
Case 2 but having in addition a radiation induced conducti- 
vity in the deposition zone. 

In the passive dielectric, charge deposited in 
trapping sites remains in the deposition layer, giving a 

70 
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charge density profile which is spatially uniform in the 
deposition zone, and which vanishes beyond. 

The additional parameters required to define Case 2 

are 

-4 2 

y “ 10 cm /volt sec 

v ** 3 x 10 3 sec" 1 

-9 3 -1 

<ov> ** 10 cm sec 

8 - .021 cm 1/2 / (volt) 1/2 

n,^ = 6.25 x 10 4 cm -3 

n to = 2 * 083 x 1q11 cm" 3 

corresponding to a dark conductivity 

—18 

cr Q ** n cQ yE “10 mho/cm 

and an effective mobility 

v -io 2 

p eff * N ~ <ov> y = 3 x 10 e* / volt sec * 

T 

Case 3 is identical to Case 2 except for addition of a 
"radiation induced conductivity" in the deposition zone. 

The ionization rate (cm -3 sec" 1 ) is taken to be 

10’ 3 j b (x,t) n t (x,t) 

—2 —1 —3 

where j fa is electron flux (cm sec ) and n fc (cm ) is the 
density of occupied traps. The induced conductivity when 

<3j»- a 10 ^ amps/cm 2 , n. * 2 x 10 11 cm 3 is approximately 

u -15 

o r s 2 x 10 mho /cm. 

The surface potentials V~ and V + as a function of 
time, plotted in Figure 3.11, are very nearly equal for all 




versus time 



three cases for t <_ 230 sec, the maximum time for which 
Cases 2 and 3 were run. The small differences between the 
dielectric potential drop in the three cases is illustrated 
by Figure 3.12. There does however appear a trend toward a 
departure from Case 1 for times larger than a few hundred 
seconds. That this might be expected is indicated by com- 
paring the electric field and charge density profiles given 
in Figures 3.13 and 3.14, respectively. The latter figure 
for Cases 2 and 3 shows a substantial charge migration away 
from the region of deposition. 



Potential drop (volts) 



Time (sec) 


Figure 3.12. Dielectric potential difference (V + - v“) 
versus time. 
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Case 



Distance into dielectric (x/t) 

Figure 3.14. Charge density versus distance at 210 sec 




4. TEST CASES FOR THE NASCAP CODE 

For the purposes of illustrating, error-checking and 
verifying the NASCAP code, test cases were run in all three 
NASCAP operating modes (see Table 4.1). The first sequence 
of problems were designed to simulate the electrostatic 
charging of material samples in a laboratory test facility. 
These ground test model test cases were done with one inch 
resolution in a computational space with dimensions in rough 
correspondence with those of the LeRC facility. An electron 
beam profile similar to one measured at LeRC was used. The 
last sequence of problems were designed to simulate the 
charging of an object in a plasma with parameters similar 
to those found in the earth's magnetosphere. These space 
model test cases, with two exceptions, were spheres in iso- 
tropic environments. The material properties were those sug- 
gested in the NASCAP User's Manual. 

4.1 GROUND TEST - FLOATING ALUMINUM PLATE 

This set of test cases were entirely successful. No 
problems were encountered concerning length of tine step, and, 
because the potential scales with total charge, repeated 
potential calculations were not required. The computer costs 
involved were quite modest. 

The electron gun emitted a current of 0.37 pA and had 
a profile similar to that supplied to S^ by NASA-LeRC (see 
Figure 4. Id). The peak flux was slightly greater than 
1.0 nA/cn . The sample was a 6 inch by 8 inch plate, 1 inch 
thick, located 40 inches from the electron source. A mag- 
netic field comparable to the earth's field (B = 0, B = 

x y 

0.52 gauss, B = -0.19 gauss, where y is vertical and z is 
z 

the beam propagation direction) was assumed present. Simu- 
lations were carried out for beam energies of 2, 5, 0 and 
20 keV. 
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Table 4.1. NASCAP Test Cases 
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The simulation proceeds at each cycle by first tracking 
particles forward in the tank (Figure 4.1a-c) and thus deter- 
mine the current density and incident angle at the sample 
(Figure 4.1d-f). The curvature due to magnetic field is 
apparent at the low energies. The charge on the plate is 
adjusted in accordance with the net current, the potentials 
scaled accordingly, and the next cycle is begun. 

For each beam energy the plate charged to its final 
value in a time roughly proportional to the beam voltage 
(Figure 4.2). The difference between the beam voltage and 
the final plate potential increased somewhat with increasing 
energy. Final current balance was achieved by a decrease in 
incident current to about half its uncharged value, and a 
substantial increase in secondary emission ratio, attributable 
in part to non-normal incidence. Decrease in incident electron 
flux was more important at higher energies (Figure 4.2, 

Table 4.2) . 

4.2 GROUND TEST - TEFLON COATED PLATE 

Another sequence of ground tests was performed with a 
5 mil teflon coating on a grounded aluminum plate. The beam 
characteristics, experimental geometry, and magnetic field 
were identical to the previous case. Since simple potential 
scaling is not appropriate, 10-30 potential iterations were 
performed each cycle following use of the "DSCALE" option. 

This case differs from the previous in that (1) because 

the charging is across a larger capacitance the time scale is 
2 

'''lO times as long, and (2) the potential varies with position 
on the insulating teflon surface. The results are given in 
Table 4.3 and Figures 4.3 and 4.4. The central area of the 

4 

sample charged to -1.8 x 10 volts for the 20 keV beam, and 
-7.8 x 10^ volts for the 10 keV beam. In both cases the 
periphery of the sample initially charged at about half the 
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Potential (Volts) 


20 . keV 


f 

2 keV 


20 keV 


,2 I / 2 keV 


Ground Test 
Aluminum Plate 


20 ,keV 


2 keV 


Tine (Seconds) 


Peak incident flux (lower curves), net current 
(upper decreasing curves) , and potential (in- 
creasing curves) versus time for an aluminum 
plate subject to an electron gun with energy 
2, 5, 8, 20 keV. 


Peak Flux (nA/cm ) Net Current (na) 
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Figure 4.3. Time development of potentials, net current, 
and peak incident flux for a 5 mil teflon 
coated plate subject to a 10 keV electron beam. 
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Figure 4.4. Time development of potentials, net current, and 
peak incident flux for a 5 mil teflon coated 
plate subject to a 20 keV electron beam. 
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rate of the center, and tended to catch up later in the simu- 
lation. The incident current was reduced to about half its 
initial value at the end of both simulations. Space poten- 
tial contours at the end of the 10 keV simulation are shown 
in Figure 4.5. 

4.3 GROUND TEST - SSPM 

The final set of ground test cases exposed to the 
electron beam a complex object (Figure 4.6) resembling the 
Spacecraft Surface Potential Monitor (Experiment SCI) 
scheduled to be flown on the SCATHA satellite. The object 
consisted of a 14 inch x 14 inch x 1 inch aluminum plate with 
four 5 inch x 5 inch material samples: teflon, kapton, Si0 2 , 

and clean magnesium. The properties of these surfaces were 
those suggested in the NASCAP User's Manual. The insulators 
were 5 mils thick, and the magnesium sample was mounted on 
an insulating spacer with a capacitance of 347 pf. These 
test cases were run in similar fashion to the teflon plate 
cases, except that a wider beam profile was used, and the 

2 

beam current was increased to 0.84 yA to maintain a 1.0 nA/cm 
flux at the beam center (Figure 4.7). The potentials and 
fluxes at the central surface cell of each sample were moni- 
tored. 

The results are shown in Table 4.4 and Figures 4.8 to 
4.15. For the 2 keV case it is seen that the teflon and Si0 2 
initially charge positive, while the kapton and magnesium 
charge negatively. This suggests the presence of a photo- 
electron current between the surfaces. NASCAP could have 
completed this case through the use of a very short time step, 
but the cost would have been excessive. 

For the 5, 8, 10 and 20 keV cases all surfaces charged 
negatively, with the magnesium surface charging most rapidly 
due to its smaller capacitance to the aluminum. Often the 


87 



POTENTIAL CONTOURS ALONG THE Y-Z PLANE OF X 



i 


Figure 4.5. Potential contours at conclusion of 10 keV teflon plate tank test case. 
The potential contours are plotted in intervals of 500 volts. 
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Table 4.4. SSPM Ground Test Case Results (Current Balance, nA/cm 
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field due to the magnesium was sufficient to reverse the sign 
of the field in front of the other samples , suppressing their 
secondary emission. The simulations were carried out well 
beyond saturation of the magnesium potential. The oscillatory 
behavior seen in many of the plots is due to having taken ex- 
cessively long time steps in order to minimize computer costs. 
Because of this problem, the runs were stopped before the 
insulating samples became fully charged. 


4.4 ALUMINUM SPHERE SPACE TEST CASES 

Several tests were performed on an aluminum sphere of 
nominal diameter 3 meters. The actual object was the smallest 
quasisphere definable by NASCAP, having one facet in each of 
the 26 symmetry directions. The effective diameter (from 
capacitance calculation) was 3.19 meters. Because of its 
symmetry and because it is a simple conductor, the aluminum 
sphere is a particularly simple object for a charging simula- 
tion. 


4.4.1 Maxwell Probe Calculation 

The Maxwell probe formulation used in NASCAP is, in 
principle, exact for a sphere in an isotropic, Maxwellian 
plasma, and thus nearly exact for a quasisphere in such an 
environment. This formulation writes the differential 
particle flux per unit area to the satellite as 

ftess ’ (, (is) 1/2 E “p'-' e + ^vi/M] £2si 

where E is the energy of incidence, q the particle charge, 
and V the surface voltage. (The formula requires 
E + qV > 0.) The incident flux is then given by 


100 



exp ( -qV/kT ) qV > 0 



(1 - qV/kT) qV < 0. 


The secondary emission and backscatter are calculated as 


00 x 2 

/ dEx2 *J a(cose) if(E,e) dEdQ 


min 


where y is the relevant coefficient. (The angular integral 
has been performed analytically, whereas the energy integral 
is done using Simpson's Rule.) 

In the environment shown in Table 4.5 the sphere 
reached an equilibrium potential of 1415 volts in a time of 
0.2 seconds (Figure 4.16). The initial and final current 
balances are shown in Table 4.6. Note the important role 
played by the cold ions and the resulting secondary 
electrons in establishing the final current balance. 

4.4.2 Reverse Trajectory Simulation — Isotropic Flux 

Three reverse trajectory simulations were carried out 
using successively finer grids of incident energy and angle. 
The environment was identical to the previous case (Table 
4.5). The potentials reached after 0.81 seconds were 
-1720 volts, -1640 volts, and -1600 volts for (n £ , n Q ) - 
(3,3), (4,4), and (5,5) respectively (Figure 4.17). The 
degree of charging is overestimated due to underestimation o 
of the secondary emission. A more nearly optimal choice 
of incident angles would serve to ameliorate this situation. 


0 
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Table 4.5. Maxwellian Environment for Space Tests 


Species 

e 

H + 

Temperature 

4.11 keV 

430 eV 

Density 

0*53 cm ^ 

0.60 cm 
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Net Flux (nA/cm 



Time (Seconds) 

Figure 4.16. Potential and net flux for 3 m aluminum sphere 
in environment of Table 4.5. 
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Table 4.6. Current Balance for Aluminum Sphere — 
Maxwell Probe Calculation 



Initial 

Final 

Potential (volts) 

0 

-1415 

2 

Incident Electrons (nA/cm ) 

-0.0909 

-0.0645 

Backscatter (nA/cm 2 ) 

0.0303 

0.0215 

Secondaries (nA/cm 2 ) 

0.0388 

0.0277 

. 2 
Incident Protons (nA/cm ) 

0.0008 

0.0033 

Secondaries (nA/cm 2 ) 

0.0018 

0.0120 

Net Flux (nA/cm 2 ) 

-0.0192 

0 
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Figure 4.17. Reverse trajectory simulations for 3 m aluminum 
sphere in environment of Table 4.5. Figures in 
parentheses indicate numbers of incident energies 
and angles used. The dashed lines are the re- 
sults of the Maxwell probe simulation. 
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4.4.3 Reverse Trajectory Simulation — Anisotropic Flux 

The final aluminum sphere space test case uses a 
reverse trajectory treatment for an anisotropic plasma based 
on ATS-5 data for hour 9.998 of day 73. The plasma tempera- 
tures were similar to those of Table 4.5, but the electron 
density was somewhat lower and the proton density much higher. 
The incident particle matrix was 5 * 5 * 5 (energies, polar 
angles, azimuthal angles). 

The results are shown in Figure 4.18. The final 
potential was -280 volts, with a nonuniform flux (averaging 
to zero) over the surface of the sphere. 

4 . 5 INSULATED SPHERES - MAXWELL PROBE CALCULATION 

Kapton and teflon both tend to achieve positive poten- 
tials in the environment of Table 4.5. Since the electric 
field just outside the surface suppresses low energy electron 

emission, the maximum potentials are only a few volts and are 

-4 

achieved xn ^10 seconds. The results of these simulations 
are given in Table 4.7 and Figure 4.19. 

4.6 COMPARISON OF DeFOREST DATA AND MAXWELLIAN 

Reverse trajectory calculations (Figure 4.20) were 

done for a 3 m teflon sphere subject to the ATS-5 data for 

hour 9.998 of day 73, and the similar environment of 

Table 4.5. Despite an overshoot in the Maxwell simulation, 

-4 

both charge to *vl.2 volts in 'vlO seconds. 
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Figure 4.18. 


Simulation of an aluminum sphere in an anisotropic 
plasma, with net flux to two representative surface 
cells. 
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Table 4 . 7 . Current Balance for Insulating Spheres in Environment of Table 
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Time (Seconds) 

Figure 4.19. Net fluxes and potentials for dielectric-coated spheres in environment 
of Table 4.5. 



Net Flux (nA/cm 



Time ( Seconds ) 

Figure 4.20. Reverse trajectory calculations for 3 m teflon 
sphere using DeForest data (hour 9.998 of day 
73) and data of Table 4.5. 
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4.7 


ISOTROPIC SUN 


A special version of NASCAP was mapped for which the 
sun shone on a 3 m teflon sphere uniformly from all direc- 
tions. (Normally/ NASCAP does a full shadowing calculation 

for photoemission.) The flux was dominated by the photo- 

2 

current of 2 nA/cm . The sphere charged to 16 volts m a 
few times 10 - '* seconds (Figure 4.21). 

The reason the sphere was charged to well above the 
2-volt characteristic energy of secondary- or photo-electrons 
is worth some explanation. This run was performed using the 
default "NOSHEATH" option to avoid tracking secondary 
electrons. Under these conditions, the low energy electron 
current emitted by an electron attracting surface is re- 
duced by a factor 

|£| < 10 
|E| > 10 

where |e| is the electric field normal to the surface in 
volts/mesh unit. Thus the satellite charged to the point 
where its surface field was approximately 10 volts/mesh unit, 
the mesh unit in this case being 1 m. Had the "SHEATH" 
option been invoked, 2 eV electrons would have been emitted 
from the sphere and tracked in the electrostatic field. 

When the surface potential exceeded 2 volts positive, 
these electrons would have returned to the sphere. NASCAP 
would then have cancelled the emitted and return currents, 
thus giving an equilibrium potential slightly in excess of 
2 volts. 


f(E) = 


exp (- | E | /2) 
0 



Potential (Volts) 



Time (Seconds) 

Figure 4.21. Potential and net flux to a 3 m teflon sphere 
exposed to sunlight uniformly over its surface 



4.8 


SPACE TEST OF SSPM 


The final test case was the SSPM (Figure 4.6) in 
the environment of Table 4.5 and mounted on an aluminum 
plate with potential fixed at -575 volts. (This was the 
potential achieved in a test run of a SCATHA-like satel- 
lite. See Appendix E.) The potentials and fluxes are 
shown in Figures 4.22 and 4.23. After 100 seconds, the 
potentials were: magnesium at -1500 volts, kapton at 

-500 volts, Si02 at -410 volts, teflon at -380 volts. 

Beyond this time, the magnesium (which eventually reached 
-2200 volts) was sufficiently negative to suppress low 
energy electron emission from the insulators, and the 
simulation became unstable to the long time step which 
was taken, again for cost considerations. 

4.9 CONCLUSIONS 

The variety of cases presented here, together with 
that presented in Appendix F, comprise the first tests 
of the NASCAP program. They were designed to exercise 
much of the modeling capability within the code and to 
point out the strengths and weaknesses of the techniques 
employed in NASCAP. These calculations were performed at 
the end of this contract period and there was little time 
remaining to modify the analytical models to improve speed 
or accuracy. No attempt was made to obtain improved output 
by re-performing a test case after analysis of the first 
simulation. 

In general the code worked remarkably well. That 
is, the collection of physical models which comprise NASCAP 
were able, without modification, to calculate the charging 
processes and yield physically reasonable results. Both 
ground test and space models worked smoothly in almost all 
instances. Flux, potential, material, and electrical 
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Time (Seconds) 

Figure 4.22. Potentials for SSPM material samples in environ 
ment of Table 4.5 and mounted on an aluminum 
plate fixed at -575 volts. 
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models functioned extremely well. Only one test case offered 
any difficulty; that was the 2 keV SSPM tank test. In that 
case the extremely large secondary currents from the di- 
electrics along with the magnesium negative potential 
created an electron sheath situation that the present 
NASCAP sheath treatment could calculate only by using ex- 
tremely short time steps. However, on account of budgetary 
constraints it was decided not to rerun the case, but to 
present it as a limitation on the current version of NASCAP, 
something that simple sheath model modifications could 
eliminate. 

These tests demonstrate that NASCAP can calculate 
successfully almost all aspects of the three-dimensional 
electrostatic charging of materials both in ground test and 
space environments* It should have great value as a design 
and analysis tool for scientific and engineering applica- 
tions. 
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A PRELIMINARY SPECIFICATION OF THE ENVIRONMENT 
AT GEOSYNCHRONOUS ALTITUDE 


Prepared by 

MAYA DEVELOPMENT CORPORATION 


for 


SYSTEMS, SCIENCE AND SOFTWARE 


August 1976 
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1.0 MAGNETOSPHERIC ENVIRONMENT 

The magnetic field of the earth is confined to a finite volume by 
the pressure of the solar wind plasma which distorts the field into a 
teardrop-like shape with a tail of indefinite length (see Figure 1.1). Within 
this volume, called the magnetosphere, satellites encounter a wide variety of 
plasma physical phenomena. During the almost two decades since the 
initiation of space exploration a truly enormous amount of information on the 
magnetospheric plasma, its dynamics and its effects upon the earth, have been 
collected. From this vast array of data we have selected a data base from 
the geosynchronous satellites ATS5,6 which will be used to describe the 
magnetospheric environment. The last ten years has seen the development of 
theoretical models of the magnetospheric plasma which also guide our 
interpretation of the data. With this data and these models we can 
adequately describe the magnetospheric "weather" which is likely to cause 
spacecraft charging to occur. 

The magnetosphere is a vast relaxation oscillator gradually distorting 
22 

and storing some 10 ergs drawn from solar wind and earth's rotation and then 
explosively releasing the stored energy into the polar atmosphere where 
visible auroras occur. This phenomenon is known as a magnetospheric substorm. 
In general the substorm reaches its peak within 15-20 minutes and then 
gradually subsides. Auroral activity follows this pattern of growth and decay 
and also indicates the concentration of most subs torn activity between local 
midnight and dawn. 

The visual aurora is but one of many manifestations of the magneto- 
spheric substorm. Observations from the ground, from balloons and from 
rockets have shown that negative bays in the horizontal magnetic component, 
enhanced absorption of cosmic radio noise, the production of millisecond 
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Fig. 1.1 The Magnetosphere of the Earth 


X-ray bursts, the occurrence of VLF hiss and chorus, the pattern of micro- 
pulsation observation, are all correlated with the onset of a magnetospheric 
substorm. Direct measurements from rockets and satellites have shown the 
visual aurora to result from fluxes of electrons and protons that occur during 
the early phases of a substorm. The ground based data and general description 
of the magnetospheric substorm are best reviewed in Akasofu (1968). 

Satellite observations provided new insight into the phenomenology 
and dynamics of the magnetospheric substorm. Using data from Ogo 1,3 and 
Vela, Vasylunas (1968) demonstrated the existence of an intense low energy 
flux of electrons (the plasma sheet). The inner boundary of this plasma 
sheet was found to move inward with the onset of a substorm as indicated by 
ground based magnetic data. Lezniakand Winkler (1971) used ATS-1 electron 
data in energy ranges 50-150 keV, 150-500 keV, and 500-1000 keV to demonstrate 
the convective injection of energetic electrons into synchronous orbit. 

These concepts were extended and the data base to support them was vastly 
enlarged by the work of DeForest and Mcl Twain (1971) who used electron and 
proton data from ATS-5 differential analyzers that measured some 64 energy 
levels between 50 eV and 50 keV. This satellite provides much of the data 
base which we propose to use. 

The theoretical understanding of plasma flow at geosynchronous orbit 
was further extended by Mcllwain (1972) who used static electric and 
magnetic fields to model the magnetosphere (see Figures 1.2 and 1.3). 

Tracing particle orbits in these fields enabled him to explain many of the 
features observed in the ATS-5 data. In particular he was able to explain 
the shapes of boundaries which are often observed between particles of 
different energies. This work also demonstrated that intense fluxes 
correspond to particles which have been convected in from regions of low 
magnetic field (<40 y) . 








Extensions of this work were presented by Mcllwain (1974) and by 
Mauk and Mcllwain (1974). In this work the concept of an "injection boundary" 
similar in many ways to the plasma sheet boundary of Vasylunas (1968) was 
added to the models of the earlier work. By tracing the trajectories of 
particles backward through the model fields one is able to determine the 
inner boundary of the fresh plasma. The injection boundary was found to be 
related to kp, the universal index of magnetic disturbance, and is given by 

122 - 10 K 

R b = $ - 7.3 P * W 

when 4>, the local time, lies between 18 and 24 hours. The concepts have 
been examined most recently by Konradi , et al (1975), who have found that the 
injection boundary and Mcllwain field model explain their observation of 1-35 
keV protons and 1-300 keV electrons from Explorer 45. 

The launching of ATS-6 again added to the measuring capability of 
satellite plasma analyzers. The instruments on ATS-6 enable one to measure 
electrons and protons over the range of 1 eV to 80 keV and can in addition 
be scanned mechanically to look at different azimuthal and pitch angles. 

With these detectors new phenomena inaccessible to study with ATS-5 are being 
examined. In particular, Mcllwain (1975) has found intense field aligned 
fluxes of electrons. These fluxes usually occur shortly after the onset of a 
magnetospheric substorm. They are characterized by an energy spectrum 
which is flat out to a break point energy usually between .1 keV and 10 keV. 
Those beams with break point energies above 2 keV seem only to occur within 
the first 10 minutes after the onset of plasma injection associated with a 
substorm. These beams will require further study and could be particularly 
effective in producing differential charging. 

The charging of spacecraft to kilovolt potentials was first 
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discovered in ATS-5 data and reported by DeForest (1972). The problem of 
spacecraft charging had been recognized early in the design of the ATS-5 
instruments. Before launch, the research group at San Diego under the 
direction of S. E. DeForest and C. E. Mcllwain had insisted that conducting 
collars be placed around the apertures of the low energy particle detectors. 
In addition some concern was expressed about the fact that the viewing cones 
of two instruments looked out through a cylinder of solar cells which could 
charge to high potential. 

At the time of the initial discovery of kilovolt charging potentials, 
only a few volts were expected. Thus the large potentials were somewhat 
surprising. DeForest (1972) was able to piece together an elementary theory 
which explained the observed high potentials. 

The pre launch worries about charging of the solar cells proved well 
founded. Differential charging, indicated by a spin modulation of the 
parallel detector fluxes, could only be caused by local electric fields. Thus 
DeForest (1973) demonstrated experimentally the possibility of kilovolt 
differential charging. Further work on this subject was presented at a joint 
AGU/AIAA special session on spacecraft charging in 1975. An especially useful 
paper, which we shall use to help select environmental data, was presented by 
Reasoner, et al (1975). It discussed the statistical relationship between 
the ATS-6 spacecraft charging events and the encounter of warm plasma clouds 
associated with the onset of a magnetospheric substorm. 

2.0 DATA SELECTION 

Environmental data obtained by plasma spectrometers on board the 

Advanced Technology Satellite 5 (ATS-5) from September 1969 through the 

vernal equinox of 1971 has been analyzed to provide input spectra for use in 

SSS spacecraft charging programs. Representative data from six days is 

presented in the following forms: 
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• 24 hr Spectrogram 

• Integrals, ne, ni , je, ji, the number densities and energy 
fluxes respectively based on 2.3 minute averages for the 
selected 24 hour periods 

• Plots of six selected spectra for each 24 hour period 

• Printouts and punched cards containing the selected spectra 
The data has been selected to typify several broadly different 

categories of magnetospheric weather which occur at geosynchronous orbit. 
While representative of the magnetospheric conditions the data is not 
extensive. It is meant to provide useful input for the development of space- 
craft charging codes. A complete meteorological survey is being sponsored 
by Air Force Cambridge Research Laboratories. Data for this report has been 
provided under Defense Nuclear Agency Contract No. DNA001-76-C-0121 . 

2.1 BACKGROUND INFORMATION FOR DATA SELECTION 

Previous studies provide background information which is useful in 
the selection of data for the study of spacecraft charging. Facts which 
should be considered include: 

2.1.1 SPACECRAFT PERFORMANCE FEATURES 

• Spin up anomalies on the DSCS-II spacecraft are well 
correlated with geomagnetic substorms (TRW SCA II, 1975) 

• There is a strong association of unexplained satellite 
performance with the midnight to dawn sector of local 
time (McPherson, et al (1975)). (See Figure 2.1) 

• The local time distribution of spacecraft charging events 
is found to maximize between local midnight and dawn 
(Reasoner, et al (1975)). (See Figure 2.2) 
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Figure 2.1 Unexplained Satellite Performance Correlates with Local Time. 
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Figure 2.2 Local Time Distribution of 



2.1.2 MAGNETOSPHERIC WEATHER FEATURES 

• Equatorial observations by the geostationary satellite ATS-5 
of charged particles on auroral lines of force reveal the 
frequent injection of plasma clouds into the magnetosphere. 
These intrusions of hot plasma are found to have a one to one 
correspondence with magnetospheric substorms. (DeForest and 
Mcllwain, 1971) 

• The elctromagnetic fields surrounding the earth act to 
separate the injected plasma clouds on the basis of both 
charge and energy. (Mcllwain, 1972). Electric fields attempt 
to bring about corotation of low energy electrons and protons 
as one moves inward toward the earth. Magnetic field 
gradients cause high energy electrons to drift toward the 
dawn side of the magnetosphere while high energy protons are 
caused to drift toward dusk. Effects of importance which 
follow are first, the spectra in the midnight to dawn sector 
are characterized by high electron energies and thus tend to 
induce spacecraft charging (DeForest, 1972) 

Second, field aligned fluxes are set up in order to maintain 
overall charge neutrality. (Mcllwain, 1975). These fluxes 
can make important contributions to differential charging of 
spacecraft surfaces. (DeForest, 1973) 

• The plasmasphere shrinks during periods of high magnetic 
activity (Chappel , 1970). At geosynchronous orbit encounters 
with the plasmasphere are concentrated in the local noon to 
local evening sector as shown in Figure 2.3. Plasmasphere 
encounters are anti -correlated with spacecraft charging. 
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Figure 2.3 Local Time Distribution of Warm 


Plasma Encounters. 
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There are two reasons for this. First, the high density 
low energy plasma provides a grounding current to the space- 
craft thus preventing large potential buildups. Second, 
plasmasphere encounters are more common during quiet times 
when substorm activity is low. 

2.2 TYPICAL CONDITIONS REPRESENTED 

With these facts in mind data from the year 1970 gathered by the 
UCSD plasma spectrometers on ATS-5 was analyzed. From this data the six 
representative days were chosen to typify the following magnetospheric weather 
conditions. 

• A quiet day with no substorm activity 

• A moderately active day with a single substorm of low intensity 

• Two days with intense localized post midnight substorms 

• A premidnight substorm 

• A day when spacecraft charging occurred in the sunlight 

2.3 SPECIAL CONDITIONS 

Several special conditions can occur which we have designed into the 
total distribution functions which are to be used in this study (See 
Appendix C). 

2.3.1 ECLIPSE OF THE SUN - the spacecraft charging phenomena was 
first discovered on ATS-5 during eclipse. The loss of the 
large photo-electron flux from the satellite surfaces allows 
the satellite to float up to high potentials (DeForest, 1972) 
FIELD ALIGNED FLUXES - usually encountered during the 
intense early phases of a substorm and of importance because 
they can cause differential charging. 
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3.0 DATA AND FORMAT 

The data presented in this section consists of four types in addition 
to sets of punched cards for easy computer use. For each day or event 
presented, a 24-hour spectrogram is used to establish the context. Following 
that, selected spectra are shown which illustrate significant events. 
Instruction for reading both the spectrograms and the average plots are 
provided in the appendix. Punched cards containing the same data are provided 
with this report with instructions for their use. The printouts of the 
detailed spectra are also provided. Note that these detailed spectra are 
produced from 6.8 minute averages. This gives good statistics without 
smoothing rapid time fluctuations. 

The final form of data presentation is a table of various integrals 
taken over 2.3 minute periods for the whole day. These might prove useful 
for studies where analytical approximations to the spectral shape is more 
useful than the actual spectra. The whole day is provided for possible 
future use in time-varying codes- 

3.1 MODERATE ACTIVITY 

2/ 1/70 The activity on this day was limited to two early morning 
injections at approximately 0110 and 0500 UT. The effect was to bathe the 
spacecraft in a moderate flux of 3000 volt electrons. From previous 
experience, we can estimate that had the spacecraft gone into eclipse on 
this day, it would have charged to approximately 1000 volts. 

The fluxes associated with this injection were insufficient to 
cause charging in the sunlight. [For purposes of this report a potential of 
less than about 50 volts will be neglected since the ATS-5 detectors do not 
sense lower energies]. Furthermore, isolated substorms of this type have 
never been seen to charge ATS-5 significantly. However, from ATS-6 data we 
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can estimate that a potential of at least -5 volts occurred and that by 
simply renormalizing the total flux by a factor of 2-10 while keeping the 
same spectral shape, we can simulate the conditions under which daylight 
charging of -100 volts would happen. 

Detailed spectra are provided for 0300 to see the pre-electron 
encounter conditions. The next spectra are at 0400 when the high-energy 
protons had been encountered, but not the associated electrons. The next 
three sets of spectra are taken at different points in the main part of the 
substorm. At 0530 ATS-5 experiences the greatest flux of high energy 
electrons. By 0630 the average energy of the electrons has fallen slightly 
. due to gradient drifting while the average energy of the ions has 

increased slightly. At 0730 the ion chasm is well developed, and a notch has 
developed in the electron spectra. This feature is common and will persist 
for the entire day. A final set of spectra taken at 1200 is provided simply 
to complete the story. The spectra at 0530 and 0630 are probably the most 

hazardous to the spacecraft. 

In summary, 2/1/70 is a good example of isolated, moderate activity 
which could be used to study the response of a spacecraft to a normal environ- 
ment. 

3.2 INTENSE LOCALIZED POST MIDNIGHT SUBSTORM 

2/11/70 On this day we were fortunate enough to find an intense 
substorm occurring right at the spacecraft location. This day is 
particularly valuable because of the lack of complicating activity at other 
times, and because no corrections for daylight charging are needed. 

The injection took place at 0850 when ATS-5 was located in the 
hazardous midnight-to-dawn sector. The total fluxes at 0900 were quite 
close (within a fraction of 2) of charging the spacecraft in sunlight. 
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The first set of spectra are taken in the quiescent plasma at 0700. 
The next spectra taken at 0900 show the first encounter with this event. 

The low-energy spike seen in both electron detectors is not due, as might be 
suspected, to charging positively, but rather is most likely the locally 
produced secondaries being reflected from a suddenly enhanced plasma sheath 
about the spacecraft (see discussion by Whipple, 1976 ). This event could 
easily have produced charging in excess of 10,000 volts if the spacecraft 
had been eclipsed at this time. Such a sharp, localized event was 
probably responsible for the main power supply failing on a non-NASA space- 
craft. (Note both ATS-5 and ATS-6 are research craft and are somewhat 
better built than operational craft. Therefore one is not surprised that 
they can sail through disturbances that would sink weaker vehicles). 

Spectra at 1000 and 1100 document the evolution of the event in the 
normal manner. The predominant spectral changes are again caused by 
gradient drifts. 

Spectra at 1400 and 1700 show the complex spectra that can result 
from the combination of injected particles, particle losses (the chasm in the 
protons, and loss of high-energy electrons), and multiple encounters with 
the high-energy particles as they circle the earth. These spectra would be 
useful for studying spacecraft conditions near noon. However, a low- 
temperature plasma must be assumed to be present in both sets (see previous 
section.) (See also the descripton of day 12/3/70 for a similar event.) 

3.3 QUIET DAY 

2/12/70 This day was chosen to illustrate a quiet period partly 
because it is one of the longest quiet times normally seen, and because it 

fortunately followed the intense event already described on 2/11/70. 

Therefore we have a single two-day period of uncommon interest for this 



project. Note that we have small data gaps at both the start and end of the 
day. These are of no consequence since the activity is so low. 

The spectra for this day is simply spaced throughout to sample 
uniformly. Any use of these spectra must assume the presence of copious cold 
(or "warm") plasma with densities of at least 30/cm. 

3.4 POST-MIDNIGHT SUBSTORM 

3/18/70 The event shown on this day actually consists of two 
closely spaced injections occurring after a large quiet period. The activity 
starts at about 1040 VI. This is sufficiently past midnight that the plasma 
response is very different than the case shown for 2/11/70. In particular, 
the main body of protons do not reach the spacecraft until about 1330 after 
travelling around the world to the west. This situation could produce 
hazardous charging. 

The first set of spectra taken at 0900 set the stage for the later 
injection. The next three sets are spaced somewhat closer than the nominal 
minimum of one hour followed in the rest of this report. This was 
necessary because of the rapid development and the desire to show all phases 
of the event. The set at 1000 shows some electron enhancement over the 
earlier spectra. The 1050 spectra shows the leading edge of the injection. 
[Note the apparent oscillations in the electron spectra in this and the 
following set are an unavoidable artifact due to the particular operating 
mode chosen that day and the relatively rapid changes taking place]. At 
1112 we see significant changes in both the low-energy electrons and the 
shape of the protons. But at 1200 we see even hotter electrons instead of 
the effective cooling we would expect normally. The explanation is easily 
seen in the spectrogram: another injection has followed this first. This 

is common and does not affect the usefulness of this day for the report. 

The final set of spectra follow the injection development. Only 
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now we see that the electrons have experienced rapid depletion (by probably 
precipitating into the atmosphere) and the arrival of the protons from 
their trip around the world has given us a spike in the distribution. 

3.5 PRE-MIDNIGHT SUBSTORM 

12/3/70 The main feature of this day is the surplus of high-energy 
protons early in the day. Although this condition is probably not hazardous 
to spacecraft from the point of view of charging, it is a common occurrence, 
and the vehicle's response should be studied. Spectra are provided for 0200, 
0400, 0600, 0700, 0800, and 1200 LT. The first four are of prime interest 
for the study of the effect of high-energy protons. The 0800 spectra can 
be used in a way similar to those of 2/11/70 for intense localized substorms. 
The main difference between the two being the higher energies, but lower 
fluxes seen on 12/3/70. The last set of spectra (1200 LT) are provided 
simply to show the time development of the storm. 

3.6 ECLIPSE AND SUNLIGHT CHARGING 

3/14/71 Although the intent of this report is to provide isolated 
examples of various types of events at synchronous altitude, we realize 
that for many purposes this is not sufficient. Therefore we also present 
an active day which has both charging in eclipse and a good example of 
charging in sunlight. As can be seen from the spectrogram, this day is 
very different from the other examples. Several distinct substorms 
follow one after the other. The plasma conditions change so quickly that 
obtaining good averages is difficult. The charging events are easily 
identified by the bright bands in the low energy protons. The eclipse is 
always centered about local midnight, and the sunlight charging on ATS-5 
is always observed between midnight and dawn. 
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The first spectra taken at 0600 is pre-eclipse. The next two 
are at different phases of the eclipse. The 0800 set in post-eclipse. The 
last two sets of spectra are preceding and during sunlight charging. 

The parallel electrons in the last two cases show the effects of 
differential charging. 

Persons using this set of data might want to correct the fluxes to 
what they would be if there were no charging. The cookbook method for 
doing this assumes that the instrument is a differential detector. Then 
by Liouvi lie's theorem, 

Jp(E) * [E 2 /(E - q*) 2 ] J m (E - q*) 

where J p is the predicted flux at energy E, J m is the measured flux, and 
4> is the potential. The sign of the charge, q, is positive for ions and 
negative for electrons. 

For the one sunlight charge case shown, <j> is -80 volts. Therefore 
corrections above a few thousands of volts are unnecessary. For the 
eclipse cases essentially all channels should be corrected. 

The lowest energy fluxes of electrons and the highest energy ions 
will be missing from the corrected spectra. 
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APPENDIX A 


* 


DESCRIPTION OF ATS -5 SPECTROGRAMS 


Format 

The spectrograms are produced in pairs: one showing the spectra 

from the perpendicular proton and the perpendicular electron analyzers 
and one showing the spectra from the parallel proton and electron 
analyzers. They are labeled by a large ]_ or J | on the middle left side. 
The proton part is always below the electron part. The day of the year 
(January 1 equals day 1) and year is given at the bottom. The month, 
day in month, and year are also given at the left just above the or | | 
label. The times at the beginnings and ends of the spectrograms can be 
arbitrarily set, and can cover any desired time span. Time scales cover- 
ing as little as 10 minutes and as great as 4 days have been used. When 
more than one day is encompassed, either negative hours or hours greater 
than 24 are used to prevent any ambiguity. Grey scales are located at 
the right. Six different integrals are plotted in grey coded bands in 
the upper part along with magnetic field quantities. At the very top 
are two data quality indicators. 


Grey Scale Interpretation 

The primary value of spectrograms is their ability to reveal, pat- 
terns in the energy-time plane. The determination of actual flux levels 
from them is of secondary importance. For this reason, and because of the 
loss in time re solution, the option which produces a coded pattern with 
which accurate flux values can be obtained is now rarely used. Color 
coding also permits accurate values to be obtained, but is more 
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expensive than grey coding. In the present case, color is reserved for 
adding another dimension: by superimposing the perpendicular and 

parallel spectrograms with color filters limiting each to one -half of 
the visible spectrum, the energy and time dependence of the pitch angle 
anisotropies are clearly displayed as patterns of different shades of 
color. 

Should one desire to estimate the flux at a given point on a 
spectrogram, first locate the corresponding level on the grey scale 
at the lower right and determine the value of "g" on the scale marked 

p 

0 to 3- The differential energy flux in ev/cm sec sr ev is then given 
by 

(icP - 1) 10 b + k ’ 367 

where b is given by "EL" in the lower left comer of the spectrogram 
for the electron fluxes or "PR" for the proton fluxes. The value of 
"ST" in the lower left corner gives the change in G between each of the 
33 discrete grey levels available. 

One option available is to let the grey scale recycle repeatedly 
instead of simply saturating. This option with a small value of "ST" 
is used to reveal small variations over a wide dynamic range of fluxes. 

Energy Scales 

The computer program which generates the spectrograms can utilize 
any arbitrary function of energy for the energy scales for exhibiting 
all or any part of the measured spectra. The entire range from 50 ev to 
50 kev is usually plotted with one of the two types of scales: 
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logarithmic with 50 ev at the bottom for both protons 
and electrons. 

proportional to l/(E +3 kev) with the electron part 
inverted and sharing the same point with the protons at 
zero energy. The bias of 3 kev was arbitrarily chosen 
to give a good presentation of the 50 ev to 50 kev energy 
range. If the scale, S, is taken to be 0.0 at infinite 
proton energy, 1.0 at zero electron and proton energy 
and 2.0 at infinite electron energy, then 

S = 3 ' kev 6V wliere ® * s t ^ ie energy in kev 

q = ± 1 depending on the sign of the particle's charge. 
Note that at low energies, S ~1 + qE/3 kev. Time tic 
marks are located at S = 0 , 1 , and 2 . The extrapolation 
of dispersion curves back to the time marks (at S s 0 or 2 
yields the time infinite energy particles would have ar- 
rived, and therefore, the time of the event responsible 
for the dispersing particles. The slopes of the high 
energy parts of dispersion curves give a measure of the 
distance of the satellite from the regions in which the 
particles were perturbed, but it is apparently necessary 
to include electric field effects to obtain useful 


accuracy 



Subsidiary Data 

A number of useful quantities are given in the lower left hand 

corner. 

The analyzers in the "master" and "mate" channels are identified 
by numbers following "MASTR" and "MATE" according to the scheme: 

1. perpendicular electron analyzer 

2. perpendicular proton analyzer 

3. parallel electron analyzer 
b. parallel proton analyzer 

TA = averaging time for the spectra in minutes. 

TS = time between spectral averages in minutes. 

TM = averaging time for the magnetic data in minutes. 

The seven bit command word is given immediately below "COMMAND". 
The first three bits give the channel assignments and are therefore 
redundant to the master and mate identifications given above. Bits 4 and 
5 specify the operating mode according to the scheme: 
bit b 5 Mode 


0 

0 

1 

1 


0 

1 

0 

1 


track- scan 

single step scan only 
track only 

double step scan only 


Bits 6 and 7 not set to zero correspond to other modes which are rarely 
used. 

"ST", "EL", and "PR" are described above. 

"PSNG" specifies the quantity being plotted in the spectrogram 
according to the scheme: 
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1. differential energy flux 

2. differential number flux 

3. ratios of the flux averaged over "TS" minutes to 

the flux averaged over the previous "TA" - "TS" 
minutes . 

U. ratios of adjacent energy steps. 

Options other than the first are used only in special studies. 

If the option to make the background black rather than white has been 
used, then PSNG will be negative. A black background is preferred for 
slides that are to be projected. 

Magnetic Field 

Data from the ATS -5 magnetometer have been kindly supplied by 
T. Ski ll man of the Goddard Space Flight Center and are plotted above the 
spectral data along with lines at 0, 50, 100 and 150 gammas. The data 
are not corrected for the effects of time changes in the spacecraft cur- 
rent systems. These perturbations can be as large as 15 gammas. The 
absolute value of the magnetic field component parallel and perpendicular 
to the spin axis is given by the darker and lighter points respectively 
(and usually the upper and lower respectively) with the spectrograms of 
the perpendic ul ar analyzers • The perpendicular component is obtained 
using only the coarse (33 gamma step size) data and is thus uncertain 
by at least ± 10 gammas. Most of the scatter in this component is due to 
using only the coarse data. 

The ma g n itude of the field and the angle of the field to the spin 
axis are given by the lighter and darker points respectively (and usually 
the upper and lower respectively) with the spectrograms of the parallel 



analyzers . The angle to the spin axis is given in degrees . Both the 
magnitude and angle are subject to the additional uncertainties in the 
perpendicular component. 

Integrals 

Above and below the magnetometer data are six strips in which 
various quantities are logarithmically encoded in a grey scale such 
that a ratio of about 2000 to 1 is covered in going from black to white. 

In the 1st, 2nd, 3rd, and 5th strips, the following integrals 
from the perpendicular and parallel analyzers are plotted with perpendicu- 
lar and parallel spectra respectively: 


Label 

Quantity 

Value at Midpoint of Grey Scale 

PR N DEN 

proton number density 

o 

1.0 proton/cnr 

EL N DEN 

electron number density 

1.0 electron/cnr 

E E FLX 

electron energy flux 

2 

1.0 erg/cm sec sr 

PR E FLX 

proton energy flux 

2 

1.0 erg/cm sec sr 


In the 4th strip labeled "PRESSURE", the total perpendicular 

electron plus proton pressure is plotted with the spectrogram of the 

—8 2 

perpendicular detectors with a midpoint value of 10 dynes/cm . In the 

4th strip with the parallel data, the magnetic field pressure is plotted 

_8 2 

with a midpoint value of 2 x 10 dynes/cm . 

In the 6th strip (near the top) labeled "PAR NFLX" the parallel 

electron number flux is plotted with the spectrogram of the perpendicular 

8 2 

detectors with a midpoint value of 10 electron/cm sec sr. In the top 

strip with the parallel data, the parallel proton number flux is plotted 

7 2 

with a midpoint value of 10 protons/cm sec sr. 
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Data Quality Indicators 


At the very top of the spectrogram is a line which increases in 
breadth with an increasing percentage of missing data. In the track-scan 
mode, about 73 percent of the potential data is usually "missing" since 
75 percent of the time is spent tracking a peak in a narrow spectral 
region. When data are not available, previous data are used unless the 
time gap is greater than 30 minutes in which case the spectrograms are 
left blank. The top line, of course, goes to its maximum width during 
gaps in the data. The magnetometer data is not plotted during such gaps. 
Care must be exercised to avoid false interpretations of spectrograms 
containing data padded in from an earlier time. 

Just below the missing data line is a line which becomes darker 
and thicker with increasing numbers of bad points. Often the quality of 
data transmission is such that over one percent of the data points are 
bad. Even the highest quality data being obtained are usually incorrect 
more than 0.1 percent of the time. This corresponds to over 800 bad 
data points per day of data. A data editing scheme has been devised 
which eliminates approximately 99 percent of the bad data and rarely 
removes data later Judged to be good. Failure to remove bad points 
usually occurs when the false data happen to form a self-consistent 
context. This type of failure to edit properly is responsible for the 
two white areas in the lower right of Figure 4. The bad data indicating 
line reaches its maximum thickness when there are more than 10 bad 
points in the four spectra measured during the time covered between 
averages (equal to "TS"). 



APPENDIX B 


DESCRIPTION OF ATS -5 SPECTRAL AVERAGE PLOTS 


Format 

The spectra from the two electron and the two proton analyzers 
are plotted in adjoining log-log plots with borders at 30 ev ©Jid 100 kev. 
The range of the vertical scale is variable and depends upon whether the 
differential energy flux or the differential number flux is being plotted. 
The parallel electron spectrum is shifted down by a factor of 100 (i.e. 
x O.Ol) and the perpendicular proton spectrum is shifted up by a factor 
of 100 (i.e. x 100). These shifts usually provide adequate separation 
and place the perpendicular spectra above the parallel spectra in each 
case . 

The universal time at the midpoint of the data being averaged 
over is given twice at the top of ^he plots. On the left hand (electron) 
side, the time is given in hours, minutes, tenths of minute, month, day 
of month, and year, and is followed by the averaging time in minutes. On 
the right hand (proton) side, the time is given in hours (to the nearest 
one thousandth of an hour), day of year (January 1 equals day l), and the 
year. The local t ime in hours and minutes is sometimes added on the left 
side. 

Also given near the top are four different integrals over each 
of the four spectra. The integrals for the perpendicular data are given 
above the integrals for the parallel data. Following two of these sets 
of integrals will be found the words "MASTER" and "MATE" to indicate 
which analyzers are occupying the two non-subcommutated data channels . 

• 
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When in the track mode, the "master" analyzer controls the peak tracking 
system. The operating mode (for example the scan only or track-scan 
modes) of the system is given on the right side. 

Error Bars 

Vertical bars which encompass the middle 68.26 percent of the 
Poisson distribution are given at each data point. At high rates, they 
correspond to plus and minus one standard deviation. The approximation 
N ± = N ± JT (1.0 - 0.17/n) is used where N is the total n umb er of 
counts accumulated at the point. 

When in the track-scan mode, there are about four times the n umb er 
of accu mu lations at the points near the energy of the peak being tracked 
than at other energies. Also the spectra from the "Master and "Mate" 
channels will have about twice the accumulation time as the other two 
( subc ommutat ed ) spectra. 

When in the single step scan only mode, every other data point 
in the subcosmutated spectra will be missing. This under-sampling of 
the spectra can lead to substantial errors in the smooth line drawn 
through the data points since structure as sharp as the instr ume nts ' 
resolution is frequently observed. 

If zero counts are obtained, then the error bar is replaced by 
a triangle pointing up to the line which is placed at one -half the 
flux corresponding to one count being accumulated. 

If no data are available for a point during the t ime period 
being averaged over, then the flux obtained during a preceding t ime 
period is inserted. In this case, the error bar is re plac ed by a 
triangle pointing down to the data point. 



Integrals Over the Spectra 

The four integrals given for each analyzer at the top of the 
plots are of course intrinsically directional quantities. The parallel 
cases correspond to pitch angles a 0 (a = the angle of the spin vector 
to the magnetic field vector) and the perpendicular cases correspond to 
averages over the pitch angle range of 90 ± « degrees. The integrals are 
taken only over the measured range of 50 ev to 50 kev and are, therefore, 

lower limits. 

The n umb er densities in particles/cm 3 are labeled "DEN" and 
correspond to 4 rr times the directional number densities in 

O 

paxticles/cnr sr* 

The particle pressures in 10' 9 dynes/cm 2 are labeled "PRES". 

They correspond to 8 tt/3 times the directional energy densities in 

ergs/cm 3 sr. The multiplication by 8 tt/ 3 simplifies computation of the 

total particle pressure perpendicular to the magnetic field vector. 

2 

The directional energy fluxes in ergs/cm sec sr are labeled 

"E FLX". 

6 2 

The directional number fluxes in 10 particles/cm sec sr are 
labeled "N FLX". 
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APPENDIX C 


Construction of Complete Spectra 

The data presented in the main report can be combined with 
experience gained in the ATS-6 program to construct a most probable set 
of total spectra. This consists of adding other components to the 
measured fluxes. 



Let dN = number density between E and E + dE, ft and ft + dft, then 
for the magnetospheric plasma at geosynchronous orbit one has 


dN = dN , . + dN . + dN*. , . . . . 

cold iso field aligned 


where 


1) COLD 


M 


- A & e 


"r 4nc~ o ~ ) 


cm ait' 

for 0 < e < 50 ev 




2) ISOTROPIC 

dKso = ff^dj 

for 50 ev < e < 50 kev 


where dj = energy flux /an -sec- s ter * data 


3) FIELD ALIGNED 


.,/ ./ / /h -[te-vmfccos* 

cINf.A ~ Nr.fi. \zTrK~fc*J ^ ^ 

for 0 < e < « , a = pitch angle 


The total number density 


t! ~ /JM — hfcoL.0 t Afi/fre+«At. d hfp^^ 

Charge neutrality requires that 
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The temperature chosen for T cold should probably be a few volts to a 
few tens of volts. The density of the cold component can be estimated from 
Reasoner's work (1975) and figure 2.3 to be about 30/cm . 

The form for the field-aligned component was derived from the 
assumption of a displaced maxwell ian plasma falling through a potential well 
of e Q . If we assume that these particles have their origins in the 
ionosphere, then we can estimate 100 <e Q < 10,000 electron volts and that 
kTp ^ is a few electron volts. The density is more difficult to estimate, 
but a few percent of the ambient would be consistant with measurements. 

Note that the field-aligned component is probably only important 
for the study of differential charging since it only influences the charge 
state at locations where the bulk of the plasma is excluded ( i . e . , in 
properly oriented cavities on the vehicle). 
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Spectrogram for 3/14/71 - Eclipse and sunlight charging. 
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Spectrogram for 2/11/70 - Intense localized post midnight substorm. 
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1.0 Introduction 


The purpose of this report is to develop plans for the 
experiments which must be conducted in order to verify the 
ground test mathematical model (GTMM) of spacecraft charging. 
Experiments will determine the ability of computer codes, 
which are being developed, to predict the electrostatic fields 
and charge distributions in the region around the GTMM. The 
verification tests consist of a matrix of experiments. 

The first experiments are simple. Later experiments are more 
complex culminating in a full scale test of an operational 
satellite. 

The philosophy of this report is to establish the overall 
objectives in depth and then to explore the details of 
implementation. While the overall objectives are well- 
defined, the physical implementation as presented is flexible 
enough to allow a variety of engineering compromises and 
optimization. 

Often during the development of this report a physical 
device or instrument was needed which was either nor available 
commercially or had never been designed. Examples are the 
Distributed Source Accelerators (Discussed in Appendices 2 
and 3 ) and the rotating sensor electric field mill illustrated 
in Figure 3 . 6 . When this problem occurred our approach was 
to provide a rough conceptual design of the needed device and 
then indicate an estimated level of effort necessary to do 
detailed design and to construct the device. While the 


conceptual designs are often sketchy they do provide a direction 
and provide a better response to the problem at hand then 
simply saying "There is no such device." 

Facilities for the performance of these tests are 
available at NASA Lewis Research Center. The 15 ft. di am eter 
tank described by Finke, et al. 1 is specifically designed 
for the testing of space packages and will, when adequately 
instrumented, make a superb facility for performing the 
FTMM verification tests. Preliminary tests of smaller test 
bodies can be performed in the 6 ft. x 6 ft. test tank which 
already contains substantial instrumentation (See Berkopec, 
et. al. ^). 

Section 2.0 contains a general description of the overall 
experimental strategy, and section 3.0 contains the description 
of specific experiments, including test set ups, procedures, 
and quantities measured. 

In summary, the experimental plan presents a series of 
experiments, each repeating a prescribed procedure to measure 
an important physical parameter, which will enable the 
predictive ability of the computer codes to be determined. 

The ground test mathematical model (GTMM) will thus be verified 
b y quantitative experiments. In addition, several extremely 
important spinoffs will result. The test facility which is 
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instrumented for an adequate GTMM verification will also be 
ideally instrumented (and calibrated) for a full scale 
operational satellite test. Tests of electrostatic properties 
can be performed. A new level of environmental simulation 
will be attained. 



2.0 Experimental Strategy 

The experimental strategy developed in this plan is as follows. 
Begin with as simple a situation as possible and attempt to 
understand it. Then in systematic steps increase the complexity 
of the situation. In this way, one can hope to progress from, 
xor instance, a flat plate made of a single conductor irradiated 
by monoenergetic electrons toward a full scale test of an 
operational satellite being charged by distributed energy beams 
of electrons, ions and photons. The plan develops a matrix of 
experiments where complexity varies along three dimensions 

representing geometry, materials and irradiating sources. 

(See Figure 2.1). 

For each experiment in the matrix the same instrumentation is 
employed, according to prescribed procedures, to measure the 
same physical parameters. Thus, as an example, a scan of an 
electrostatic voltmeter from the outer boundary to the surface 
of the test body along a prescribed path is a procedure which 
would be repeated for a number of different test bodies, each 
representing a different point in the experimental matrix. The 
data generated is then compared with that calculated using the 
computer codes, thus determining the predictive ability of the 
codes. A serious attempt has been made in this plan to develop 
experimental procedures which are a sensitive test of the 
computer code capabilities. 
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Fiaure 2.1. EXPERIMENTAL MATRIX 


All tests are assigned an experimental code (S, G, M) 
which defines the tests position within the experimental 
matrix. Large values of S, G, or K signify more complex 
experimental situations. 


[•] 


A flow diagram indicating the overall experimental approach 
is shown in Figure 2.2. One begins by specifying a point 
in the experimental matrix, i.e. one must define the three 
parameters S. G. M. which specify the particle source, the 
geometry and the materials of the test object, respectively. 
Then one runs the computer codes and performs the experiments. 
Next the results are compared. We note that these experiments 
generate a large amount of data, and that computer aided, 
on line, interactive reduction and presentation of this data 
should be seriously considered. The data output should, 
whenever possible, be formatted in such a way as to enable 
direct comparison with the output of the computer codes. 

After comparison of the data, it will likely be necessary 
to modify either the experiment or the computer code. One 
then repeats, as necessary, the process until the results of 
experimental and computer code agree. 

Finally, one selects a new point in the experimental matrix 
and repeats the entire procedure. 
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Figure 2.2. 

Experimental approach to computer code verification 
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As implied by the center box in Figure 2.2, a detailed 
examination of the comparison procedure is necessary. The 
comparison should take into account the structure of the 
computer codes and the difficulties of the experimental 
situation. Points to where the code predictions are sensitive 
and which can be experimentally examined with ease are sought 
out. Results of this examination follow as we present the 
details of the comparison procedure. 

A review of the computer code development as presented in the 
Systems, Science and Software interim report (SSS-R- 77-3124) 
provides the flow diagram shown in Figure 2.3. The code 
development has proceded in "top down" fashion with each 
block being filled in with increasing detail. The experiments 
are similarly developed in a top down fashion with the plan 
structure paralleling the code structure. Both the computer 
code and the experiments accept as input data the specifications 
of the test object and particle source characteristics. 

Thus to specify the input data to either the computer code or 
the code verification plan one must specify a point (S, G, M) 
in the experimental matrix. Verification consists of comparing 
the results of experiment and computer code then modifying 
one or the other until any differences are resolved. 
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Figure 2.3. Flow diagram of the main program. 
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xf possible, individual subroutines should be verified 
individually by experiments . The computer codes under 
consideration lend themselves well to this approach. 

The subroutine POTEITT, which calculates the electrostatic 
potential that results from a given geometric configuration 
of boundaries and charges, can be verified individually by 
experiments aimed at measuring potentials or alternatively 
capacitances. Figure 2.4a illustrates the approach. Trajectory 
measurements will provide a test of the capability of the 
subroutine PUSHER as indicated in Figure 2.4b. The verification 
of POTENT and PUSHER should be relatively straightforward. 

MATRIL, the subroutine which treats charged particle 
interactions with the surface will be the most difficult to 
verify. The source of this difficulty lies in the difficulty 
of performing the experiments and in the wide range of variability 
of material surface properties. The approach is indicated 
in Figure 2.4c. An important result which should result from 
the implementation of this experimental plan is an improvement 
in the state of the art of materials measurements. The 
combination of computer codes with experiments will, once 

verified, create a powerful new tool for the study of material 
surface properties. 
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1 

Proceed to SELECT 


Figure 2.4a. 

Testing of subroutine POTENT. Potentials or alternatively capacitances 
can be measured and compared with the results predicted by the code. 
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Figure 2.4b. 

Testing of subroutine PUSHER. Charged particle trajectories in fields which are 
known from either measurement or calculation can be experimentally determined 
Comparison with PUSHER trajectory calculations leads to verification? 
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Fiaure 2.4c. 

Testing of subroutine MATRIL. PUSHER is used to compute trajectories for 
materials experiments. Results are compared to verify MATRIL. 
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After verification of the individual subroutines the experiments 
should focus upon the verification of the overall code. (See 
Figure 2.2.) As a diagnostic procedure one may revert at any 
point during the overall verification of the GTMM to the 
procedures for testing an individual subroutine. 

In summary, the experimental plan presents a series of 
experiments, each repeating a prescribed procedure to 
measure an important physical parameter, which will enable the 
predictive ability of the computer codes to be determined. 

The ground test mathematical model (GTMM) will thus be 
verified by quantitative experiments. 
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3-0 Ground Test Set UP 


The purpose of this section is to describe a test facility 
which will be used to verify the Ground Test Mathematical 
Model (GTMM). Specifications for this test facility should 
determine the electrostatic boundary conditions and the source 
particle characteristics. To test the GTMM it is necessary 
to measure the electric charge distribution, the electric 
field intensity, the electric potential, the total current to 
the test object, (and leakage current) as well as current fluxes 
throughout the test space. These quantities are then to be 
compared with the results predicted by the GTMM. 

In Figure 3.1 we indicate schematically the basic elements 
in an experimental facility designed to assess the GTMM. The 
facility supplies ion, electron and photon sources capable of 
simulating magnetospheric fluxes, establishes the appropriate 
electrostatic boundary conditions and provides the instrumentation 
required to measure the quantities stated in the previous 
paragraph. Data processing and display is incorporated into 
the facility in order to enable the large amounts of data 
which will be produced to be examined intelligently. 
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3.1 Test Facilities 


Tests will be performed in two facilities which exist at 
NASA Lewis Research Center. The smaller of the two, a tank 
which is 6 ft. diameter by 6 ft. long, is described in 
NASA TM X -73602 by Berkopec, et al.^, is instrumented for 
use as a substorm simulation facility. This facility has 
been used by Stevens, et al. to measure the response of small 
samples of material to fluxes of electrons and photons 
comparable to those found at Geosynchronous Orbit (GSO). A 
larger facility measuring 15 ft. diameter by 63 ft. long 
is described in NASA TM D- 277 1 * by Finke, et al. 1 This tank 
is designed primarily for environmental testing of space 
packages and plasma thrusters. It is large enough to perform 
full scale tests of an operational satellite. In this report 
we outline an instrument complement which, when installed in 
the large tank, will enable the experimenter to obtain the 
data necessary to verify a ground test mathematical model (GTMM) 
and to perform a full scale spacecraft charging test on an 
operational satellite. 

The large LeRC tank is capable of high pumping speeds and 
ultimate vacuums approaching 10“ 8 Torr thus it is suitable 
for use in this work. The high pumping speeds (e.g. Atmospheric 
to 10 6 Torr in about 2 hours.) are important during the 
initial setting up of a test, when it may be necessary to 
pump down from atmospheric several times in order to work out 
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the details of a test set up. The ultimate high vacuums are 
important m the simulation of the hard vacuum of space. We 
note that a pressure of 4 x 10 7 Torr is obtainable in the 
empty tank without coolant in the pump traps, thus early 

tests could be run without using liquid nitrogen at a potential 
cost savings. 
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3.2 Boundary Conditions 

Lf :> 

An earlier MAYA report (Appendix E of the Interim Report 
3SS-R- 77-31 2*0 specified boundary conditions and source 
particle characteristics for the ground test environment. An 
amplified version of this work, targeted for use in the large 
LeRC tank is presented here. 

It is of great importance that the electrostatic boundary 
conditions be accurately defined in the ground test environment. 
If effort is not expended to do so, then stray fields can 
invalidate the test data. In order to obtain an accurate 
definition of the electrostatic potential at a large distance 
from the test object a Faraday cage should be constructed 
to contain a test volume within the vacuum tank. One might 
consider as an alternative simply allowing the vacuum tank 
to act as the Faraday cage. This alternative should be 
rejected. Despite the simplicity of implementation, using 
the vacuum tank as a boundary results in a geometrically 
and electrostatically irregular surface which will complicate 
test interpretation. Surrounding the test object with a 
Faraday cage will enable the experimenter to specify accurately 
the outer boundary condition for the electrostatic potential. 

The satellite model or test object is suspended in the 
center of the cage. The supporting strut requires some care 
in its design in order to avoid seriously perturbing the 



potentials. In Figure 3.2 we show an approach to the design. 
The strut is made of a number (say 10) of conducting 
cylinders, each of which is held at well-defined potential, 
chosen so as to minimize the perturbation caused by the strut. 
An example of the choice of strut potentials for a cylindrical 
test object at vacuum potential is given in Appendix 1. 
Voltages, currents and signals to and from the test object 
are conducted on wires contained within a shielding cable 
which runs through the strut. 

A cross section through a typical test object is shown in 
Figure 3-3. The test object consists of two parts, an outer 
shell which represents the satellite surface and an instrument 
module contained within a Faraday cage. This Faraday cage 
establishes a solid reference potential and shields the outer 
shell from extraneous fields produced by the instrumentation. 
The outer shell # can be changed to test different satellite 
configurations. The instrumentation module will be used with 
the different outer shells to establish potentials on the 
surface of the satellite and to operate those diagnostics 
which are located on the satellite model. 

The inner and outer boundary conditions on the electrostatic 
potential are thus well defined by the experimental set up. 
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3-3 Particle Source Characteristics 

The plasmas present in the magnetosphere at geosynchronous 
orbit are too hot and tenuous to be produced directly in a 
ground test environment. (See Appendix for a comparison of 
magnetospheric and ground test environment.) While it is 
possible to produce plasmas with the requisite number density 
(n -v 1 -100/cm 3 ) it does not appear possible to heat these 
plasmas to kilovolt temperatures. We have therefore elected 
to simulate the magnetospheric environment with particle 
accelerators. The accelerators provide a directed, rather 
than a random, flux of energetic particles but have the 
advantages of ease of design, operation and characterization. 
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3 - 3-1 Accelerator Characteristics 

Hie particle accelerators chosen for use in the ground test 
environment should be capable of producing fluxes which are 
comparable in intensity and in energy spectrum to the fluxes 
present in the magnetospheric environment. Data representing 
this environment have been previously presented in the MAYA 

report entitled "A Preliminary Specification of the Environment 
at Geosynchronous Orbit" (SSS-R- 76-2996) . ^ 

The electron data given in this report are characterized by 
a high intensity number flux re % 1 . 5x10^ o/ cm 2 _ sec with 
mean electron energy of 4-5 Kev seen on hour 10.8 of 3/18/70. 
This situation is typical of a post-midnight substorm and could 
produce hazardous charging. High proton fluxes occur on the 
day 2 /11/70 during hour 9. These fluxes, ri * 2. 5x10 8 /cm 2 -sec 
are characterized by a mean particle energy greater than 
10 Kev. The particle accelerators should be capable of 

producing fluxes with these properties over some large part 
of the test chamber* 

With this background, after extensive review of the data, we 
have chosen parameters for the particle accelerators which 
will bracket the magnet ospheric variations. These parameter 
choices are shown in Table 3.1. 
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Parameter 

Directed Flpc Intensity 
(cm* 2 sec" 1 ) 

Current Density 
(amps/cm 2 ) 

Particle Energy (Kev) 

Energy Spread 

(Full Width at Half 
Maximum) 

Table 3 . 1 


Electron Beam 
5x10 6 <re<5x10 1 0 

SxlO" 1 3 <Je<8x10" 9 

1 <V <40 
e 

10 


Ion Beam 
1x10 6 <ri<1x10 1 0 

1 .6x10' 13 <Ji<1 .6x10 

1 <v i < 1 +o 

10 $< ^1 < 50 ^ 

v i 

Particle Accelerators 


Characteristic Parameters for 


flux densities are assumed to be delivered to an area of 
approximately one square meter. Thus the maximum total 
currents are approximately 8 x 10" 5 amps and 1.6 x 10" 5 amps 
for the electrons and ions respectively. • The accelerators 
to be used should be low current large area devices. The 
beam area at the exit aperture of the gun should be at least 
0 . 05 m 2 , and preferably larger, so that an area of lm 2 can be 
covered at a distance of approximately 1 meter from the gun 
without excessive beam divergence. A particle source of this 
size could be tested and debugged in the small 6 ft. vacuum 
tank with existing instrumentation. Note that the large 
vacuum tank is 15 ft. in diameter, hence has an area of 16.7 m 2 
so that multiple particle sources must be considered. 

The energy spread indicated is optional but obviously 
desirable. The accelerators indicated for this work are of 
a rather peculiar sort. One usually seeks to make an 
accelerator with as monoenergetic a beam 



as possible and often one seeks high current densities. For 
this application we are operating in the opposite limit, 
seeking a beam broadly distributed in energy with low current 
densities. Such beams are not currently commercially available, 
nor have they to the best of our knowledge designed or built 
in other research programs. We have thus undertaken to outline 
approaches to the design of such accelerators. In Appendix 2 
"Distributed Source Ion Accelerator" we provide simple 
calculations which indicate that an ion accelerator based 
upon the use of crossed electron and neutral beams would 
provide the requisite currents. Such an accelerator has 
never been built and would require perhaps two man years 
to develop. A "Distributed Source Electron Accelerator" 
is discussed in Appendix 3* Again such a device has never 
been built, although MAYA personnel have used photoelectron 
source mono-energetic accelerators for the calibration and 
testing of auroral particle experiments on Advanced Technology 
Satellites 5 and 6. A distributed source electron accelerator 
would require more development, perhaps on the order of one 
man year by a qualified scientist. 

3*3*2 Experimental Determination of the Source Particle 
Characteristics 

The stability and accuracy of the particle accelerators may 
well be good but should not be trusted. Actual source particle 
characteristics should be measured experimentally, by examining 
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the bean: with a retarding potential analyzer (RPA). The 
particle energy, variation of energy and current density across 
the beam and the angular spread of the beam should be defined. 

A series of scans with an RPA across the beam at different 
distances from the source will provide the requisite information. 
Using a MAYA RPA (See Section 3*5*4) energy analysis can be 
performed with an accuracy of at least 1 % of the total beam 
energy and a spatial resolution of 0.1 cm. 

3*3*3 Programming of the Source Particle Characteristics 

Substorm simulation will be achieved by programming the 
temporal variation of the source intensity, mean energy and 
energy spread. This programming which should be controlled 
by the developing computer data link can be based upon the 
integrals computed as 2.3 minute averages from the ATS- 5 data 
presented in the report SSS-R-76-2996. ^ 

3*3*4 Complications Caused by Electron Cyclotron Radius 

As has been noted in the Introduction, the electron cyclotron 
radius in the terrestrial experiments can be comparable to 
the dimensions of the test object. (r '2.5m for a 1 Kev 

V 

particle in a magnetic field of 0.3 gauss, See Figure 3.4a . 

Thus as a minimum provision should be made for pointing the 
electron accelerator through a range of angles. 

For a monoenergetic accelerator the beam may be deflected 



electrostatically to compensate for the deflection by the 
magnetic field. Alternatively, the accelerator may be 
attached to a gimballed mount and thus pointed mechanically 
to compensate for the magnetic deflection. 

An accelerator producing a beam distributed in energy will 
exhibit worse problems. Electrons with different energies 
will be deflected by different amounts. Neither electrostatic 
deflection nor mechanical pointing of a gimballed accelerator 

will compensate for this effect. Figure 3.4b illustrates 
the problem. 

Photoelectrons are characterized by electron cyclotron radii 
of perhaps 10 cm. (See Figure 3 . 4. a.) The photoelectric 
sheath around a test object much larger than 10 cm would thus 
be substantially modified by the presence of the earth's field. 
'While the computer codes can handle the effect of the magnetic 
field on the particle trajectories, the physics of the 
sheath may be changed profoundly. Some numerical and analytical 
work is needed to determine the impact of finite electron, 
cyclotron radii. Serious thought must therefore be given to 
constructing a Helmholtz coil system for the purposes of 
reducing the field within the test chamber. 
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3*3-5 Photon Source 

Photons should be provided by a solar simulator. Since 
photoemission is very important in the determination of 
overall satellite potentials both the intensity and spectrum 
of the illuminating photons should be carefully matched to 
that of the sun. ORIEL Corporation of America, 15 Market St., 
Stamford, Conn. 06902 is a comprehensive source of solar 
simulation equipment. 



3 A Measured Quantities 


The computer codes use as input the following: 

1 ) specification of the boundary conditions far from 
the satellite 

2) specification of the source particle characteristics 
far from the satellite 

3) detailed geometry and materials of the satellite 

The computer codes then calculate the following: 

1) particle trajectories, interactions at the satellite 
surface, and resulting electric charge distributions 

2) electric field intensity 

3) electric potential 

4) total current to the body and total leakages through 
the body 

5) current fluxes 

This experimental plan develops the methods of measuring 
items 1 ) to 5) so that direct comparison can be made between 
experimental measurements and the results of the computer codes. 
To do so we define an array of diagnostic instrumentation, 
a series of procedures for making the measurements and a method 
of reducing the resulting data. Of these parameters, the surface 
electrostatic potentials and the currents are perhaps the 
most important. These are also the most easily measured and 
so should provide a sensitive test of the computer codes. 
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MEASURED 

PARAMETERS 


3-5 Diagnostic Instrumentation 


An array of instrumentation is needed in order to perform the 
measurements defined in the preceding section. These 
instruments should be interfaced to a data acquisition system 
(DAS) which is spelled out in more detail in Section 3.7. 

The instruments which are needed and the measurement which 
they perform are indicated in Table 3.2. 


Diagnostic Instrumen ta ti on 

T r 



Field Mill 

Electrostatic 
Vol tmeter 

E Beam/RPA 

Electro- 

meters 

1 

Magneto- 

meters 

Charge/Trajectories 

X 





E Field 

X 






Potential 


X 

X 



Total Current 


j 


— ■ 1 

X 

pT" 

B Field 

1 

1 



' 

X 


Table 3.2 

Description of the instrumentation follows. 
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3.5-1 Field Mills 

The electric field is best measured with a time varying 
capacitive probe, i.e. a field mill. The basic equation 
describing such a device is 
Q = CV 


whence the current 


I 



C 


dV ® 
dt 


+ V 


dC ® 
dt 


Term 0 is the usual term which one finds in ordinary 
circuit analysis of constant capacitance devices. Term 0 
is important whenever the capacitance varies in time. In 
order to understand how such a device can be used to sense 
electric fields we consider the situation illustrated in 
Figure 3.5. 


For a plate of area A separated by a distance d one has a 

A 

capacitance C = t Q ^ 


And an electric field 
E = V/d 


where V is the potential difference across the 
plates 


Thus the current associated with the time varying capacitor is 


t - * dE , „ dA 

I - £ A ~rr + e E XT 
0 dt o dt 


A (t) varies periodically in time with a fundamental frequency 

to thus A(t) = aZa r e inu>t 

•• n 
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ideal rectifier 


Figure 3.5. Elementary Field Mill 





For purposes of illustration we take 
A(t) = 7 ? a (1+sinait) 
and thus 

I(t) = e o a dt + e 0 E \ awcoswt 

For a static electric field *• o 

a b 

and the voltage at point x in Figure 3.5 is given by 

v«) = k«)z<«) = bo 

For an ideal integrator. 

I Z | — ► for u)RC> > 1 

lot 

thus 



The field mill therefore provides a simple way of measuring 
the electric field at a surface. In actual instruments one 
can use substantially more sophisticated filtering and phase 
sensitive detection techniques to improve the sensitivity and 
accuracy of the field mill. Also in practice one would 
calibrate the field mill using known fields. The analysis 
above is thus meant only to illustrate the principles involved 
and would not be used to determine actual field strengths. 

A simple analysis of field mills has been provided by Seeker. ^ 
Field mills are commercially available from Monroe Electronics, 
Inc., 100 Housel Avenue, Lyndonville, New York, 14098. The 
Monroe units use a rather large sensing element and so it 
may well be necessary to design and fabricate a smaller sensor. 
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The device thus far discussed is useful for measuring the 
electric field strength near a flat surface. It must be 
modified however to measure the local electric field at a 
point in space well away from a surface. Our survey did not 
turn up a commercially available device which would perform 
this measurement,, thus we have again, as in the case of the 
distributed source accelerators, been forced to propose a 
device which will require some development. 

The device used to measure the electric field at a point in 
space well away from a surface is illustrated in Figure 3.6. 

A variation on tne earlier theme is used to measure the electric 
field at a point in space well away from a surface. A 
rotating sensor consisting of two small spheres of radius a 
separated by a distance d is used to drive an integrating 
differential amplifier. The rotating sensor would be 
mounted on the end of an extensible shaft. This shaft could 
be extended from the top of the vacuum system downward 
into the measuring area. Alternatively the sensor would be 
extended outward from the wall. The mechanics of this one 
dimensional movement should not be too difficult. 

As is shown in Figure 3.6 the potential on the spheres varies 
as they are rotated. Thus one senses a term like 0 in the 
current equation. By integrating this term one finds a* and 
hence can directly infer E = A*/d. This device and the 


237 



conventional field mill previously discussed are both simple 
theoretically and easy to construct and calibrate. The 
rotating sensor device will require perhaps six months of 
developmental work. 
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Figure 3.6 Rotating Sensor Field Mill - The potential A(J) impressed 
between two small spheres rotating with frequency is measured 
with an electrometer. The local electric field E = A<j)/d is thus 
determined. 
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3-5-2 Electrostatic Voltmeter (ESVM) 

A field mill can be combined with a phase sensitive detector 
and a feedback amplifier to make a sensitive electrostatic 


ROTOR 



Figure 3- 7- Block Diagram of Electrostatic Voltmeter 
voltmeter. A block diagram of such a system is shown in Figure 3.7. 
The voltage out of the field mill is proportional to the electric 
field E = (V - Vref)/d as has been demonstrated in the previous 
section. The voltage Vref is served through the integrating 
amplifier to match V so as to reduce E to zero. Thus Vref 
approaches V. Typical commercially available systems have an 
accuracy of 0.1# and will settle to within this accuracy in 
10 millisecs . Instruments of this type are manufactured by 
Monroe Electronics, Inc., 100 Housel Avenue, Lyndonville, New 
York and by Trek, Inc., 8460 Ridge Road, Gasport, N.Y. 14067. 
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The Trek unit uses a probe based upon a mechanically resonant 
reed which leads to a smaller sensor than is available from 
Monroe . 

3*5*3 Diagnostic Electron Beam (DEB) and Beam Locator (BL) 

In addition to the primary electron sources which simulate 
the magnetospheric fluxes an electron beam should be available 
for diagnostic purposes. Simple inexpensive guns designed 
for use in CRT’s are available. These guns produce a small 
diameter beam of monoenergetic electrons and work well in the 
100 ev to 20kev range. The CRT guns are provided with plates 
for x-y deflection of the beam and lenses for focusing. The 
Beam Locator consists of a four quadrant electrometer as 
shown in Figure 3*8* The beam is located when signals from 
the four quadrants are all equal. The electronic Beam Locator 
is not commercially available but should be relatively simple 
to develop. We estimate that three man months of engineering 
effort should be sufficient to produce a working unit. An 
alternative approach which is of particular use in debugging, 
is to use a photoemissive screen. With such a screen one 
simply observes a glowing spot where the beam is impacting 
on the screen. 



INDICATING 

LAMPS 


LAMP SENSING 

DRIVERS ELECTROMETERS 

Figure 3.8. Block Diagram of Beam Locator 
The sensing electrometers produce signals indicating which quadrant 
the beam is in and drive the lamps proportionately. Provision 
can easily be made for handling a large variation in beam current. 
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3-5.4 Retarding Potential Analyzer (RPA) 

Retarding potential analyzers have many uses in this program. 
A unit based upon development work done by MAYA personnel is 
shewn in Figure 3-9- The RPA measures the flux of particles 
with energies E>V^ which enter the entrance apperture of the 
sensor. Used in conjunction with a diagnostic electron beam 
and beam locator, the RPA will provide detailed information 
on the particle trajectories and electrostatic fields. 

An RPA must use a current sensor to detect the flux of charge 
particles. The MAYA RPA is based on a spiral electron 
multiplier which counts individual particles, and hence is 
capable of sensing very low fluxes of particles. An RPA can 
also be designed around an electrometer. (See Section 3 . 5 . 5 .) 
The RPA based on the spiral electron multiplier will operate 
at approximately 10 6 times as low a current density as will 
an instrument based on a conventional electrometer. 



INSULATING SPACER 

n APERTURE PLATE . MOUNT !NG 



i- 

ront of a spiral 
at have energies 
building most 
cts. 

CD 

4- -C >> *3 

rsi 

r— L- 

>> 

0) *M +J 


-C c C 

03 

+-> c <d o 

C 

O OO o 

<c 

O *M CD 

r* 

4-> O S~ CO 
L. Ql 3 

03 

*o a o 

*r- 

<U CO *r— 


*-> S- 

c 

O ro 

CD 

< > 


4- LO >- 

O 

c <£ i_ 

Q_ 

co o a> 


•«“ S- "C 

4-j c 

c 

>i U *3 

*r* 

*— <U r— 

■o 

-Q i— co L- 

L. 

E 0) 

03 

a> 4-> ni 


CO CO C >1 

03 

CO CD CD r— 

Qd 

03 tO 4-> <TJ 

a> 

c o c 
<d a. 03 

03 on 

* 

•*“ Die 

ro 

-m x: c: 03 

<u 

c U *r- 

CD *»“ X3 x: 


-M -C S- U 

3 

O i 03 3 

cn 

CL 4-> CO 

•r— 

L- CD 

Lu 

CJ) <D 4- 


c •*“ o 

— a; 

T3 Q-X: CO 
S- -r- 4-> -M 
03 -M C 

4-> r— C CD 
CD 3 03 E 
w E jz <D 

-M r— 

03 C 01 

•— o u 

Ql i- CD CD 
E+J4J c 
•r U (OD 
CO CD CD 
«— 4- 

< CD CO O 


244 



1 - 5-5 electrometer (EM) 

The direct measurement of current fluxes and of total leakage 
currents to surfaces on the test object is probably best 
accomplished with standard electrometers. An electrometer is 
simply a high input impedance amplifier which measures very 
low currents (less than 10 9 amps) by the voltage which they 
produce in very large resistors (greater than 10 9 ohms). A 
number of electrometers suitable for specially designed 
electrometers, such as might be used in the test module, 
are also readily available. 

3-5*6 Magnetometer (MM) 

A survey of the magnetic field within the vacuum tank should 
be performed. An accuracy of 1^ or 3 milligauss should be 
quite adequate. A number of commercially available magnetometers 
are capable of this accuracy. 
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3.6 Procedures 


The diagnostic instrumentation outlined in Section 3-5 is 
used in standard procedures to make those measurements 

called for in Section 3 A. These procedures are outlined in 
this section. 

Procedure 1 - Trajectory measurements 

This procedure is used to investigate the charge distribution 
and to verify the predictions of subroutine PUSHER. This 
procedure makes use of the Diagnostic Electron Beam (DEB) 
and Beam Locator (BL). A Retarding Potential Analyzer (RPA) 
may also be used for detailed measurements of the electro- 
static potential. 

The procedure is: 

1) Locate the DEB on an outer surface element. By an 
outer surface element we mean a place on the 
outer boundary of the region of interest, say at 
the faraday cage surrounding the test region. The 
nomenclature here is chosen to match that reported 
by S^ in the interim tech report SSS-R-3i24 lf . 

2) Define the DEB energy and direction thus determining 

v (9,$, v g ), the velocity vector characterizing the 
electrons from the DEB. 

3) Place BL directly in front of DEB. 

Trace trajectory outward from DEB with the BL. 
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The last two steps are conceptually straight forward, but may 
in practice, prove to be quite difficult. The Beam Locator 
must move along the three dimensional locus of a particle 
trajectory if the particle trajectory is to be traced. The 
mechanism required to make this movement within a vacuum system 
does not appear to be simple. What is required here is 
remote manipulator technology. Solution of this problem will 
require further developmental work. 

To determine potentials using the RPA: 

1) Perform the previous procedure. 

2) Replace the BL with the RPA. 

3) Perform energy analysis on DEB with RPA. 

*+) Process data to extract RPA charging effects. (This 
might be difficult and will require more study.) 

Again the mechanical movement may well be a major problem, 
otherwise the approach is conceptually simple. 


247 



Procedure 2. electric Field Measurements 


This procedure is used to investigate the electric field 

distribution and hence to determine the charge and potential 
distributions . 

2.1 Surface Fields. The field mill discussed in Section 3-5.1. 
can be used to measure surface electric fields and hence to 
infer the surface charge. The procedure consists simply of 
placing the field sensor near the surface and making the 
measurement. To avoid shadowing of the bombarding particle 
beams it will be necessary to move the field sensor across 
the surface more rapidly than the charging time. 

2.2 Electric Field Distribution. The direct measurement of 
the electric field distribution at points well away from 
surfaces requires an extension of the field mill techniques 
to make them applicable to free space. The rotating capacitive 
sensor described in Section 3-5-1 and illustrated in Figure 3.6 
measures local electric fields. The rotating sensor probe 
would move on an extensible shaft through the measurement 
region of interest. This motion is essentially one -dimensional 
and thus should be far less difficult to implement than the 
troublesome three-dimensional movement required for the Beam 
Locator. The procedure would consist of simply extending the 
sensor to the appropriate point and reading its output. 
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Procedure 3 . Direct Measurement of Potential 


This procedure systematizes the use of electrometer probes 
to measure the electrostatic potential. 

3.1 Surface Potential Measurement. The procedure is 
essentially the same as for the measurement of surface 
electric fields, except the feedback amplifier system described 
in Section 3*5*2 is employed with the sensing element. Again 
care must be taken to avoid shadowing effects. Here the small 
probe provided by the Trek, Inc., Model 340, electrostatic 
voltmeter will help. 



Procedure k. Measurement of Total Current 

The measurement of total current to the test body is 
accomplished by an electrometer which is connected to the 
test body through a variable resistor R. The procedure is 
analogous to that which would be applied to measure the 
characteristics of a vacuum diode with the particle source 
playing the role of the cathode and the test object that of 
the anode. (See Figure 3.10.) 

The procedures described below will not work if the test 
object is a perfect insulator. If however the test object is 
a conductor or if leakage currents (10" 10 amps) are present 
then the instrumentation should be relatively easy to 
implement. The resistor R should be attached to a conducting 
part of the test object. Insulators and other parts of the 
test object will each charge to equilibrium potentials which 
may be different from the potential at the point where the 
resistor R is attached, but which should be predicted by the 
S 3 computer codes. Similarly the total current collected 
through the resistor R should be predicted by the computer 
code. Care should be exercised in the choice of test objects 
so as to avoid absurd demands on the instrumentation. For 
instance, a conducting sphere completely coated with a very 
good insulator would yield too low a leakage current for useful 
measurements to be performed. Nonetheless the computer codes 
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should be aDle to predict, and the experimentalist reasonably 
measure the resulting leakage current through R for simple 
test objects consisting of mixtures of conductors and 
insulators. We note in passing that this technique is the 
generalization of existing work at Lewis as described in 
Berkopec, et. al.^ and in Purvis, et. al.^ 

The procedures are: 

A) Equilibrium I-V curve. 

1) Define energy V Q of particle source (s). 

2) Vary R from Rmin to Rmax and measure total 
current with an electrometer. 

3) Plot I-V curve. V here is the voltage between 
the point of attachment of the resistor R to 
the test object and the outer faraday cage. 

4) Repeat for new V Q . 

B) Dynamic Response for small signals. 

Several approaches are possible. 

1 ) Define energy by V Q = V Q + V 1 sin u t 

2) Measure oscilliatory component of current as 
function of u>. 

3) Plot normalized frequency response vs. w in order 
to estimate lumped parameter values. 
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PARTICLE OVERALL IMPEDANCE OF MEASURES 

SOURCE TEST CONFIGURATION LOAD CURRENT 



(c) 


0 Olfist set-up for measuring and analyzing total currents a) A flux 
of charged particles is incident upon the test object b)I-V ?urve fo? 
test object in beam is analogous to that of a vacuum diode c) Equivalen 

ba ? ed upon e1th ® r Jumped parameter model (Inouye) or three 
dimensional computer code (S 3 ). 
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Alternatively J or R may be varied in an oscilliatory 
fashion to obtain frequency responses. From these responses 
one can determine the magnitude of the various circuit 
elements in a lumped parameter model of the system. Note 
also that the small signal dynamic response techniques can 
be extended to the study of small variations in parameters 
defining trajectories by employing the diagnostic electron 
beam and the retarding potential analyzer. For example, 
consider an instance where the test object is a metallic 
sphere charged by an electron beam to a negative 10 Kilovolts. 
Electrons from a diagnostic beam with energies greater than 
10 Kev would penetrate to the test object. The electrons 
with energies less than 10 Kev would be reflected before 
reaching the test object. Detailed trajectories would be 
predicted by the computer code. If the Diagnostic 
Electron Beam current were modulated at frequency ui and 
then the DEB energy was slowly swept through 10 Kev a 
signal at frequency w would suddenly appear on the oscilloscope 
when the electron energy was sufficient to enable them to 
reach the test object. An RPA could be positioned to detect 
the reflected DEB and thus provide another independent 
measurement of the test object potential and of the electron 
trajectories. Variations on these themes are endless and 
only the experimenters ingenuity need limit the possibilities. 
The equilibrium and dynamic I-V curves should be a sensitive 
test of the capabilities of the computer codes. 
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3*7 Data Processing and Display System 

These experiments generate large volumes of data. If the 
data are to be interpreted meaningfully then the data must be 
processed intelligently and displayed in an easily interpretable 
form. The raw data from instruments should also be preserved 
for the purposes of checking any sophisticated analysis 
routines which might be built into the Data Processing and 
Display System (DPDS). Where possible the DPDS should 
present data from the experiments in a form which is directly 
comparable with the data from the computer codes. For 
example, a goal of the DPDS design would be to present 
trajectory plots derived from the experimental data in a 
format similar or identical to that of the spacecraft 
charging computer models . 

The instruments should also have the ability to stand alone 
and be used meaningfully without the DPDS. A system which 
fulfills all of these requirements and which would be cost 
effective could be based upon the recently established 
general purpose instrumentation bus standardized by the IEEE. 

This standard IEEE-488^ has been adopted by a number of 
the major U.S.A. instrument manufacturers, who are beginning 
to provide a wide variety of instruments and controllers 
compatible with IEEE-488. A block diagram of a DPDS based on 
the IEEE-488 standard is indicated in Figure 3.11. 
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Should more versatility and higher data acquisition rates be 

desired for other experimental purposes then a DPDS system 

based upon the Computer Automated Measurement and Control 

or CAMAC standard (IEEE-583) 8 " 10 . As described in the CAMAC 

Tutorial Issue of the IEEE Transactions on Nuclear Science: 10 

"The CAMAC system is a standardized assemblage of modular 

units and the dataway which interconnects them. Dimensions 

are specified for the "crates" (the containers for the 

modules) and the plug-in modules which supply the various 

logic functions contained within the system. Also detailed 

are the interconnection arrangement, including the sockets 

and the interconnecting "highway". These standards thus 

permit mechanical and electrical compatibility between 

equipments supplied by different sources." CAMAC has 

continued to grow in popularity and has gained wide international 

acceptance as a modular instrumentation system for the transmission 

of digital data between instruments and between instruments 

and computers, and computer peripherals. 

The CAMAC system has much to recommend it. CAMAC has been used 
by all of the ERDA national laboratories for almost a decade. 

As a result of this use instrumentation for CAMAC is well 
developed and widely available. Another point of major 
importance for LeRC applications is the increasing interest 
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in CAMAC as a standard data system for use with Space 

Shuttle payloads. Trainer, et al . 11 discusses the use of 

CAMAC and NIM (Nuclear Instrumentation Modules which are 

CAMAC compatible systems developed under ERDA auspices) in the 

space program. Standardized equipment for space shuttle 

experiment payloads based upon CAMAC are discussed by More 
?? 

and Ebert . NASA has let contracts to develop CAMAC 
hardware suitable for space flight thus there is reason to 
believe that the future will see the use of CAMAC in space. 

An environmental test facility which was compatible with 
the CAMAC standard would interface easily with experiments 
from any other group using the standard. All of these 
advantages suggest that CAMAC must be seriously considered for 
while initial costs might be higher than with IEEE-4-88 the 
operational costs could be substantially lower. 

Whichever choice is made MAYA does most strongly recommend that 

either IEEE-V88 or CAMAC be used. Developing a special purpose 
interface will almost certainly drive up costs, increase 
system design and fabrication efforts, result in reduced 
capability, and be incompatible with most other systems. 



APPENDIX 1 


STRUT DESIGN 


The strut which supports the satellite model or other test 
object in the ground test environmental chamber should be 
designed to minimize electrostatic perturbations. An insulating 
rod would be a good choice were it not for the beams or charged 
particles which are present. These charged particle fluxes 
would charge up the insulator and thus create potential 
distortions which would be difficult to predict. So an 
insulating strut is ruled out. 

A conducting strut is called for. If the strut were a single - 
conducting rod, then it would create a potential surface which 
would not match closely to either the inner or the outer 
potential surface. We overcome this problem by making a strut 
of a series of conducting sleeves, each of which is run at a 
difl* e rent potential. While this technique is somewhat more 
trouble, it will reduce substantially the inevitable perturbations 
of the electric potential by the strut. 


As an example, and as a probable good approximation for use 
in more complex geometries, consider a strut with N segments 
which is designed for use in a cylindrical test system. The 
potential is held at ground on the inner cylinder and at $= V 
on the outer cylinder. The resulting potential is: 
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The radius at which the potential is equal to is 



4 > 

where a = ~ 

v 

Say *£ = n V where both integers 

then A J At 

r. = «'l» J b W 


The segments are chosen to extend between r and r 

£ £+j 

the potential chosen as an average of $>! and $ 
example with 10 segments is shown in the accompanying 


, with 
An 

figure. 


More generally when the general shape of the potential is 
known beforehand one chooses the strut segments to match the 
equipotential contours. In practice, for most purposes, the 
logarithmic strut indicated above should be quite adequate. 


* 
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APPENDIX 2 - DISTRIBUTED SOURCE ION ACCELERATOR 

An ion accelerator which produced a bean of particles with a 
finite and variable spread in energy would simulate the 
nagnetospheric particle flux much more accurately than does 
a monoenergetic beam. In this appendix we outline a possible 
approach to the problem of designing such an accelerator. 

The accompanying figure shows the approach. Ions are produced 
by electron impact on a crossed neutral beam. The ions are 
then accelerated (at right angles to the crossed' electron 
and neutral beams) down a drift tube. Each modular stage in 
such an accelerator produces ions and accelerates ions from 
all of the preceding stages. The final output of the device 
would be a beam made up of a number of ion beams of differing 
energy. It should be possible by varying the electron and 
neutral beam currents to vary the number of particles at each 

energy. A final stage would post accelerate and diverge the 
entire beam- 

Now let us analyze a single stage of the module. First 
consider the ionization process: 

an. 

It = r e« n n 

where n^ = number density of ions 

r g = number flux of electrons 
Q = ionization cross-section 
n n = neutral beam number density 
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The ionization region, within which the crossed beams collide, 
is l b wide on the electron side, l fl wide on the neutral side 
and d thick. So the volume is V = l^d and the rate of 
production of ions within the stage is 

~ = ~ V = Vn^^ 

We consider the emission limited current from the ionization 

region for which: 

3N i 


is the ion current 


X i = q i "It 
where 

is the ion charge 

and hence the ion current density is: 


J i - f 


A = lwl„ 
d n 


q i 

= ~ (Qnd)J e 

where J Q = e r e is the electron current density. 

A free jet can produce beams of neutrals with pressures p - lO’-torr 
corresponding to number densities n-IO^/cm^ (Anderson, 197 ^ 3 ). 

From McDaniel 1 ^, we note that Q-IO'iW for lOOev electrons 
bombarding molecular hydrogen. 

An electron number flux of 10>«/W-sec corresponding to a 
current density of 1.6 ma/cm 2 is easily obtained. 1 ^ 
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Thus Qnd'10 5 for d - .1cm, which should be easy to construct. 

So that for singly ionized particles one has 

J, = 10" 5 J 
l e 

The desired ion current densities (See Table 3.1) should not 
be difficult to obtain. 

Finally we note that the beam intensity might well be subject 
to drift induced by variations in the electron or neutral 
beam currents. An electrostatic retarding potential analyzer 
should be used as a sensor in a feedback control loop to 
stabilize the beam intensity at each energy. Such a stabilization 
loop should not be difficult to design and construct. 

No Distributed Source Ion Accelerator of the type discussed 
in this appendix has ever been built. The physical principles 
are straightforward and well understood but the detailed 
engineering of such a device would be a substantial undertaking. 

At least two man years of scientific engineering effort would 
probably be required in order to build a working unit. 
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APPENDIX 3 - DISTRIBUTED SOURCE ELECTRON ACCELERATOR 


An electron accelerator which produces a beam of particles 
with a distribution of energies is needed to simulate the 
magnetospheric electron flux. Such a gun could be constructed 
by reversing the extraction potential applied to the 
distributed source ion accelerator. We sketch here another 
alternative, based on the use of photoemission as a source 
of electrons. 


U-V lamps are used to illuminate the back side of a wire grid. 
(See accompanying figure.) Photoelectrons produced in each 
stage are accelerated through each of the subsequent stages. 

A composite beam is thus made up of a number of beams of 
different energy. A final stage of post acceleration shifts 
the entire distribution to some mean energy. 


Consider a single stage. The processes involved are: 


Stage 1 


^ - electron flux from i u stage 


Photon Flux 


0 

1 


Grid Wire 


th 


Jj - electron flux from j stage 


One has : 

where 

<ij is the effective photoyield of stage j 

is the effective secondary emission 
** yield at the jth stage from electrons 
which are emitted at the ith stage 


265 




Figure A-3 Distributed Source Electron Accelerator 




The second term in this expression represents the effects of 
electron multiplication. These effects should be kept small. 
The beam current will then be controlled primarily by the 
photon fluxes which are in turn easily controlled. 


Now let us analyze the effect of electron multiplication. The 
worst case results from the final stage, N. The electron 
multiplication current is 

n ' j <n n3 J j 


* 6 


(n-1) 


Require that 


n*« « <n - 1) 


Note that 

8 


= (u) 


where 


aA 


is the ratio of the grid area to the 
total area 


6 is the secondary emission yield 


thus 


’t- [ft )•]“"’ 


for gold & < 6 TnQV . = 1.45 and ^ ^0.1 should be easily 


max 

realizable so that 


j£ = (0.l45) n_1 = 2.8x10’ 8 for a system with 
i 10 stages 

Electron multiplication is not a serious problem 
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We note that the Distributed Source Electron Accelerator 
discussed here has never been built. MAYA personnel have 
built mono-energetic electron sources using photoemission. 
These sources were used for calibration of the auroral 
particles experiments on the Advanced Technology Satellites 
5 and 6. The photoemissive sources were easy to use and 
worked well. Should this concept be extended to build the 
Distributed Source Electron Accelerator discussed in this 
Appendix then approximately one man year of engineering 
science effort should be allocated to design, construct 
and test a prototype unit. 
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APPEND: 4 - COMPARISON BETWEEN THE GRQjD TEST AND THE MAGNETOSPHERIC 


ENVIRONMENT 

The ground test environment is designed to simulate the 
magnetospheric environment and to test the predictive abilities 
of the computer codes which describe spacecraft charging. 

The first task, simulation of the magnetospheric environment, 
is a formidable problem, and admits to only a partial solution. 
The testing of the computer codes is a substantially more 
tractable problem. 

The difficulty with simulation of the magnetospheric 
environment is one of scaling. The plasma surrounding a 
satellite at geosynchronous orbit has four important components: 
the background magnetospheric electrons (me), background 
magnetospheric ions (mi), electrons produced by secondary 
emission from the satellite surface (se), and the photoelectrons 
(pe) which result from solar illumination of the satellite 
surface. The Debye lengths associated with each of these 
components typically stand in the ratio 

L d (mi): L d (me): L^tse): L^tpe) = 1000:^00:2:1 
where typical values are as indicated in Table 1. Note 
that the photoelectric Debye length is by far the shortest 
(^(pe)^ meters) and is the only one which is comparable to 
the dimensions of most present satellites. Thus for most 
purposes one can assume that except for the photoelectrons 


the Debye length is much larger than the typical dimensions of 
the satellite. Thus it is not practical to attempt to simulate 
the Debye shielding of any of the particles except for the 
photoelectrons. Nor fortunately is it particularly important 
to do so. 


Lab 

(B*3 gauss) 

1 ) Background Electrons 

n^/cm 3 ’L^(me)'v1x10 4 cm 

T^5 Kev R(me) / v560cm 

r^l . 8x1 0 1 °/cm 2 -sec 

2) Background Ions 

n^/cm 3 L ( j(mi)'v5* 3x10 4 cm 

T ^30 Kev R(miH5- 8x10 4 cm 

r^1x10 2 /cm 2 -sec 

3) Secondary Emission 
Electrons 

n , ''125/cm 3 L^(se)^1l5 cm 

T^3 ev RCseHlif cm 

r^O^xlO 1 0 /cm 2 -sec 
*0 Photoelectrons 

n-vi+50/cm 3 L d (pe)^50 cm 

T-v-2 ev R(pe)*v11 cm 

r^2 . 6x1 0 1 °/cm 2 -sec 


Magnetosphere 
(B^ICf gauss 

L ( j(me)^1x10 4 cm 
R(me)^1 . 7x10 5 cm 


( mi ) ^ 5 • 3x1 0 4 cm 
R(mi) , v1 . 7x1 0 7 cm 


L^Cse)^ 1 5 cm 
R(se)^. 1x10 3 cm 


I' ( j(pe)' v '50 cm 
R(pe)y3Ax10 3 cm 


Table A-1 Characteristic Plasma Parameters and Scale Lengths 
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More severe problems arise when one considers the scaling 
with cyclotron radius, R. For the magnet ospheric plasma one 
has L d <<R thus for all practical purposes the sheath may be 
considered unmagnetized. In the lab, if no provision is made 
for bucking out magnetic field of the earth, then one has 
L d >>R. Indeed for the photoelectrons, the cyclotron radius 
will be comparable to the size of the test object. This 
condition only occurs in the magnetosphere when the satellite 
becomes as large as one hundred meters, and even then, the 
Debye sheath should dominate. Thus if the background magnetic 
field of the earth is not reduced within the test chamger, 
simulation of the magnetospheric environment will be problemati 

We note though that the computer codes are being developed 
to include a magnetic field. Tests of the predictive ability 
of the codes do not require bucking out the magnetic field. 

The codes should be capable of computing the magnetized sheath 


structure. 
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APPENDIX C 


A VARIATIONAL/FINITE ELEMENT APPROACH FOR SATELLITE 
PROBLEMS OF ARBITRARY GEOMETRY 

The following will describe a variational formulation 
for a certain class of charged satellite problems. Using a 
finite element approach to the solution of the variational 
problem, a symmetric system of linear equations results. 

Consider a conductor with total surface charge Q and 
« local surface charge density a and at constant potential <t> , 

C 

comprising surface C^, surrounded by volume V with fixed 
charge density p and dielectric constant e, which is in turn 
bounded by surface C2 on which the outward normal component 
*> E r of the electric field is specified (see figure) . The 

only requirement on p, e and E n is that they be integrable; 
thus, p may have an (integrable) singularity representing a 
fixed surface charged distribution. After the equations are 

• formulated, "geometric" boundary conditions — i.e., fixed 
potentials — may be applied on surface C2 by a simple modi- 
fication of the equations. 

Note that a dielectric layer on the conductor is con- 
sxdered here as part of volume V in which the potential <f> 
is unknown and to be solved for. Obviously, the volume 
integrals which must be evaluated in any finite element 
scheme will reduce to surface integrals in the elements com- 

* 

prising dielectric layer. 

If one assumes a variational principle for the system, 
of the form 

= 0 (C.l) 


Preceding page blank 
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where 



f f(x,0,V0) 

V 


dV 



g(x,<j>,a,i|> c ) dS 


h (x, 0) dS 


(C.2) 


then one obtains 


■*•/••!» -fer 


Bf 


V 


3x^ 3 0 , ^ 


dV 


/ { 6 *[lf + fi i Ht- 


+ so If + 6* c |f^ j as 


/ S*(fr + n, 14— i as - 0 


3<f “i 3*^ 


(C.3) 


where n is the unit normal vector (directed out of V) and 


*'± 5 Hr 


(C. 4 ) 


Now 0 and a may be varied independently in their 

0 respective domains; however, since 0 is constant, the 

c 

variation 60 c is constant over C^. Thus, the variational 
equations are 
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U A. u o A- n 

in V 

(C.5a) 

3g , 3f 

tt + n . « =0 

o $ 1 

on 

(C. 5b) 

If ■ ° 

on 

(C. 5c) 

/ Ur ds - 0 

n ^ 


(C.5d) 


ah - 3f 

39 n i 9$7~ * 0 ° n C 2 (C.5e) 


Choosing 


f » | V<J>.V<f> - p<J> 

g = (2 . a ) (*-4, c > - 2 * 


(C. 6a) 
(C.6b) 


h = E <p 
n 


(C. 6c) 


where A is the area of surface the variational equations 
(C.5a) - (C.5e) become 


V* (-eV<{>) 

= P 

in V 

(C.7a) 

-n • (-eV0) 

* a 

on 

(C.7b) 

* ■ *C 


on 

(C. 7c) 

J a dS 

= Q 


(C.7d) 

C 1 

n • (— e ) 

= E 

n 

on C 2 

(C.7e) 
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which are the desired governing equations of the system 
(in rationalized MKS units) . 


Suppose the volume V is divided into elements with 
nodes i and interpolation functions N 1 at each node. The 
"surface" elements on C 1 and C 2 will be faces of elements 
of V. Expanding <J> and a in the N 1 , i.e.. 


<t> = <p x N 1 


a = a 1 N 1 


there results 


<£ = J § <P Z <p j v^-vn^ - p^V dv 

v' ' 

+ O 1 N 1 ~ (P 1 a j N 1 N j - | <J> c | dS 


+ j E n <t > ** N 1 


dS 


(C. 8a) 
(C.8b) 


(C.9) 


Defining 


: u ; / £ 71,1 


•VN 3 dV 


(C.lOa) 


Pi - f PN 1 dV 
V 

s ±j 5 / 1,1 Nj 


dS 


(C. 10b) 


(C.lOc) 
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(C.lOd) 


= f N 1 dS 
C 1 

E i- f E n 1,1 ds 


one has 


(C.lOe) 


<£= j ^ - P i (p 1 + 


+ E. <p ' 

i 


Extremizing £ with respect 
suiting equations are 


£ ij * 3 - s ij ° 3 - p i * E i 
s i *c - S ij * 3 - 0 


a =» Q 


s i ol * C “ S ij * X ’ Q ♦c 

(C.ll) 

to <p* f a 1 and <p , the re- 

C 


(C.12a) 

(C.12b) 

(C.12c) 


Of course, matrix elements S i; . can be nonzero only if 
nodes i and j are both on C^; s^ will be nonzero only if 
node i is on C^; and E^ can be nonzero only if node i is on 
^ 2 * Furthermore, c is defined only on C. , so the degrees 
°f freedom a , where node i is not on C^, are to be removed 
from the system. Equation (C.12a) is the discretized analogue 
of the expression for the charge on C^. 

From the definitions of Equation (C.lOa) and (C.lOe), 
and the matrices e.. and S. . are symmetric. 
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The finite element equations (C.12a-c) can be re- 
formulated in matrix-vector form; the writer begs the 
readers’ indulgence with "civil engineering" notation. 

The equations may be expressed as 

[K] {q} « {R} (C.13) 

The matrix [K] is symmetric, but has some zero diagonal 
elements and, in general, is indefinite; the solution is 
not a minimum of , but a saddle point. However, this 
circumstance is easy to remedy, while at the same time re- 
ducing the number of degrees of freedom to be solved for. 

Assume the interpolation functions obey the usual rule 

En 1 - 1 (C.14) 

i 

From Equations (C.lOc-d) 

S. . = s. (C.15) 

ID i 

Now the variational expression <£ in Equation (C.ll) may be 
written 

+ {s) T {a} <J> - {<j>} T [S]{a> - Q <f> + {E) T {<|>} (C.16) 

^ W 

where the dimensionality of vectors and matrices is N and 
N x N, respectively, where N is the total number of nodes. 

From Equation (C.12b), assuming [S^l is nonsingular, 

{<t»l> * fS n ]“ 1 { Sl } <f> c (C.17) 
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where the vector {<f>} has been decomposed as 



{(|> 1 } is the vector of the i on and {<f> p } the vector 
of the <J> X , i not on C^; the decomposition of [S] and {s> is 
analogous to that of {$}. Equation (C.15) states that 

{s l } " [S 11 3 {1} (C 

where {1} is the vector, with all unit entries, of length 
N^, where is the number of nodes on C^. Thus, 

{4» 1 > - {1} <*> c ( C 

i.e», all 4> , i on C^, are equal to <p^ in the finite element 
solution, as might have been anticipated. Since the solu- 
tion obeys Equation (C.20), we may impose it on £ as a set 
of constraints: 

«C- I {1) T !e 11 ]{1}^ + {1) T U 1P H* F ) * c 

+ i {♦ F ) T [e pF ){* F ) 

- {e 1 } T {l) - fp F > T U p } 

- Q * c + {E p } T {$ p } 

The a and <|>* do not appear in Equation (C.21). 
tional equations for £ , Equation (C.21) are 

e* 4> c + { e p > T C 4> p } » p* + Q (c. 


(C 

The varia- 
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.18) 


.19) 


20 ) 


21 ) 


22a) 


I 



{e F } *C + [e FF ]{ *F } = {p F } " {E F } 


(C. 22b) 


where 


e* = {l> T [e 1 , ] {1} = 2 e.. 

11 i,j«C 1 


{e F } " [e lF ] {1} 


X • 6 * 


£ v 


£pi j f C, ij 


p* = {p 1 ) T {l} = Pj_ 


The system of equations (C.22) is symmetric, and presumably 
positive definite and diagonal dominant, making for ease of 
solution. 

After the system of Equation (C.22) has been solved 

for 0 and the <t>i, the o 1 may be obtained from Equation (C.12a) 

C r 

(S 11 ]{a 1 > = + [e 1F ]{<J> F > - (p 1 > (C. 27 


{ei> * [e 11 Hl> 


i • 6 • $ 


- • k ■« 
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APPENDIX D. POTENTIAL INTERPOLATION IN NASCAP 


D.l LINEARLY BLENDED ELEMENTS 


In the finite element solution of Poisson’s equation, 
the requirement that the potential <j> be continuous is met by 
establishing unique conventions for the interpolation func- 
tions at element boundaries. For our case, we require that <p 
be linear on any edge, bilinear on a square face, and linear 
on a triangular face. Subject to these restrictions, we must 
find the interpolation functions N and weights W, defined 
by 

4>(?) * £nV 

i 

f | V<{> | 2 dr = E w ±j (|> i 4^ 

cell 

for any cell. 

An instructive example is the simple cube 0 < x,y,z < 1 
having an extra edge drawn from (100) to (010). The poten- 
tials on the z = 0 and z = 1 planes are 


$ (x,y,0) 


(1-x-y) <f> Q00 

+ x *i00 + ^010 

x+y<l 

(<x+y-l)* u0 


x+y>l 


<j>(x,y,l) - (1-x) (l-y)<|) 001 +x(l-y)4) 101 +y(l-x)<|) 011 +xy(f> 111 , 


the total potential being 

<J> (x,y, z) = (l-z)4) (x,y,0) + z<|>(x,y,l) 

This last step is an example of "linear blending". 


Preceding page blank 
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Note that the electric field is discontinuous across 
the plane x+y = 1. However, within each of the regions, 
x+y<l and x+y>l, the functions N* and VN 1 are simple. 

In general, the problem is handled as follows: 


1. Divide the cell into regions R q in which a 
differentiable potential function may be con- 
structed by linear blending. This may be done 
in an asymmetric manner, the final result 
being obtained by averaging over equivalent 
functions. 

2. In each R^ write 


VN* (r) 


-? a i 


av 


f (r) 
av ' 


3. Calculate 


(f ay' f av ) 


f 


d r f f 

ay av 


a 


4. Finally, we have 


M[{a}]<j»(r) = E L N* (r) 
{a} i a 


M[ {a}] 


E E E (f 

{a} a yv 


ay 


,f ) a . a . 
av lay jav 


where {a} is the set of equivalent partitions of the cell, 
and M[{a>] is the number of members in the set. We have 
found it convenient to perform steps 1 - 3 by hand, using 
a computer program to perform step 4. 
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D.2 ELEMENT TABLE CODES 

Code for Element Table [LTBL (NX,NY,NZ) ] . 


54321 0987654321 098 765 432109876 5 43210 

t 1 I I I I I I 

E D C B A 


Field Bits 

A 4-0 

B 14-6 

C 18 

D 19 

E 30-21 


Cell-type code (see Appendix D) 

Orientation code (see Appendix D) 

Set if cell is completely filled 
(interior) 

Set for an empty special cell 

Index used to reference PHOJ array 
to determine low energy electron 
currents 


ORIENTATION CODE 


3x3 bits. Each group of 3 contains 1, 2 or 3 in the 
lower 2 bits, with the high bit set for negative. 


m. m_ m. 

Code (-) (-) 2 i 2 , (-) 3 i 3 


takes (r^,r 2 ,r.) to 


r m i m 2 m 3 

L i 2 


e.g. r the following codes take a point to (x,y,z) : 


Octal Code 

123 

365 

532 

176 

567 

617 


Dec . Code 

1 , 2,3 
3, -2,-1 
-1,3,2 
1,-3, -2 
-1,-2, -3 
-2, 1,-3 


Point 

(x,y,z) 

(-z,-y,x) 

(-x,z,y) 

(x,z,-y) 

(-x,-y,-z) 

(y,-x,-z) 
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D.3 


STANDARD VOLUME CELLS 


(Format) 

Description 
Standard Orientation 


Potential Function = £ • 


Weight Matrix, W. fdfi | V<j> | 2 = £ w. .* . d> . 

^ J ij J-ZTi*: 


Point Index 

1 

2 

3 

4 

5 

6 

7 

8 


Cube Corner 

0 0 0 
10 0 
0 10 
110 
0 0 1 
10 1 
Oil 
111 
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Standard Cell 0 

Empty trilinear cube 

Orientation: Arbitrary 

Potential Function: 

„i 


(1-x) (1-y) (l-z) 
(1-z) (1-y) x 
(l-x)y (l-z) 
(l-z) yx 
z (1-y) (1-x) 
x(l-y) (z 
zy (1-x) 
xyz 


0 

1/3 





0 

- 1/12 

1/3 




- 1/12 

0 

0 

1/3 



0 

- 1/12 

- 1/12 

- 1/12 

1/3 


- 1/12 

0 

- 1/12 

- 1/12 

0 

1/3 

- 1/12 

- 1/12 

0 

- 1/12 

0 

- 1/12 1/3 

- 1/12 

- 1/12 

- 1/12 

0 

- 1/12 

0 0 


289 




Standard Cell 1 


Half-Empty Wedge 

1 < x + y < 2 
0 < z < 1 


Orientation : 


Right angle along 
line 7-8 


Potential function: 
,i 


x 

1 

2 

3 

4 

5 

6 

7 

8 


N 
0 

(1-y) (1-z) 
(1-x) (1-z) 
(x+y-1) (1-z) 
0 

(1-y) z 
(1-x) z 
(x+y-1) z 







S?*® 



v.v.v 

'.V.V 

m 



>v>5 

v!v! 

*X*Xv 



.sv.'.'i 




o 


0 

1/4 


0 

1/24 

1/4 

0 

-1/8 

-1/8 

0 

0 

0 

0 

0 

-1/24 

0 

-1/24 

0 

0 

-1/8 

-1/8 


5/12 



0 

0 


-1/8 

0 

1/4 

-1/8 

0 

1/24 

-1/12 

0 

-1/8 


1/4 

- 1/8 
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Standard Cell 2 

Cube with diagonal line on one face 
Orientation: Line from 2 to 3 



Potential Function: 


i 


N 


i 


1 (1-x-y) (1-z) 0 (1-x-y) 

2 [x0 (1-x-y) + (1-y) 0 (x+y-1) ] (1-z) 

3 [y0 (1-x-y) + (1-x) 0 (x+y-1) ] (1-z) 

4 (x+y-1) (1-z) 0 (x+y-1) 

5 (1-x) (1-y) z 

6 x(l-y)z 

7 (1-x) yz 

8 xyz 


5/12 

- 1/8 1/2 


w 

-1/8 1/12 1/2 






0 -1/8 -1/8 

5/12 





7/360 -37/360 -37/360 

-23/360 

1/3 



* 

-11/180 -1/45 -19/180 

-11/180 

0 

1/3 



-11/180 -19/180 -1/45 

-11/180 

0 

-1/12 

1/3 


-23/360 -37/360 -37/360 

7/360 

-1/12 

0 

0 


1/3 
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Standard Cell 3 


Tetrahedron 

2<x + y+ z< 3 

Orientation: Empty corner at 

point 8 

Potential Function: 

1 N 1 

1 0 

2 0 

3 0 

4 1-z 

5 0 

6 1-y 

7 1-x 

8 • x+y+z-2 




0 

0 

0 

0 

0 

0 

0 

0 


0 

0 0 

o 0 1/6 

0 0 0 0 

0 0 0 0 1/6 

0 0 0 0 0 1/6 

0 0 - 1/6 0 - 1/6 - 1/6 1/2 
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Standard Cell 4 


Truncated Cube 


Orientation: 000 corner (point 1) 

missing 


Potential Function: 


1 
1 

2 

3 

4 

5 

6 

7 

8 


N 


0 


exercise 

for 

reader 



0 


0 5/12 


0 

1/72 

5/12 





* 0 

-11/120 

-37/360 

13/3 6 




0 

1/72 

1/72 

-1/9 

5/12 



0 

- 37/360 

-1/9 

-7/180 

-11/120 

13/36 


0 

-1/9 

-11/12 0 

-7/180 

-37/360 

-7/180 

13/36 

* 0 

-5/3 6 

-5/36 

1/45 

-5/36 

1/4 5 

1/45 


0 


7/20 
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APPENDIX E 


The following paper was presented at the IEEE Annual 
Conference on Nuclear and Space Radiation Effects, Williamsburg, 
Virginia, July 1977. (Note that the material properties used 
in this calculation differ somewhat from those used elsewhere 
in this report.) 
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NASCAP, A THREE-DIMENSIONAL CHARGING ANALYZER 
PROGRAM FOR COMPLEX SPACECRAFT 


I. Katz 

D. E. Parks 
M. J. Mandell 

J. M. Harvey 
S . S . Wang 

Systems, Science and Software 
La Jolla, California 

J . C . Roche 

National Aeronautics and Space Administration 
Lewis Research Center 
Cleveland, Ohio 


ABSTRACT 

A computer code, NASCAP (NASA Charging Analyzer Pro- 
gram) , has been developed by Systems, Science and Software 
under contract to NASA-LeRC to simulate the charging of a 
complex spacecraft in geosynchronous orbit. The capabilities 
of the NASCAP code include a fully three-dimensional solution 
of Poisson's equation about an object having considerable 
geometrical and material complexity, particle tracking, 
shadowing in sunlight, calculation of secondary emission, 
backscatter and photoemission, and graphical output. A 
model calculation shows how the NASCAP 'code may be used to 
improve our understanding of the spacecraft-plasma inter- 
action. 


Preceding page blank 
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1 . 


INTRODUCTION 


Several anomalies observed on geosynchronous spacecraft 
have been attributed to differential charging caused by magneto- 
spheric substorms. These substorms consist of hot plasmas with 
mean energies ranging from several thousands to tens of 
thousands of electron volts. The several physical processes 
taking place during a charging event include backscattering 
of incident particle fluxes, emission of secondary electrons, 
charge redistribution on the spacecraft, as well as modification 
of incident fluxes due to the charging itself. However, if 
spacecraft are to be optimally designed for operation under all 
magnetospheric conditions, it is important to understand how 
the different processes interact. To accomplish this, a 
charging analyzer program, NASCAP, has been developed. This 
code self -consistently simulates the three-dimensional charging 
of complex model spacecraft subject to magnetospheric conditions. 
NASCAP is also capable of treating objects in a ground test 
environment . 

In the next section of this paper we will discuss very 
briefly the physical models employed in the code. Following 
that we will present a calculation of a moderately complex 
asymmetric spacecraft similar in size and materials to the 
body of the experimental satellite, SCATHA, which is being 
constructed to study Satellite Charging At High Altitudes. 
Conclusions drawn from this calculation concerning surface 
material properties on satellites are discussed in the last 
section. 

2. COMPUTATIONAL TECHNIQUE 

Magnetospheric substorms are encountered rather fre- 
quently by satellites in geosynchronous orbit. Under these 
conditions the hot plasma environment may be characterized by 
a temperature of 10 4 eV or greater, and a density of order 
1/cm . Such a plasma may charge a satellite negatively to 

298 


i 
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many kilovolts; also kilovolt potentials may develop between 
different parts of the satellite or across dielectric coatings, 
causing dielectric breakdown and various electronic malfunc- 
tions. 

To simulate the charging phenomenon requires self- 

consistent calculation of the electronic and ionic fluxes to 

each surface element of the satellite, and the material 

response to the incident flux. Since the shortest times 

1/2 -3 

characterizing satellite charging l(V/4irRne) (m/kT) ^ 10 

seconds] are significantly longer than plasma dynamic response 
times [(m/4ime 2 ) 1 ^ 2 * 10 -5 seconds] the plasma responds quasi- 

statically to charge accumulation on the satellite. NASCAP 
follows, therefore, an explicit timestepping procedure, each 
step consisting of a charge accumulation calculation and an 
electrostatic potential calculation. A flow chart of the 
NASCAP code is shown in Figure 1. 

The net charge accumulation by each surface cell of the 
satellite is calculated in the presence of the electrostatic 
and magnetostatic fields about the satellite, and specified en- 
vironmental characteristics. In the ground test case, the 
incident flux is provided by a monoenergetic electron gun of 
specified beam profile. In the space case the incident flux 
of electrons and ions at surfaces is determined using the 
reverse trajectory sampling method. The ambient plasma may 
be isotropic and Maxwellian, or may be represented by any of 
several sets of data from ATS— 5 and —6 prepared for S by 
MAYA Development Corporation. Alternatively, a spherical 
probe approximation may be used. Electron backscatter and 
secondary emission due to electron and ion impact, as well as 
photoemission and shadowing, are taken into account. The 
formulation of the various physical processes take into ac- 
count their energy and angle dependence. Details may be found 
in Reference 1. Optionally, a first-order photosheath 
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Figure 1. Flow chart of NASCAP code. 











calculation may be performed, but in general, because magneto- 
spheric Debye lengths are large compared with spacecraft 
dimensions, space charge is neglected. 

The electrostatic potential about the satellite or in 

the test tank is calculated by NASCAP using a finite element 

formulation of Poisson's equation. Under magnetospheric con- 

2 1/2 

ditions the Debye length X D = (kT/4irne ) is typically 
hundreds of meters, so that space charge can be ignored, ex- 
cept for a positively charged satellite which may develop a 
photoelectron sheath. The computational space consists of an 

arbitrarily large number of nested cubic meshes. The resulting 

4 

set of several times 10 linear equations is solved using the 
Conjugate Gradient technique. The satellite or test object 
is defined within the innermost mesh, and may have surfaces 
normal to any of the twenty-six cubic symmetry directions. 

It consists of one or more conductors which may be covered 
with thin dielectric layers. The conductors may be floating, 
held at fixed potentials, or biased relative to one another. 

3. CALCULATION OF DIFFERENTIAL CHARGING OF A MODEL SATELLITE 

The calculation presented here is of a model satellite 
similar in size, shape and surface material to the experi- 
mental SCATHA satellite presently under construction. The 
purpose of this calculation is to examine how the various 
surfaces react to a fairly intense substorm environment. To 
the best of our knowledge this is the first such calculation 
to self-consistently incorporate three-dimensional geometry, 
complex material electrical boundaries, and realistic surface 
response effects. For computational ease the incident cur- 
rents were assumed to be from an isotropic Maxwellian whose 
mean ion and electron densities and energies were chosen to 
be similar to those actually encountered by ATS-5 on March 14, 
1971, during a magnetospheric substorm. At the end of the 
calculation, the self-consistent particle fluxes incident on 



a material test patch on the top of the spacecraft were cal- 
culated using the reverse orbit technique and the reduced 
ATS-5 data. The fluxes agreed with analytical formulas 
within a factor of two. 

The model satellite is shown in Figure 2 and its 
dimensions are given in Table 1. The satellite model con- 
sists of an aluminum right octagonal parallelepiped which 
forms spacecraft ground. All of the sides of the object are 
covered with a four mil layer of silicon dioxide, with the 
exception of a 35 cm band on which the aluminum is exposed. 

On opposite sides and on the top are kapton covered magnesium 
plates. The kapton is 4 mils thick while the plates are 
capacitatively coupled to spacecraft ground by 200 picofarads. 
The remainder of the top surface and all of the bottom surface 
are bare aluminum. Electrically we have four separate con- 
ducting structures and hundreds of separate dielectric sur- 
faces. Coupling is purely capacitive except for the very 
small dielectric conductivity. In this calculation the di- 
electric conductivity was independent of the field strength 
and had a numerical value of 10~ 14 mho/m. 

The calculation of backscatter and secondary coef- 
ficients was performed as described in Reference 1. 

Initially, the satellite is assumed to be in eclipse 
with all surfaces at zero potential. Table 2 shows a break- 
down of the particle fluxes both incident and emanating from 
the four different types of material surfaces. It is most 
interesting to note that even though there is no photoemis- 
sion, secondary emission from the insulators, in particular 
from SiC^, exceeds the incident electron flux. Notice also 
the high electron secondary emission caused by incident pro- 
tons. As such then, a satellite completely covered with SiO, 

mm 

would actually charge a few volts positively with respect 
to the plasma. It is the aluminum surfaces that are charging 
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Figure 2. Model satellite used in sample calculation showing 
surface resolution and materials. 
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TABLE 1 


Dimensions of Model Satellite and Comparison 
with SCATHA Spacecraft 


SCALE — 1 Zone * 11.5 cm 


Dimensions in Meters 




MODEL 

SCATHA 

Octagon Mean Diameter 

S 

1.68 

1.71 

Height 

s 

1.73 

1.77 

Aluminum Band 

= 

0.35 

0.25 

top 

35 

.35 x .35 


SSPM 



.32 x . 

sides 

= 

.35 x .33 


(For this calculation all SSPM 

are Kapton 

over Mg) 


32 
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TABLE 2 

Charging at Early Time 
(Satellite Potential = -80 V) 
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negatively and driving the overall satellite potential nega- 
tive. 

As one would expect/ the overall body goes negative at 
about the same rate as an equivalent sphere with the same total 
net charge. However, the maximum potential achieved is only 
575 volts, just over 10 percent of the mean plasma energy. The 
potential about the satellite is illustrated in the upper part 
of Figure 3. After developing this voltage, which takes less 
than 300 milliseconds (see Figure 4), the surfaces begin to 
develop differences of potential. The ratio in time scales 
between differential charging and overall charging is approxi- 
mately the same as the ratio of the overall body size to the 
dielectric thickness. In this case that ratio is the order of 
10 «. 

Figure 4 shows the highest differential charging rate 
is across the 4 mil Si0 2 surface. This rate is only about 
75 vo It s/cm- sec , or less than a volt per second. The actual 
potentials grow more quickly on the kapton surfaces in spite 
of a lower differential current. This is due to having less 
capacitance per unit area relative to the underlying metal. 

It is important to note that the Si0 2 surfaces charge positive 
with respect to spacecraft ground but the kapton goes negative. 
This occurs because the net positive current to the Si0 2 tends 
to "soak up" negative current to the aluminum. As a result 
the kapton has an effective higher negative current than the 
aluminum. 

This calculation was run to 56 seconds. The fields in 
the dielectric have been increasing linearly to over 3 x 10^ 
volts/meter with only small changes in Si0 2 fluxes. The flux 
to the top kapton surface had dropped to about two-»thirds its 
initial value with the voltage differential less than 300 volts 
from spacecraft ground. The current balance at this tine is 
indicated in Table 3, and the potentials are in the lower part 
of Figure 3. 
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Early potentials 
No differential charging 



Late time potentials (56 seconds) 

Si 02 covering charged positive 
Aluminum band charged to -575 volts 
SSPM in upper right charged to -650 volts 



Figure 3. Potential contours about the satellite before (upper) 
and after (lower) differential charging. 
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Volts 


^Potential of 
Aluminum Satellite Body 


Potential of SC1-3 


Field Across SiO^ 
Covering (volta/cm) 

d E 

g£ a 75 volts/cm-sec 


Time (sec) 

Time development of potentials on aluminum satellite 
body, on a kapton-coated sample, and the field across 
the SiC^ cover cells. 


It is useful to contrast the results presented here with 

those that would be deduced from a simplified theory based on 

complete current balance on every material of the satellite's 

surface. The simplified theory has been invoked, for example, 
[21 

by DeForest to analyze charging observed on the ATS-5 ex- 
perimental spacecraft. Complete current balance for each sur- 
face material, however, is expected only after charging in the 
given environment has proceeded for a sufficiently long time, 
t g g, and assuming that the plasma environment in the vicinity 
of the satellite does not change for times to £ t . 

In the present application, considering the secondary 
yields of aluminum, kapton and SiC^, one should expect an equi- 
librium in which both SiC >2 and kapton are positive with respect 
to spacecraft ground (aluminum) , and in which SiOj is slightly 
positive with respect to infinity. After approximately one 
minute of charging, the SiC^ is positive relative to the alumi- 
num; the kapton, however, is negative with respect to spacecraft 
ground . 

Based on the net charging currents in Table 3, this in- 
version of the relative potentials of kapton and aluminum is 
expected to persist for several minutes. 

In the absence of dielectric conduction, complete equi- 

—7 2 

libration, based on net charging current j of order 10 amps/m 
and a iV 'v 600 volts is attained in a time of order 


t 


eq 



600 x 200 x 10" 12 x 10 4 «\, 

— —m mxn 

4* x 10~' 


Apart from dielectric relaxation times, t is the longest time 
scale associated with differential charging. In many circum- 
stances, it is longer even than the characteristic time for 
variations in the magnetospheric environment, so that full equi- 
libration of differential potentials may not be achieved within 
the duration of a substorm. Such conditions clearly require a 
quasi-static as opposed to a fully equilibrium treatment of 
spacecraft charging. 
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4. CONCLUSIONS 

The calculation presented here illustrates several 
features that are common to most magneto spheric charging 
situations. First is that all the surface potentials charge 
initially according to the overall net charge on the object, 
regardless of their local currents. For example, the silicon 
dioxide surfaces always have a net positive current, but 
at first their potentials rise along with the rest of the 
body to a negative 575 volts. In general, the potentials 
on all surfaces will initially remain constant with respect 
to each other, but change with respect to infinity in 
response to a sudden change in environment. Then, on a 
longer (minutes) time scale, the surface potentials will re- 
adjust with respect to each other while the integral net cur- 
rent to the satellite is almost zero. 

Second, the insulators, in particular silicon dioxide, 
resist charging due to their high secondary yields and high 
secondary crossover points. This is very important since so 
much of the surface of satellites, particularly those which 
are spin stabilized, is covered with solar cell coverplates 
which are made of SiO^. These materials keep the body 
potentials thousands of volts below a high plasma tem- 
perature. Indeed, it is the exposed conductors that col- 
lect the most charge. This implies that a high degree of 
active control may be possible just by discharging 
the spacecraft ground if the insulating surfaces are chosen 
for high secondary yield coefficients. While this will not 
solve all charging problems, it can certainly lower potential 
differences on the spacecraft. 

The third point is that not only is the magnitude of 
the potential difference between two surfaces dependent upon 
what is happening on the entire vehicle, but even the direc- 
tion of the difference depends on non-local effects. Com- 
paring the current fluxes to the aluminum and kapton surfaces. 



one would assume that kapton would not become as negative as 
the aluminum. While in the limit of long times this may be 
true, the very opposite occurs on a time scale of minutes. 
Since the environment can change in minutes, both situations 
are important. Indeed the system may not even approach the 
long time limit before substantial changes in the environment 
occur. 

The calculation presented here has used just a few of 
the capabilities of the NASCAP code. By performing detailed 
simulations of a variety of both space and ground charging 
situations, we hope to better understand the interaction of 
geometry, environment, and material properties. It is this 
knowledge that should lead to simple design criteria to pre- 
vent charging induced anomalies from interfering with the 
operation of future satellites. 


312 


i 



REFERENCES 


1. Katz, I., et al . , "3 -Dimensional Dynamic Study of 
Electrostatic Charging in Materials," Systems, 
Science and Software Interim Report, SSS-R-77-3124 , 
January 19, 1977. 

2. DeForest, Sherman E. , "Spacecraft Charging at Syn- 
chronous Orbit," JGR, 77, 681, 1972. 



V 




it * 


314 


I 



PRECEDING PAGE BLANK NOT FrT4MED 


APPENDIX F.l 
ELECTRON BACKS CATTER 

Electron backscatter plays an important role in estab- 
lishing current balance to a spacecraft. Quantities of interest 
are: 

1. Backscattering coefficient for normal 
incidence. 

2. Energy distribution (or, at least, mean 
energy) of backscattered electrons. 

3. Angular distribution of backscattered 
electrons. 

4.-6. Dependence of (1-3) on angle of 
incidence. 

r is i 

We have developed a modification of the theory of Everhart 1 1 
and McAfee which gives good agreement with available data 
for 1, 2, 4 and 5 above for 0 < 60° and energies in the range 
10 keV to 100 keV. The theory can be made to yield (3,6) and 
can be modified to yield results in the 1 keV to 10 keV energy 
range . 

The backscattering theory of Everhart, which was extended 

by McAfee to yield an energy distribution, assumes (1) a single 

scattering in accordance with the Rutherford cross-section, 

dE —1 

and (2) the Thomson-Widdington slowing down law, « E 
(valid for most metals for E > 10 keV) . The single-scattering 
approximation is illustrated in Figure F.la. All scatterings 
toward the surface are considered to deplete the beam, though 
only those in region I have sufficient energy to escape. For 
normal incidence, the backscattering coefficient is 

n _ a-l + (1/2) a 
n a + 1 

where a is (somewhat arbitrarily) taken to be 0.045 Z. 


Preceding page Wank 
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For non-normal incidence, it becomes painfully apparent 
that scattering into region II is best ignored; these electrons 
are more likely to escape than those remaining in the beam. 

We adopt the model shown in Figure F.lb: the beam is depleted 

by escaping electrons and by electrons scattered inward into 
a cone of angle 0 c < 9. For normal incidence, this gives 


n = l - 



a 


where we take a = 0.037 Z to give backscatter coefficients 
whose agreement with experiment is similar to that of Everhart. 
The ELTRAN Monte Carlo results for backscattering coefficients 
appear to be somewhat lower than the experimental consensus. 

Darlington and Cosslett^®^ (to be referred to as DC) 
indicate that the angular dependence of the backscattering 
coefficient is well fit by 


&n[n (8)/n (0) 1 = 

1 - cos 6 


(F.l) 


where the right hand side depends only on the material and 
may be assumed independent of the incident energy if the 
Thomson-Widdington law holds. Figure F.la shows this quantity 
as a function of angle for A1(Z = 13) and for kapton (Z = 5). 
For the theory curves, we have taken 


0 

c 


0 0 < 30° 

0-30° 30° < 0 < 60° 

30° 60° < 0 


It is apparent that Equation (F.l) provides an adequate repre- 
sentation of the results, although the ELTRAN data show a syste- 
matic decrease of F with angle. Values of F are shown in 
Table F.l. It is apparent that F decreases with Z, and that 
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Table F.l. 


Angular Dependence of Backscattering Coefficient 


Substance 

z 

E (keV) 

— 

4 

10 - 100 

Be 

4 

25 

Kapton 

5 

10 - 100 

— 

5 

10 - 100 

C 

6 

1-3 

Teflon 

8 

10 - 100 

— 

13 

10 - 100 

A1 

13 

10 - 100 

Cu 

29 

25 

Cu 

29 

1.5 

— 

29 

10 - 100 


p <a) 

Method ^ 

2.54 

Theory 

3.4 

Experiment 

3.75 

ELTRAN 

2.5 

Theory 

1.8 - 2.3 

Experiment 

3.3 

ELTRAN 

2.12 

Theory 

2.24 

ELTRAN 

1.1 

Experiment 

0.75 

Experiment 

1.5 

Theory 


(a) For theory and ELTRAN, F(Z,£) is measured at cos9 » 0.75. 
For experiment, a line is "eyeballed" through the graph of 
fn[n (8)/n (0) ] and its slope is quoted. 


"Theory" is the modified Everhart theory described herein. 
"Experiment" data is taken from Darlington and Cosslett. I 19] 



for small Z the theory gives values for F which are too low, 
while the ELTRAN values are somewhat high. 

Results for the mean energy of backscattered electrons 
are shown in Figure F.2. It is apparent that the theory does 
poorly near normal incidence, where it gives too small a 
dependence of mean energy on angle. 

For E < 10 keV, the Thomson -Widdington law does not 
hold and experiment (DC) indicates increased backscattering 
for low energies, at least for low-Z materials. DC claim 
that the Everhart theory can be modified to give reasonable 
agreement for backscattering at normal incidence in the 1-10 
keV range . 

It is apparent that the theoretical framework can be 
made to yield the angular distribution of backscattered elec- 
trons. Such a calculation would, however, require a non- 
trivial effort, as this function is now integrated out at an 
early stage and the quality of such results would be question- 
able due to the single-scattering approximation. 




0 0.5 1.0 


cos 6 


Figure F.2. Results for angular dependence of backscattering 
coefficients and mean backscattered energy com- 
pared with ELTRAN data. 
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APPENDIX F . 2 

DEPOSITION OF ELECTRONS IN MATERIALS 


We have developed a simple theory for deposition pro- 
files and albedo (net reflection) for electrons characterized 
by a plasma temperature of 1-100 keV. We find that the 
deposition profile can be adequately represented by a simple 
exponential. 

The underlying assumptions of our theory are: 

1. The electron plasma is isotropic and character- 
ized by a temperature, T. 

2. An incident electron is either backscattered 
or travels in a straight line to its "maximum 
practical range." 

3. The backscattering probability is taken to be 
energy- independent and given by (Appendix F.l) 

t) = n Q exp [n^ (l - cose)] 


n 


o 


l 



0.037Z 


n 1 = -log n Q 


It follows that the albedo A q is given by 

a q = (2n 0 /n 1 2 ) [expCn^ - (1 + n x ) 3 

f4 ] 

4. The range is given by Feldman's J expression 
R = bE n 

where, for R in angstroms and E in keV, 

n = 1. 2/ (1 - 0.126 log Z) 
b = 250 A/pZ n/2 
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3 

A is the mean atomic weight and p the density in gm/cm . 

Numerical integration leads to a charge deposition 

profile 

N(x) = 2ir J du(l - n (u) ) f [e (x/u) ]/R' [e (x/u) ] 
o 

where u = cos 0, e is the inverse function to the range R, 

R' is the derivative of R, and f (E) is the particle flux. 

The result appears to be adequately represented by 

1 - A 

N(x) = — — exp(-x/x) 

x 

Because of the simple expressions for range and back- 
scatter, x scales with temperature as T n . Parameters are 
shown for various materials in Table 3.1. It should be stated 
that the deposition profile is not exactly represented by an 
exponential. In particular, low-Z materials (kapton) have a 
slight minimum in N (x) at x - 0. However, such behavior 
occurs only within a few hundred angstroms of the surface 
and will, in any case, be overwhelmed by charge depletion due 
to electron emission. 

While it is our belief that the average over energies 
and angles ameliorates the crudeness of the assumptions of 
this theory, several improvements might be made: Assumption 

(2) could be relaxed to take account of range straggling; 
account could be taken of increased backscatter at low energy; 
and a more accurate range expression might be used. 
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APPENDIX F.3 

ANGULAR DEPENDENCE OF PHOTOEMISSION 

It has been known for some time (Weissler and ref- 
erences therein) that there could be a strong dependence of 
photoelectric yield on the incidence angle of the incident 
light, and that, at least in some cases, this dependence could 
be well explained by the optical (or dielectric) properties of 
the material. We suppose that the yield per photon incident 
relative to normal incidence is given by 


Y(9) 

yToT 


sec 0 


E(9) 2 „ 
2 x 

E(0) 


1 + 2a(0)L 
1 + 2a(0)L 


Here, E(0) is the electric field just inside the solid cal- 
culated using the appropriate electromagnetic boundary con- 
ditions, a (8) is the E-field attenuation coefficient normal 
to the surface, and L is the escape depth for photoelectrons 
which may be estimated from low energy electron transport 
studies. ^ 

We have applied this theory to aluminum using the known 
dielectric constant and compared the results to the measure- 
ments of Samson and Cairns^ 4 ^on 94 percent aluminum alloy 
having undergone "routine polishing and deeming" (see Figure 
F.3).. Below the plasma frequency, the agreement is excellent, 
indicating little angular dependence until fairly glancing 
angles are reached. Above the plasma frequency, we predict 
a much sharper angular dependence them was measured, although 
we agree on the general form of the curve. The difference may 
be due to surface roughness and/or contamination and/or dif- 
ferences in the optical constants of the alloy from those of 
pure aluminum. 
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Yield («) 



Figure F.3. Angular dependence of photoemission compared with 
data for 94 percent A1 (Ref. 43). [NOTE: ..Data 
for pure, clean aluminum (Ref. 44) is in good 
agreement with the theoretical curve . ] 
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